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Chapter 1 

Preliminaries 



Mathematical Preliminaries 

A fundamental principle of QM is the superposition principle of dynamic states: The possible dy- 
namical states of a given quantized system may be superposed linearly and therefore may be 
represented by vectors of a certain linear space (to be defined). In this space, dynamic variable's 
are associated with linear operators. 



1.1 Vector Spaces, Deal Spaces, and Scalar Products 
1.1.1 Vector Spaces and "kets" 

Definition 1: A vector space, or linear space: 

y = {|i),|2),...»,|6),...} 

has the following properties: 

1. Is a group under addition 

• closure (\a) + \b) E V) 

• exists identity: (\a) + |0) = \a) V|a}) 

• 3 inverse: (|a) + |— a) = |0)) 

• associative: (\a) + (\b) + |c)) = (\a) + \b)) + |c)) 

2. Commutative under addition 

3. We can also define a multiplication of elements of V with num- 
bers "a, f3, 7, . . ." so that V \a) , \b) G V: 

• a | a) G V 

• a({3 \a)) = (a(3) \a) 

• (a + (5) \a) = a\a) + (3\a) 

• a(|a) + \b)) = a \a) + a \b) 

• 1 • |a) = |a) V |a) 
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Remark 1: \a) , | b) , . . . are vectors 

a, /3, . . . are scalars that define a field over which the vector space is defined 

If a, (3, . . . 6K, then we have a real vector space 

If a, j3, ... € C, then we have a complex vector space 



Remark 2: \a) is a "ket" in Dirac notation 

It is an abstraction of the more familiar column vector 

Definition 2: The vectors \i) for i = 1 . . . n are linearly independent if: 

n 

aj |i) = =^ aj = Vi 

i=l 

Definition 3: Let n be the number of linearly independent vectors in V. Then n = 
the dimension of the vector space, denoted V n , and the linearly independent 
vectors form a basis. 

1.1.2 Dual Spaces and "Bras" 

Definition 1: The adjoint of a vector \a) in the vector space defines a "bra" (a\, which 
is a member of the dual space. The dual space is antilinear . 

If \ c ) = a = \a) +/3\b), then (c| = a* (a\ + [3* (b\ 
Remark 1: "bras" are abstractions of row vectors 

Remark 2: The adjoint operator takes the "ket" from the vector space to the corre- 
sponding "bra" in dual-space and vice versa 

\ay = (a\ and (ay = \a) 

Remark 3: If there 3 any type of scalar multipliers, they are replaced by complex 
conjugates. Adjoints of complicated expressions are written in reverse order. 

Ex: [a \a) (b\ j5 \c)] j = (c| 5*7* |6) (a\ a* 

1.1.3 Scalar Products 

Definition 1: For any vectors \a) , \b) £ V n , we denote (a\b) to be the inner product, 
or scalar product. (a\b) G C, and has the properties: 
= (b\a)* 
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• (o|a) > 0, = 4^ |o) = |0) 

• Linear in ket and antilinear in bra 

Remark 1: The properties of the scalar product imply the Schwartz inequality: 

\(a\b)\ 2 < (a\a) (b\b) 



Postulate: The space of ket vectors (and the dual space of bra vectors) form a 
Hilbert Space. It is complete and separable. Here, complete refers to the 
existence of a basis in which any vector can be expressed as a superposition of 
the basis vectors. 

Definition 2: \a) ,\b) G V n are orthogonal iff (a\b) = 



Definition 3: W (a\a) = \a\ is the norm of the vector \a) 

Definition 4: The set of basis vectors, |1) , . . . , |n), that have = 5ij,Vi,j, is an 
orthonormal basis. 

Theorem (Gram-Schmidt): An orthonormal basis may be formed from a linear 
combination of basis vectors. 

Remark 2: We will write vectors as expansions of orthonormal basis so that 

n n 

\a) = ^2ai\i) and \b) = ^bj\j) 
i=i j=i 

n n 

(a\b) = a * b j(j\j) =Y.< hi 



Definition 5: A subspace is a subset of the vector space V that is also a vector space. 



Definition 6: £i and £2 & re orthogonal subspaces if every vector in £1 is orthogonal 
to every vector in £2- 

Definition 7: The vectors orthogonal to £1 form a subspace £f which is the complementary subspace 

to 6- 
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Theorem: Any vector \a) € V may be written 

| a) = \ai) + \af) 
where \a\) € £1 and |af) € £f and £i,£f C V. 

Definition 8: \a±) is the projection of \a) on the subspace of £i. 

1.1.4 Tensor Products of 2 Vector Spaces 

To represent a system described by products of wave functions in the bra-ket notation, 
we use Tensor Products . 

Example: Let ja 1 ) represent the state of particle 1. Let |a 2 ) represent the state of 
particle 2. Then the 2 particle system is represented by the tensor product: 

la 1 )®^ 2 ) = (a 1 ) \a 2 ) = |aV> 
Remark 1: An operator Cl l (Cl 2 ) only operates on [a 1 ) (|a 2 )), so that, if: 

n 1 ^ 1 ) = lb 1 ) 
n 2 \a 2 ) = \b 2 ) 

Then 

Cl 1 \a 1 a 2 ) = |&V) 
£l 2 \a 1 a 2 ) = \a 1 b 2 ) 
Furthermore, every operator Cl 1 commutes with Cl 2 . 

Ex: n 1 ^ 2 \a l a 2 ) = Cl 2 Cl l \a l a 2 ) 



1.2 Linear Operators 

1.2.1 Products of Operators 

Consider an operator Q that transforms one vector \a) into another \a'): 

n\a) = \a') 

Definition 1: a linear operator obeys the following: 
• f2a | a) = ail \ a) 
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• Cl [a I a) + \b)] = aCl \a) + /30 \b) 

• (a| af2 = (a| fia 

• [(a| a + <6| /3] A = a (a\ Cl + [5 (b\ Cl 

Remark 1: the products of 2 operators are just 2 operations carried out in sequence: 

(ClA)\a) = fi(A|a» 
In general, the product does not commute. Define the commutator : 

ClA - ACl = [Cl,A] 

Definition 2: The inverse of a operator Cl is denoted Cl -1 , and f2 _1 f2 = Of) -1 = I = 
identity operator. 

Remark 2: The inverse of a product of operators is: 

(Ciky 1 = A -1 fi -1 

1.2.2 Matrix Representation of Operators 

Consider a vector \a) and operator Q such that 

Cl\a) = \b) 

The matrix form of Cl has elements f2« = (i\ Cl \ so that for: 

| a) = \i) 



Then 



i&> = £Mi> 



bj = 0"l&> = <il^k> 



01 ^ I 
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And, in matrix and vector representation: 



Examples: 



V 




"(i|fi|i) <i|n|2) • 


•• <llnl»l>" 




ai 






<2|n|i) 






02 






.(n|0|l) 


(n\Q \n) 




_a n _ 



The identity operator I 
I 



• The Projection operator P: Consider a set of basis vectors 
|1) , 1 2) , . . . , |n). The projection of any vector onto the i basis 
vector is obtained by: 



So: 



Pj \a) = \i) {i\a) = (i\a) \i) 



Pi = \i) (i\ = projector 



Remark 1: Because the basis vectors form a complete set (in Hilbert Space), we 
have 



a,- I) 



\i) 

at 

i 



Since this is true for all la) 



1 = J» (i 



a completeness relation 



Remark 2: Projectors have the property that: 

P 2 = P 

PiPj = \i) (j| = ^|i) {j\=SijPi 



1.2. LINEAR OPERATORS 



15 



In fact, any Hermitian operator (denned below) that satisfies the above is a 
projector. 

Remark 3: The projector is the first example of an operator expressed as an outer product 
of vectors: \i) {i\. The jk element of this operator in matrix form is: 

(i\k) = Sijdik 

In matrix representation it is: 






l(i th ) 



[0,0,0,- •• ,1(1^),- ,0] 



:th 



:th 



This concept can be generalized to other operators. 



1.2.3 Hermitian and Unitary Operators 

Definition 1: Consider an operator Cl. 

• Cl = Cl* Hermitian 

• Cl = — Antihermitian 



Remark 1: Any operator Q may be decomposed into it's Hermitian and Antihermi- 
tian parts: 

(l = fo H + Cl A 



Where 



Remark 2: A hermitian operator is positive-definite: 

(a\fl\a)>0 V|a) 

Definition 2: An operator U is unitary if: 



uu^ = t)W = i (trt = u~ 1 ) 
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Remark 3: \iU\ and U2 are unitary, then so is the product of IJ-jJi. 

Theorem: The inner product of vectors | a) , \ b) is invariant under a unitary transfor- 
mation. 

U\a) = | a') 
U\b) = \b>) 

Then: 

(a'\b') = (a\&U\b) = (a\b) 

4 Transformations (Unitary) 

Consider a transformation T. 

Definition 1: An active transformation transforms all vectors in a subspace: 

\a) —> a') = T | a) 

Definition 2: A passive transformation transforms all operators on that space, leav- 
ing the vectors unchanged. 

Ci -> n' = fnf- 1 

Remark 1: In general, T does not conseve Hermitian conjugation. For this, we 
require that if: 

Q' = fnf' 1 and (ftt)' = f&f- 1 

Then: 

(0 / ) t = (fnf' 1 )^ = (f- x )Wf* must = (ffitf- 1 ) 

^ f f = f" 1 or ff' 1 = f- x f = I 
=^ T is unitary 

Remark 2: If we have unitary transformations, active and passive transformations 
are equivalent: 

(a\ O 1 6) -» (a'| A |&'} = (a| C/^C/ |b) where |a'> = 17 |o) 
Remark 3: In addition of Hermitian conjugation, unitary transformations conserve: 
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• Trace 

• Determinant 

• Algebraic equations between matrices and vectors 

Remark 4- Any Hermitian matrix may be diagonalized by a unitary transformation. 

UflW = &d where fop is diagonal. 

The diagonal elements are called eigenvalues of Cl. They are real, and they are 
solutions of the secular or characteristic equation. 

det(0 - A J) = (See below) 

Remark 5: A single unitary matrix may diagonalize 2 (or more) Hermitian matrices 
O, A if: [Cl, A] = 0. 

Remark 6: If the vector space is real, then 

C/ _1 = — > U~ l = U T =>• orthogonal transformation 
which is a rotation . 

1.2.5 The Eigenvalue Problem 

Definition 1: Consider a linear operator A and \u). Then a £ C is an eigenvalue of 
A if: 

A \u) = a \u) 
and \u) is an eigenvector or eigenket. 

Remark 1: If \u) is an eigenket of A, then any multiple of \u) is also an eigenket be- 
longing to the eigenvalue a. These eigenkets form a subspace of the eigenvalue a. 
If n equals the dimension of the subspace, then: 

• n = 1 => single or non-degenerate eigenvalues 

• n > 1 degenerate eigenvalues, where n equals the order of 
degeneracy 

Remark 2: If {v\ A = a' (v\, then (v\ is an eigenbra of A. If A is Hermitian, then: 

• The eigenvalues for bras and kets are the same 

• The eigenvalues are real 
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• The subspace of eigenbra's is dual to the subspace of eigenkets 
for the same eigenvalue. 

Remark 3: Eigenkets make up the columns of the unitary matrix that diagonalizes 
a Hermitian matrix. 

Remark 4 : In some cases, the unitary matrix is time dependent, and gives the 
time evolution of a vector in some initial state. 



x(t)) = U(t) |x(0)> 



In these cases, it is called the time-evolution operator, or propagator. 



.6 Functions of Operators, Generalizations to Infinite Dimensions, 
and More 

In some cases, we will encounter functions of operators. (In fact, the propagator can 
be expressed as e lHt / h with H the Hamiltonian, see Section 2) 

In order to evaluate the action of these on vectors, we expand functions in a Taylor 
Series. 

oo 
n=0 



Remark 1: Expressing /(0) as a Taylor Series in fine so long as the series converges. 
Remark 2: It is often helpful to transform into an eigenbasis. For example, if 



W) = e 
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Then 



3 O" ^ 1 



n! ' — ' n! 

n=0 n=0 



Ei £ ^ », 

A? 



^ n! 

n=0 



V 



A-2 



e Xn J 



V 



Remark 3: A power series expansion is also useful for taking derivatives with respect 
to a parameter. For instance, consider: 



Then 



da 



d ^a n Cl r 



E 



da — ' nl 

n=0 



E 



n=0 



E 



nl 



1(^1)! 



E 

Oe 1 



a- 



n=0 



nl 



Remark 4 : Here ft behaves as if it was a complex number. The important distinction 
is that multiple operators in general do not commute. 



Generalization to Infinite Vector Space 

Most of the the formalism we've encountered deals with finite vector spaces. 
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Sometimes we will need to consider infinite vector spaces which contain vectors having 
infinite norm. The eigenvalues in this case are part of a continuous spectrum, and 
most of the formalism still holds. (For instance, Section 1.2.5, remark 2 remains 
valid.) Other generalizations are: 

(x\y) = 5{x — y) — > Delta Dirac Function 

where 5(x — y) = x / y 

rb 

I 5(x — y)dy = l Va < x < b 

J a 

Also, 

I = / \x) (x\ dx 

J a 

expresses completeness for a continuous spectrum. 

Remark 5: One caveat concerning operators in a continuous vector space: 

A f = h 

is not sufficient to be Hermitian. 



Example: Consider the operator: 

^ d 
p = —in— 

dx 

acting in real space representation. The matrix elements in this representation: 

(x\p\y) = -ihS(x - y)^j- 
where x,y are eigenvectors. Checking to see if p is Hermitian: 

is (9\p\i>) = &\p\e)* 

(e\p[if;) = [ [ (e\x)(k\p\y) (yty) dxdy 



a J a 
b r b 



a J a 
b 



(9\x) (-ih)5(x - y)^- (y\ip) dxdy 
ay 

x) (—ih)-^- (x\ip) dx 
dx 



b d 
(x\9)* {-in)— {x\rf>) dx 



dx 

(-ih) (x\9)* (x\if)) \ b a - I (-ih) (x\il>) ^- (x\9)* dx 



b d 



(-ih) (x\9)* (xty) \ b a - I l(^\x)(-ih^)(x\9)) dx 



dx 
d 

dx' 
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which will equal what we want if: 

i-ih) (x\ey (x^)\ b a = o 

So we require surface terms to vanish for continuous Hermitian operators. 

Remark 6: Restricting surface terms to vanish is equivalent to considering the 
Physical Hilbert Space. 

Remark 7: If (a, b) = (—00,00), then the typical functions 

(a#> = i>{x) ^ x ^ a ,b or e ikx 



End of Chapter Asides: 



gifcXg— ik'x 



00 



-00 



if k = k' 



If k + k>, then: 

1 rL+A 

hm e ik *e- ik ' x = lim - / e^ k ~ k > dx = 



A— >oo 

Some other properties: 



p = hk 

x\x) = x\x) or k \k) = k\k) 

t . d 
=^ k = —1 — in position 

dx 

t . d 

=>• k = 1— m momentum 

dx 

— / e ik( - x - y 1 dk = 5(x - y) 
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Chapter 2 

Classical Mechanics Review 



2.1 Lagrange's and Hamilton's Formulations of Classical Theory 

2.1.1 Lagrangian 

Consider a system with N degrees of freedom. In classical theory, the dynamical state 
of the system is completely specified by the generalized coordinates qi,q2, ■ ■ ■ ,Qn an d 
generalized velocities q\ , q\ , . . . , q~N ■ 

Remark 1: The description gives the evolution of the system in configuration space. 

Definition 1; The Lagrangian of the system: 

C = C(q 1 ,q 2 ,.. ■ ,qN,q\,Q2, ■■■ ,9iv) 

is a characteristic function from which the equations of motion can be derived: 

d_ /d£\ _d£ . _ 

dt \dqi J dqi 

These are called Lagrange's Equations of Motion. 
Remark 2: 

dC 

oqi 

These are called Lagrange conjugate momentum. 



Remark 3: If all forces in the system are derivable from a static potential (the 
situation we usually encounter), then the Lagrangian takes the form: 

C= T - V 

Kinetic Potential 
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Remark 4 : The Lagrange Equations of Motion may be derived by minimizing the 
action of the system: 

r>t2 

Action = S = I Cdt 



5S = 5 I Cdt — ► Principle of Least Action 



.2 Hamiltonian 

An alternative but equivalent formulation of mechanics is in terms of the Hamiltonian. 
Here, the dynamical state of the system is represented by a point in 2n dimensional 
phase space, whose coordinates are the N position coordinates, q±, q 2 , ■ ■ ■ , qN, and the 
N corresponding conjugate momenta pi,p 2 , ■ ■ ■ ,PN- 

Definition 1; The Hamiltonian is defined as the function H: 

N 

dq, 



dC 

H = H(q 1 ,q 2 , . . . ,qN,Pi,P2, ■ ■ ■ ,Pn) =X^* t T" ~ jC 



The equations of motion are: 

dH 



,i = l...N 

OPi 

Pi = — t = l...N 
oqi 



Remark 1: If C = T — V, then H = T + V = the total energy. In more general 
situations, we still consider the Hamiltonian to be the total energy of the 
system. 

Remark 2: 

jj_dH _dH 
= lit ~ ~dt 

This can be seen by using the chain rule and the Hamiltonian Equations of 
Motion. This implies that if the Hamiltonian is not an explicit function of 
time, then energy is conserved and we say that the system is conservative. 



Remark 3: The Hamiltonian formulation is particularly important for QM since the 
generalized coordinates and the conjugate momenta become operators in the 
HUbert Space. 



2.1. LAGRANGE'S AND HAMILTON'S FORMULATIONS OF CLASSICAL THEORY 
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2.1.3 Cyclic Coordinates, Poisson Brackets, and Canonical Transforma- 
tions 

Definition 1: A coordinate qi is cyclic if: 

dH 

Pi = = 
dqi 

which is a statement of conservation of conjugate momentum. 



Definition 2: The Poisson Bracket between two variable uo(p, q),X(p, q) is given by 

^ / duo dX duo dX \ 
~ ^ V% dpi dpi dpi J 



Remark 1: Consider u (p, q). Then 



duo / duo duo 



dt \dpi dqi 



/ duo dH duo dH 



E 



K dpi dpi dqi dqi 
= {uo,H} 

So, if a variable's Poisson Bracket with the Hamiltonian vanishes, then that 
variable is conserved. 

Remark 2: Note the similarity between the Poisson bracket and the commutator. 

Remark 3: 

Ui,Qj} = {Pi,Pj} = 
{QhPj} = Sij 

Definition 3: A transformation of coordinates and momenta 

q -> q'(q,p) 
p^p'(q,p) 

that preserves Hamilton's equations of motion is canonical . 

., dH , .. dH 

dp. dq[ 
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Remark 4 : Poisson brackets are invariant under canonical transformations. 
.1.4 Symmetries 

Theorem (Nother's Theorem): Symmetry properties in the Hamiltonian imply 
the existence of conserved quantities (conservation laws). 
Proof: See Goldstein (p587, Section 13.7) 

Remark 1: In classical mechanics, this is equivalent to finding a variable whose 
Poisson bracket with the Hamiltonian vanishes. In QM, we'll see that an 
operator whose commutator with the Hamiltonian vanishes, then this implies 
that there exists a conserved quantity. 

Examples: 

• Translational Invariance =>• Conservation of Linear Momentum 

• Rotational Invariance => Conservation of Angular Momentum 

• Time Translational Invariance =>- Conservation of Energy 
Some examples closely related to QM: 

• Reflection Invariance => Conservation of Parity 

• Invariance under Permutation =^ Conservation of Symmetry 
(symmetric, anti-symmetric) 

Remark 2: Symmetries imply degeneracy in the eigenspectrum. Breaking the sym- 
metry often lifts the degeneracy of the system. (Eg: Zeeman Effect) 



Chapter 3 

Postulates of QM 



3.1 Postulates 



I. The state of a particle is represented by a vector \tp(t)) in Hilbert Space. 

II. Every observable in classical mechanics corresponds to a linear Hermitian 
operator in QM. In particular: 

x -> X 
P 



T=P 2 /2m 
V=V{X) 



The particular action of the operators depends on the choice of basis. 

III. Any measurement of the observable associated with the operator f2 will 
result only in values, uj, which satisfy 

IV. For a system described by the normalized wave function \ip), the average 
value of the observable corresponding to ^ is 

This is called the expectation of £1. Note: if is not normalized, one need 
only divide the above equation by (tp\ip)- 

V. The time evolution of the state vector is governed by the Schrodinger Equa- 
tion: 

ihf t \m) = Hm)) 

where H is the Hamiltonian. 
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Remark 1: In measurements, \ip) may not initially be an eigenstate of Cl (Postulate 
III). But one can expand \ip) in an eigenbasis: 

i 

The coefficients in the expansion {oJi\ip) give the probability of finding the 
system in eigenstate |o»j) according to: 

P(u) = | <u#) | 2 



Remark 2: Once an eigenvalue is measured, the system has probability of 1 of being 
in the corresponding eigenstate. The measurement (even if "idea" ) has changed 
the state of the system; this is known as the collapse of the wave function. 

Example: 

in 1 I » 1 , . 1 . . 

m = 2 \ u v + 2 + ^/f l^ 3 ^ 

( 1/4 of being in state w\ 

The system has probability < 1/4 of being in state u 2 

(1/2 of being in state C03 

Say that a measurement gives the eigenvalues of then the state after the 
measurement is 

W) = 1^3) 



Remark 3: Because coefficients correspond to probabilities, expansions of normalized 
wave functions must have 

i 



Remark 4 : We must be careful with the consequences of Postulate II! Operators, 
unlike observables, must obey commutation relations. In general, we will take 
the symmetrized form (though this is not a universal prescription). 

XP + PX 

Eg: xp -> 

This also makes a Hermitian operator! 



Remark 5: For degenerate eigenvalues, probabilities represent projection onto the 
corresponding "eigenspace" . Probability is for measuring a certain eigenvalue. 



3.2. COMPATABILITY, INCOMPATIBILITY AND UNCERTAINTY 
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Remark 6: For continuous eigenspectra, we require that 



/ 



P(u) dco 




And P(lo) is called the probability density. P(u>) doj is the probability that 
the measurement will give a value in the range uj to to + alio. 

Remark 7: Not all QM observables have classical counterparts. (Eg. Spin) 

Remark 8: The probabilistic nature of QM is a fundamental aspect of the theory. 
It is not the same as the classical probabilistic description of a many-particle 
system where the probabilistic approach is for practical purposes and where, 
in principle, know the exact q(0),p(0) and the nature of the interaction would 
give exact time evolution solutions. 



3.2 Compatability, Incompatibility and Uncertainty 
3.2.1 Compatible Variables 



Definition 1; Two variables are compatible if their corresponding operators com- 



Remark 1: If two (Hermitian) operators commute, one can find a unitary transfor- 
mation which simultaneously diagonalizes them. Also, they will have a simul- 
taneous eigenbasis. (See, for example, HW 2, Problem 2c, Shankar 1.8.10) 

Definition 2: Two variables are incompatible if their corresponding operators do not 



Remark 2: Incompatible variables do have have any simultaneous eigenkets. 



mute: 



[ft, A] = 



commute. 



[ft,A]/0 



Remark 3: Two variables do not have to be either compatible or incompatible; they 
may have some simultaneous eigenkets, but still do not commute. Thus they 
fall into neither classification. 
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Definition 3: For any number of variables that mutually commute, one can find a si- 
multaneous eigenbasis which simultaneously diagonalizes each observable. The 
set of these commuting variables is called a complete set of commuting observables 
or a complete set of compatible variables if the eigenbasis is unique (up to an 
overall phase); ie: is degenerate. 

Remark 4 : Consider an observable f2 and state vector \tp). One can expand in 
an eigenbasis of fl. The basis of f2 is unique (up to an overall phase) if none of 
it's eigenvalues are degenerate. If there is a degenerate eigenvalue, and if there 
is an observable A that is compatible with fl, then the simultaneous eigenbasis 
will be unique (up to an overall phase). If there are no longer any degenerate 
eigenvalues, then f2 and A form a complete set of commuting observables. 

Remark 5: Compatible variables may be measured simultaneously and precisely. A 
simultaneous measurement of all compatible variables will completely define 
the state vector of the physical system (and furthermore, that state will be 
unique). Also, the dynamical state of the system is completely specified by 
the quantum numbers associated with the eigenbasis. 

Remark 6: A simultaneous measurement of a complete set of compatible variables 
is akin to preparing the system at the time to, after which the system evolves 
according to the Schrodinger Equation (See 3.1, Postulate V) 

Definition 4: A pure state is one that has been prepared by measuring all values of 
compatible variables so that the state vector is know exactly. 

Remark 7: In practice, a complete preparation is rare, and the dynamical state is not 
specified exactly. In this case, the system is in a mixed state , and is described 
by a statistical measure of states. 

Remark 8: Mixed states include the statistical description in a classical sense, as 
well as the inherent probabilistic nature of QM. We will visit this shortly in 
Section 3.3 below). 

.2.2 Incompatible Variables and Uncertainty Relationships 

Recall Section 3.1, Postulate IV ((A) = U \tp)). 
Definition 1: The uncertainty in f2 is given by: 




1/2 
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In a statistical description, this is also called the standard deviation or root- 
mean-square, and is a measure of the uncertainty due to fluctuations. 

Note: 

1 /2 

An = ((£i - (n\) 2 \ 




Remark 1: If the state \ip) is an eigenstate of ft, then there is no uncertainty. 
(AO = 0) 

Uncertainty Derivation: Consider two incompatible Hermitian operators f2 and A 
that have 

[fi, A] = if 

where V is also Hermitian. Let [Cl, A]+ = j)A + Af2 denote the anticommutator, 
and let lib) be a normalized state so that: 



A 



A| 



Then 



(AO) 2 (AA) 2 = (^|(fi- (0 



(A — (A 



If we let 5ti = Cl — (Cl) and 5 A = A — ( A ) , then we have: 



(AO) 2 (AA) 2 = (V>|(<m) 2 |V>> <Vl(5A) 2 |V> 

= ^\(5Cl)(6Cl)\^) (^\(5A)(5A)\^) 

= (tl>\ (5Clf(5Cl) \if>) {i>\ (<5A) t (<5A) \i>) since Hermitian. 
5Cl^ 5 A | ijj) by Schwartz Inequality 

5m AW 



> 
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But now, we also have that 



, rAwr . 5Q5A + 6A6Q 5ttd A - 5A5Q 

(sn)(SA) = + 

= ±[8Cl,6k] + + ±[6Cl,6k] 



And so we get that: 

(Af2) 2 (AA) 2 > 



snsA 
i 



[sn,SA]- 



+ 2^|r| 



This is the most general uncertainty relation between any two Hermitian op- 
erators. 



Remark 2: If £1 and A are canonically conjugate, then 



r = M 



And thus 



(Af2) (AA) > - (tp\ [5n,5A]. 



h 2 

+ T 



> 



AQAA > - 
~ 2 



THE uncertainty relation 



Remark 3: The equality holds for 

• <5fi |V>) = (constant) • 5A 

. (ij;\[Sn,6A] + \^=0 

Remark 4 : The most popular statements are 

• AxAp > | <— holds for each component 

• AEAt > § 



3.3. PURE STATES AND MIXTURES: DENSITY MATRICES 
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3.3 Pure States and Mixtures: Density Matrices 

Recall from Section 3.2.1 that: 

• A pure state is one where we know the dynamical state of the system exactly 

• A mixed state is one where the exact state of the system is unknown. 

Suppose that we want to measure the quantity O in a mixed system. Then there is a probability 
P m of obtaining the average ( Cl\ = (m\ Cl \ m), and 

\ / m 

(^) = J2 P m( m \&\ m ) 

m 

where Vt indicates an average in the classical sense. 

Remark 1: In the above equation, the |m)'s are normalized but not necessarily orthogonal, and 

P m > and Yl P ™ = 1 

m 

Remark 2: There are two kinds of averaging in the equation above: 

• Quantum average {m\ ft \m) for each system in \m) 

• Classical average over systems with different \m) 

3.3.1 Density Matrix 

Definition 1; A density matrix (or operator) is a convenient way to describe a mixed 
state: 

9 = / .Pm \ m ) ( m \ 
m 

This is called the density operator. 

Remark 1: The average value of in a mixed state is given by 

(Clj = Tr (pCl) 

To see this: 

Tr (pCl) = ( n \ 9& \ n ) 

n 

= Pm (n,\m) (m\ O \n) 

n,m 

= P m (m\ Cl \m) 

m 
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Remark 2: It is sufficient to know p in order to calculate all measurable physical 
quantities, and thus every quantized, statistical mixture is exactly and com- 
pletely defined by it's density operator. 

Remark 3: p = is a special case: it is the density matrix formulation of a 

pure state (all P m = except 1). In this case: 



= Tr(V>|0|V>) 
= Tr (|V) (V| Cl) 
= Tr (p£l) 

which is the usual quantum average. 

Properties of the Density Matrix: 

1. p is positive definite: p\ip) > Vip G 



(vi p iv) = Pm ^\ m ) h^) = Yl p -i M m > i 2 ^ 

m m 

2. p is Hermitian 

P ] = \ ^2 P m \m) (m| = Pm = ^ 

\ m / m 



3. p has unit trace 



Tr (p) = Tr (^P m |m)(m|) = ^P m Tr (|m) (m|) = ^P m = 1 

\ m / m * m 

4. Tr (p 2 ) < 1 

Tr (p 2 ) = P m Pn |m) (m|n) (m| = P 2 . < 1 



3.3.2 Density Matrix of a Pure State 

If the density operator represents a pure state, ie. p = \m) (m\, then we also have 

• p 2 = p — ► (projector) 

• Tr (p 2 ) = 1 — > (see properties above) 



PURE STATES AND MIXTURES: DENSITY MATRICES 
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Suppose the density operator is expressed in its matrix form in a particular represen- 
tation: 

/ pu P12 ■■■ Pln\ 
P21 ' ' ■ '■ 



\Pnl PnnJ 



with p nm = (n\p\m). 



Question: Given the explicit matrix form, how do we determine whether it represents 
a pure state or a mixed state? 



Answer: Since the pure state can be written p 



Then we have that: 



a, 



i- 



a, 



\a±\ a±a 2 
a n a\ 



and 



a\a r 



Now, if p satisfies a pure state, then its elements must satisfy: 

Pijpji = pijPij = \Pij\ 2 = Pupjj 

If this holds for all i,j, then the density matrix is a pure state; otherwise it is a mixed 
state. This can be seen by doing the matrix multiplication. 



3.3.3 Density Matrix of Mixed States 

Given that p = ^ • Pj < — it is not required that the ^-'s be orthogonal. 

Example: Consider a two state system where |1) and |2) form the orthogonal basis 
for the 2D Hilbert space. One possible density operator is: 

=^ 50% mixture of state |1) and nd-*-) 1^)) wn i cn are n °t orthogonal. 
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Remark 1: The density operator in the above is a non-unique representation of 
non-orthogonal mixtures. 

Remark 2: It is always possible to write p as a mixture of orthogonal states. For ex- 
ample, consider the general density operator p = Pj Let R be the 
matrix representation of p in the orthonormal basis defined by 1 1) , |2) , . . . , \N). 
Then: 



Pn P12 • • • Pin 

P21 



PNl 



Pnn 



Pmn = (m\ p \n) = } Pj {rn\ipj) (tpj 



n) 



Then we write: 



p=[|l),|2),... 



(N\ 



ID, 12 



J / J * ' ' J 



[\4i),\<h),-,\4>N) 



\N)]&- URU^ 
Pi 



PN. 



u 



(If 

(N[ 



'>N 



P ~ ' y A 3 

J 



Remark 3: If |1) , |2) , . . . , \N) form a orthonormal basis, then so does \4>i) , {fo) , • • • , 
because U is unitary. 
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Recall Example: 



# = 5l 1 >< 1 l + 5{^ (|1> + |a>) ^ ((1| + (a|) 

= 1 11> <l| + j|l) <2| + i|2)(l| + i|2) (2| 



"3/4 1/4" 




"(If 


1/4 1/4 




.(2|_ 



(II 

(2| 



= [|1>,|2>] 
= [|1>,|2>]H 

Now we would just diagonalize R to create a diagonalized p in some basis. 
Remark 4 : In some cases, the non- uniqueness of p goes even further. 
/5=[|1),|2)] 
= [|l),|2)]f/[/t 



"1/2 " 




"(If 


1/2 




.(21 



1/2 
1/2 



(1| 
(2| 



\1 VC7 



1/2 
1/2 



(1| 
(2 1 



[|1>,|2>]J7 



[|01> , |^2>] 



In this case, p can be expressed in an infinite number of forms depending on 
the unitary transformation U used. 



"1/2 " 




r 


l/2_ 







3.3.4 Time Evolution of the Density Matrix 

Regression: QM can be expressed in equivalent "pictures" depending on whether 
the transformations which govern time evolution are passive or active. (Recall 
§1.2.4) 

• The Schrodinger picture (active) has state vectors evolve with 
time while operators remain constant. 

• The Heisenberg picture (passive) has states remain constant 
while operators evolve. 

In addition, there is 

• Interaction Picture : Here, the non-interacting part of the Hamil- 
tonian is in the Schrodinger picture while the interacting piece 
has time-dependent operators. 
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Time Evolution of p: Recall §3.1 - Postulate V -> ih\ip(t)) = H What then 

is the time evolution of p? 



d 



dt 



l^>l^<^l) + l S l^>)<V'i 



d 



3 

= Y. P > 

j 



in 



Remark 1: This is the quantum version of the Liouville Equation from classical 
statistical mechanics, (see §2.1.3, remark 1 with p = the density of phase 
space.) 

Remark 2: We'll see later that the Heisenberg equation of motion for an operator, 
Ah, is given by: 

A H = hA H ,H] 

in 

Witness that the equation of motion for p has the commutator in reverse order. 

Remark 3: The density matrix is a special operator since it changes in time in the 
Schrodinger picture. 



3.4 Schrodinger Equation 

Recall §3.1, Postulate V, which said that the Schrodinger equation is: 

ihj t \m) = H\m) 

where H is the Hamiltonian of the system. 

Remark 1 : In general, one may use Postulate II to construct the Hamiltonian H from the classical 
Hamiltonian. Modifications to this are required if there is no classical analog (eg spin). 

Remark 2: Though we can only postulate about the form of the time evolution of a vector in 
Hilbert space, we do know that it must be linear and homogeneous (so that we can have 
linear superposition). We also know that it must be a DE of the 1st order with respect to 
time (so that the evolution is uniquely determined by the state at some initial time). 



3.4. SCHRODINGER EQUATION 
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Remark 3: Predictions of Quantum Theory must coincide with classical mechanics in the regime 
where classical mechanics is valid (so there is some formal analogy with classical mechanics). 
This is called the correspondence principle. 

3.4.1 Solving the Schrodinger Equation 

For now, we will consider H = H{x,p), where the Hamiltonian is a function of po- 
sition and momentum, but not explicitly time. In general, the time dependence (if 
present) is treated perturbatively, so we will discuss Hamiltonians with explicit time 
dependence later. 

Solving the equation formally gives a solution of the form: 
Remark 1: One can write this as: 

\m) = u(t)\m) 

where U(t) = e ~ lHt / h [ s called the propagator or time evolution operator. 
Remark 2: Since H is Hermitian, U is unitary. 

Remark 3: Time evolution may be thought of as a rotation of a state in Hilbert 
space. 

To evaluate (*), it is useful to expand |^(0)) in an energy eigenbasis. This will diago- 
nalize the Hamiltonian and allows one to write: 

-iHt/H 




To determine the energy eigenbasis: 

1. First solve the time- independent Schrodinger equation: 

H\E) = E\E) 

2. The expand the wavefunction in terms of the energy eigenbasis: 

a E {t) 
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3. Now substitute this into ih\i])) = H\ip) to obtain: 

d 

in— a E {t) = Ea E 
=► a E {t) = a E (0)e- iEt / h = (E\i/;(0)) e~ lEt l h 

=> \m) = E \ E ) w(°)> e ~ iEt/h = \m) 

E 

This, by comparison, gives a propagator of the form: 

(j(t) = E \E) {E\ e~ iEt/h 

E 

Remark 4 : If the eigenspectrum is degenerate, an additional index is required to 
specify the state: 

= EE l^«> (E,a\e- iEt / h 

a E 

If the eigenspace is continuous, then ^2e ~^ I ^E 

Remark 5: \E(t)) = \E) e. lEt / h are the normal modes called stationary states since 
the probability distribution for any variable is time independent in these states. 



Chapter 4 

Simple Problems in 1 Dimension 



In general, we will need to chose a basis to solve Schrodinger's equation. Though one can choose 
any basis, some choices are better than others depending on the form of the Hamiltonian. The 
following are some examples of 1 dimensional problems. 



Remark 1: The x basis (position basis) is the basis for the alternative form for QM based on the 

id., and H^ihg. 



wave equation. In this basis, x — ► x, P — > —ih4^, and H — ► ihf 



4.1 Free Particle 

4.1.1 Solving the Schrodinger Equation 

In the free particle case, we have 

1 ~ 9 

H = P 2 

2m 

Solving the time dependent Schrodinger equation: 

which gives us general solutions of 

= \E) e- iEt/h 

Substituting this back into (*) gets us the time independent Schrodinger equation 

p2 

— E) =E\E) 
41 
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Remark 1: Here, we may choose the momentum basis since \P) is also an eigenstate 
of P 2 

p2 

IP) = E IP) 

2m 1 1 17 

2m 



=4> p = ±V2mE 
Thus the two orthogonal eigenstates are: 

|P+> 



V2mE 



\E- 



P 



-V2mE 



Remark 2: The energy eigenstate is degenerate since a particle can be moving to 
the right or the left with the same energy. Using the compatible operator P, 
we can determine the possible outcomes so that a measurement of p, and thus 
E, will give us a complete description of the dynamical state of the system. 

Remark 3: The most general state is given by 

\i/>) = a\E + )+P\E„) 
which is a superposition of the particle moving to the left and right. 

Remark 4 : Note that measuring the momentum gives us no information about the 
position of the particle since Ap = Ax = oo by Heisenberg's uncertainty 
principle. 

.1.2 Propagator for a Free Particle in the Position Basis 

The matrix elements for the propagator in the x-basis are 

/oo 
dp{x\e-W 2mh \P) (P\y) 
-00 

(x ) 

Now consider: 



1 f 

( x \y) = 5(x-y) = — / dke ik ^ 

/oo 
dp(x\P) (P\y) inserting a complete set 
-oo 
poo 

= H / dk (x\k = p/h) (k = p/h\y) 



oo 



PARTICLE IN A BOX 
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where the last step follows since we know p = Hk. Now, by comparison, we can see 
that 

(k\y) = ^ lk{x - y) 

or 

(x\k = p/H) = —t= e ipx / h = ( x |p) 

Now, using this to evaluate U(x,t;y), we get 
1 r 00 

U(x,t;y) = — I e *p(*-»)/R e -* a */ 2mR dp 

-it ( m(x — y) x 





poo 




'— oo 




POO 




/ exp 


1 


2mh-K 



from completing the square 

im(x— y) 2 /2ht 



2irh V it 



m c im(x-y) 2 /2M 



2iirKt 



Remark 1: For a free particle whose state is known at some time to, the state at 
some t > to is given by 



ip(x,t) = J U(x,t-y,t )ifj(y,t )dy 

where 

U{x,t;y,t ) = (x\U(t-t )\y) 

Note that this depends only on t — to- This is because there is no explicit time 
dependence in the Hamiltonian. 

Remark 2: If we consider a free particle whose position is known at time to: 

ip{t Q ) = 5(x - xq) 

Then 



. 2irhit 

Physically, this corresponds to the delta function broadening as a function of 
time. 
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-L/2 




L/2 



Figure 4.1: Diagram of the "Box" Potential 

4.2 Particle in a Box 

Consider a particle under the influence of a potential V(x) where: 



V(x) = { 



\x\<L/2 
oo otherwise 



We will choose the position basis for the time-independent Schrodinger equation: 

S2 



« = k +v H 



E\ 



(x\H\ 



(.r 



1 



-ih- 



d 



+ V 



E(x\i)) = Ei/>(x 



2m dx 2 



2m V dx 
Combining the last two, we have 

d 2 2m 

Now in Regions I and III, ipi = ipni = (we'll see why in §4.4) 
In Region II: 



ip(x) + Vijj{x) 



V = 



dx 2 



2mE 



ipu = Aex.p i\j2mE /h 2 



h 2 



+ B exp 



l v / 2mE/h 2 c 



Ae tkx + Be %kx where k 



2mE 



n 2 



Remark 1: We require that the wavefunction must be continuous 

= i>n{-L/2) and ipn(L/2) = ^ m {L/2) 



4.2. PARTICLE IN A BOX 
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So then, applying the boundary conditions, we have that 

= Ae~ ikL/2 + Be ikL/2 
= Ae ikL/2 + Be~ ikL/2 

Which implies 

/ c -ikL/2 JkL/2 \ / A \ 

|0) 



e -ikL/2 e ikL/2 \ / A 
e ikL/2 e -ikL/2 J i B 

The non-trivial solutions to this problem are found by finding the zeroes of the determinant. 

e~ ikL -e ikL = = -2ism(kL) ^k = — n = 0, ±1, ±2, • • • 

L 

Now we substitute k back into (*) to get the allowed values of A and B: 

Ae -inn/2 + Be inn/2 = q 

_l /< R c * n7r R / 1 if n is even 

^A — -&e — — a ■ if n is odd 

Thus we have that 



A{e ikx +e~ lkx ) n odd 

^4^ e ifea;_ e — !fea;^ n even 



Remark 2: Choose A so that J ip*ipdx = 1 (normalize): 



2 cos (2ipi) n is odd 

j- sin^SJ 2 ) n is even 



Remark 3: If n = 0, then ?/> = everywhere, which is not very exciting. Also, ± values of n are 
physically equivalent, so we only deal with positive n to make life simpler. 

Remark 4- Comparing our expressions for k: 



2mE nir 



H 2 L 

We find that the energy in the box in quantized 

n n 2 7T 2 h 2 

En = ^^ n = l,2,3,-.. 
Remark 5: As L goes to infinity, the energy spectrum becomes continuous. 

Remark 6: The ground state, n = 1, has E = tt 2 H 2 /2mL 2 > which in a necessary consequence 
of the uncertainty principle. 

Ax = L/2 Ap > h/2 
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since 



= Y (p 2 ) — (p) 2 and (p) = 

>(2™(^}) 1/2 = j 

L irh h 

=> AxAp = > - 

2 L 2 



4.3 Step Potentials: Reflection and Transmission 

Consider the potential: 



x = 



Here, the time-independent Schrodinger equation is 

h 2 d 2 



2m dx 2 
2m dx 2 



ij) = (E- Vi)tp , x < 
ifj = (E- V 2 )^ , x>0 



-v 2 



(4.1) 



Remark 1: We must consider two different energy regions: V\ < E < Vi and E > V2 



Vi < E < V 2 



We expect the wavefunctions to have the following form: 



e ±i v /2m(E-V 1 ) [ X<Q 
e ±^/2m(V 2 -E) X>Q 



(4.2) 



Remark 2: The actual solution is constructed as linear combinations of the forms in 
(4.2) 



4.3. STEP POTENTIALS: REFLECTION AND TRANSMISSION 
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Remark 3: The solution oscillates in the region where E < V2, and decays (or grows) 
exponentially when E > V<j. 

Remark 4- To satisfy the boundary conditions, we impose the following conditions: 

1. ip{x) must remain finite as x — > 00 to be physically realizable 

2. ip(x) and ijj'(x) must be continuous at the discontinuity at x = 0. 
This gives a continuous probability in space (no source or sink 
of particles). 

Taking these considerations into account, the general solution will have the following 
form: 

, f A + e ikx + A_e~ ikx , x<0 
Mx) = { 4.3) 
^ V ' \ B + e KX + B_e~ KX , x>0 V J 

where 



hk = v / 2m(E - Vt) ; Hk = ^2m{V 2 - E) (4.4) 

Now, from Remark 4, we know that: 

• B + = (from condition 1) 

• A + — A- = B_ (from condition 2, continuity of ip) 

• ik{A + — A_) = —kB_ (from condition 2, continuity of ip') 
Combining these and solving: 

ik(A + - A_) = -k(A + + A_) 
(ik + n)A + = (ik - k)A^ 



where 



(k 2 + k 2 )e-^ = -(k 2 + fc 2 )e^A_ (4.5) 
k 

tan(0) = 

K 

^A_ = -e~ 2i<i> A+ 

Remark 5: A + and A_ are not necessarily complex conjugates of one another, but 
we can introduce an overall phase. We multiply (4.5) by i so that 

A + = -e l<t> and A_ = --e"^ 

where A is some real constant. 
Then, we have 

J Asin(ki + 4>), x<0 
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Remark 6: In the region x < 0, we have a wave traveling toward the step potential 
and a reflected component of equal magnitude traveling away. We see this by 
the time-dependent Schrodinger equation: 

sm{kx + <f>) 

_ ^_ e -i(Et/h-kx-<j)) _^_ e -i{Et/h+kx-<t>) 

2i 2i 



Wave going right Wave going left phase shifted by tt 



Remark 7: Note that the wavefunction doesn't vanish for x > 0; there is some 
barrier penetration that decays more rapidly for larger step sizes. If V2 — > 00, 
there is no barrier penetration. This is the justification for setting ip(x) = 
outside of the infinite square well in the previous section. 



Remark 8: In wave phenomena, this penetration is a familiar part of reflection and 
transmission. In E&M, the penetrating part of light at an interface is called 
an evanescent wave. 



E > V 2 : 



In this case, the solution of the Schrodinger equation is: 

A +e ikx + A-e~ ikx , x<0 



^( x ) = S „ „ _^ „ ( 4 - 7 ) 



B + e ikx + B_e~ ikx , x>0 



with 



hk' = yj2m(E - V 2 ) (4.8) 
Imposing boundary conditions, we have 

A+ + A_ = B + + B_ (4.9) 
ik(A+ - A_) =ik'(B + - BJ) (4.10) 

Remark 9: Note that we only have 2 equations to specify 4 unknowns. This is 
because there are 2 linearly independent solutions for each energy (energy is 
2-fold degenerate, or doubly degenerate, or has a degeneracy of order 2. 



The 2 possible scenarios corresponding to this degeneracy are: 

So we see that the wave is "launched" from either the right or the left and is partially 
reflected (or transmitted) at the step boundary. 
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Vi 



x = 
Case A 



x = 
Case B 



Consider now scenario A above, where the wave is launched from the left. In this case, 
we have 2?_ = 0, which, from Eq (4.9) above implies that 

B + = A+ + A_ 

Using Eq (4.10) above, we also have that 

ik(A+ + A_) = ik'(A+ + A-) 

A + (k - k') = A^(k + k') (4.11) 
Thus we get that the solution for a wave originating from the left is 



A(e lkx + 



k-k' 
k + k 1 ' 



ikx\ 



A 2k c ik'x 



x < 
x > 



(4.12) 



k + k' 

where A + = A for convenience. Now, including the time dependence for x < 0: 



k + k' 

Remark 10: As E — > oo, the amplitude A ^ + y - of the reflected wave 



(4.13) 



while the 

amplitude of the transmitted wave goes to the amplitude of the incident 
wave (eg, the step became negligible). As the energy E — ► V2, then k' — ► 
and the amplitude of the reflected wave goes to the amplitude of the incident 
wave and the transmitted wave amplitude — > since it is preparing to become 
exponentially damped barrier penetration. 



4.4 Finite Square Well: Discrete Energies and Resonances 



Consider the potential 
So we have 

n*) = {£; 

From our previous two examples, we expect the following: 



\x\<L/2 
\x\>L/2 
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v 



• Discrete energies for V\ < E < Vi 

• Barrier penetration for V\ < E < V2 where |x| > L/2 

• For E > V2, transmission and reflection at discontinuities in potential. We also expect 
2-fold degeneracy since we can have waves coming from the left or the right. 

• Resonances as a function of energy in the transmission and reflection of waves. 
To actually solve the problem, we will again consider each energy range separately: 



Vi < E < V 2 : 



The time independent Schrodinger equation is: 

-h 2 d 2 



2m dx 2 
Then the solutions are 



\-{V 2 -E), 



ip(x) 



A + e ikx + A_e ikx , 
B+e KX + B^e~ KX , 



\x\<L/2 
\x\>L/2 



\ X \ < L/2 

\x\ > L/2 



where 



hk = v / 2m(E - Vi 



and 



hn = v / 2m(V 2 - E) 



Applying boundary conditions, we have 
• ^(x — > ±00) < 00 



B + = for x > L/2 
B_ = for x < -L/2 



A+e ikx + A_e ikx , 



ip(x) 



\x\ < L/2 
x < -L/2 
x > L/2 



B + e KX , 
B_e~ KX , 

Continuity of ip at x = ±L/2 

A+e~ ikL/2 + A_e ikL/2 = B + e~ KL/2 

A e ikL/2 + A _ e ~ikL/2 = B _ e -KL/2 



(1) 



(2) 



(3) 



(4) 



(5) 
(6) 
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Continuity of ip' at x = ±L/2 

=> ik [A + e~ ikL/2 - A_e ikL/2 ^j = KB + e~ KL/2 (7) 
ik (A + e ikL ' 2 - A_e- ikL / 2 ) = -KB_e~ KL l 2 (8) 



Now we solve the system: 

1. take k ■ (5) i (7) 

(k + i/fc)A + e^ ifeL/2 + (k - ^)^_e ifcL/2 = 2KB + e~ KL ' 2 (9) 
(k - ik)A + e~ ikL / 2 + (k + ik)A_e lkL ' 2 = (10) 

2. Similarly, take k • (6) ± (8) 

(k + ik)A + e lkL/2 + {k- ki)A_e~ lkL/2 = (11) 
(k - iA;)^ + e ifci/2 + (k + H)A_e" ifci/2 = 2KB—e~ K ^^ 2 (12) 

3. From (10) and (11), we have 

= Al => = ±A_ (13) 

Remark 1: The "plus" in (13) corresponds to an even parity solution (cosines in 
box)' the negative sign corresponds to an odd parity solution (sines in box). 



Even Parity Solutions: Let 



2 

Then from (10) or (11): 

kL 



A + = A„ = ^ (14) 



^ k 2 + k 2 cos [ — + (f>) =0 (15) 



where eft = arctan (&/«;). 



kL f k\ nir 

h arctan I — I = — where n = odd (16) 

2 \ hi J 2 



Also, from (9) and (12), we get: 



Vk 2 + k 2 Acos{kL/2 -</>) = 2KB + e~ KL/2 
Vk 2 + k 2 Acos{kL/2 -cj)) = 2kB^ rL/2 



(17) 



5+ = ^ = ^?co S (f-0)e-/ 2 (18) 
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So the wave function for the even parity solutions is 



Acos(kx), 



A^k 2 + k 2 
2k 

A\/k 2 + k 2 
2k 



cos 



COS 



(kL 

fkL 
\~2 



\x\ < L/2 
^ e KL/2 e KX , x < -L/2 

ij e KL / 2 e~ KX , x>L/2 



(19) 



Odd Parity Solutions: Let 



Aa 



2i 



(20) 



Then, following the same procedure as in the even case, we get the energy 
constraints: 

nir 
~2~ 



kL (k 
h arctan — 

2 \K 



And the coefficients B+ and i?_ are: 



-B- 



. Vk 2 + k 2 . (kL 

-A sin — 

2k 2 



— -rh\ P * L I 2 



This gives us odd parity solutions of: 



(21) 



(22) 



ip(x) 



Asin(fcr), 



\x\ < L/2 



A-J k 2 + k 2 ( kL 



(23) 



Remark 2: The general form of the solutions looks like: 



Remark 3: The allowed energies are obtained from Eq (16) and (21) for the even 
and odd parities (respectively). Consider the even parity energy values. It is 
helpful to eliminate the k so that all energy dependence is in k. As a reminder, 
Eq (16) was: 

nir kL ( k 

— = h arctan — 

2 2 \ k 



Looking at just the trig, component, we have the triangle: 
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(a) Even case: Note that the expo- (b) Odd case: Note that the ex- 
nential decay parts have the same ponential decay parts have opposite 
sign signs 



s/r<? + k 2 




arctan 



Which means that we have: 



arctan 



arcsm 



arcsm 



arcsm 



k 



,Vk 2 + k 2 



hk 



2m(E - Vi) + 2m{V 2 
hk \ 



E) 



y/2m{V 2 -Vx) 
Thus, Eq (16) becomes: 

kL . / 



nir 
~~2 



hk 



arcsm 



\y/2m(V 2 -Vi) 



by Eq (3) 



(24) 



(25) 



To find the actual energy levels is very difficult, as the above is a transcendental 
equation and thus must be solved either graphically or numerically. We can 
look at it graphically by plotting the left and right hand sides of the equation 
and looking for their intersections. 
So solutions for n > 1 require: 



1 



V 



-L 



7T 

1 - < - 

2 ~ 2 



where 



■ A /2m(V 2 -V -1 



L 
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Ihs 
rh s 




► k 



k = ± s /2m(V 2 -V 1 ) 



is the largest k allowed by the rhs of Eq (25). 



^ V2m(y 2 -y!) 

=> — - r > (n — 1)it 

And thus there are a finite number of intersection points since the maximum 

value allowed for k is V 2m ^ YlL (which corresponds to the kinetic energy 
being equal to the depth of the well). Eg: in the figure below there are 7 
allowed values (energies/intersection points). 

, , v 2 



Vi 



E > V 2 : 

Remark 4 : I n this range, there is 2-fold degeneracy. We will construct solutions for 
the wave originating form — oo. 



_> > ) 

< < — 
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Thus, we will have a solution of the form: 



' A + e ik ' x + A_e~ ik ' x , x < -L/2 
B+e ikx + B_e' ikx , \x\ < L/2 

x > L/2 



Ce 



ik'x 



where 



hk = yj2m(E - Vi) and Hk' = yj2m(E - V 2 ) 
Solving the boundary conditions, we have continuity of ip at x = ±L/2: 

A+e -ik'L/2 + A _ e ik'L/2 = B+e -ikL/2 + B _ e ikL/2 
B+e ikL/2 + B _ e -ikL/2 = Ce ik'L/2 

And by the continuity of ij/ at x = ±L/2: 

ik' (A + e- ik ' L l 2 - A„e lk ' L ' 2 } = ik ( B+e~ ikL / 2 - B„e ikL ' 2 



ik (B + e ikL ' 2 - B_e~ ikL ' 2 ) = ik'Ce lk ' L l 2 
Solving these, if we go k' ■ (4) ± (6), we get: 

(k' + k)B + e ikL/2 + (k' - k)B.e ikL/2 = 2k'Ce ik ' L/2 
(k' - k)B + e ikL/2 + (k' + k)B_e~ ikL/2 = 

From (8), we have 



B_ = B. 



k + k'' 



ikL 



Putting (9) — ► (6), gives C in terms of B^ 



k f B + e lkL / 2 - B+^—^e^A = k'Ce lk ' L / 2 
k + k' I 



C 



-2kB + e^ k - k 'W 2 
k + k' 



Now, relate these to A + and Take k' ■ (3) + (5): 

2k'A +e - ik ' L/2 = (k' + k)B + e~ ikL/2 + (k 1 - k)B.e ikL/2 



B JkL/2 



(k + k')e~ ikL 



(k - k') ikL 
(k + hf) 



by (9) 



B ikL/2 

— t — — \Akk' cos(feL) - 2i(k 2 + (A/) 2 ) sin(feL)! 
[k + k ) 
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A i(k+k')L/2 



WpCOs(kL) 



(12) 



Remark 5: Note that 



A + B< C 
> > > 



<J- < 

A- B- 



where A±, B±, C are the amplitudes of the waves. 



Remark 6: We may choose A + to be the overall scale factor of the wave, and thus we get the 
other constants in terms of A + . The most important coefficients are C and A- . 



Now from (11) and (12): 



And we also have: 



c 



-ik'L 



A + cos(fcL) - i(fc2 +g ,)2) sm(kL) 
lk , L (k^-(k'f) sm(kL) 



(13) 



1kk< cos(fcL) - i(fc2 +g ,)2) sin(fcL) 



Remark 7: One can show that 



c 


2 






+ 









(14) 



(15) 



Physically, A + , as illustrated in Remark 5, represents an incoming beam of particles which 
scatter off the square potential. Transmitted particles are represented by the beam with 
amplitude C, and the reflected particles are represented by the beam with amplitude A_. 
So, by (15), the incident particle flux (oc |A + | 2 ) = the transmitted flux flux (oc |C| 2 ) + the 
reflected flux (oc |yl_| 2 ). 



Remark 8: From (13) and (14), we know that if: 

k n L = nir 



c 



1 and the reflection coefficient 



then the transmission coefficient 
gives a series of transmission resonances with resonance energies: 

E„ 



(16) 



0. This 



2m V L ) 



(17) 
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At these energies, the entire wave is transmitted. 

Remark 9: The energies in (17) are not discrete in the sense of the energies confined in the 
well. All the energies in the neighborhood of the resonance energies are allowed, but the 
transmitted wave amplitude is diminished if not at resonance. 



Remark 10: The sharpness of resonances is determined by the ratio between k and k' . Consider 
the case where 



, k , 
•P 



Then 



E-V 2 
E-Vx 



C 



< 1 



Vi < V 2 (a deep well) 



-4, 



cos 2 (kL) + (Jfa) sin 2 (>L) 



(18) 
(19) 



In between resonances, where sin 2 (A;L) ~ 1, we have 



C 



A, 



(20) 



This implies that transmission nearly vanishes between resonances (but not completely!). 



k± Ak 




To approximate resonance widths, let's expand the cos(kL) and sin(feL) in the vicinity of the 
resonance. 

cos(kL) = cos(k n L + Ak n L) = cos(k n L) cos(Ak n L) — sm(k n L) sin(Ak n L) 

w (±1) cos(Ak n L) + (0) sin(A/c n L) 

~ ±1 for small Ak n 
sin(A;L) = sin(A; n L + Ak n L) = cos(k n L) sm(Ak n L) + cos(AA; n L) s'm(k n L) 

« (±1) sm(Ak n L) + cos(AA; n L)(0) 

~ ±A/c r ,L for small Ak„ 
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With this in mind, in the vicinity of the resonance, (19) becomes: 



C 



!+(#-) (Ak n ) 2 L 2 



(21) 



In order to get the width at the half-max point of the peak, we set the above to 1/2. This 
implies that: 



2k' n 

kn 



Ak n L 

<^ 1 for deep wells 

And so we can define a half-width in k by 

2k„ 



Ak r , 



k n L 



(22) 



Remark 11: We can convert to the half- width in energy by substituting (3) into (22): 

2hk' 

y/2m(E - Vx) - hk n = 

k n L 

y / 2m(E-V l ) + 2m{E n -V 1 ) - hk n = — p. 

k n L 

IE-E n \ _ 2hk' 



V2m(E n - V r ) ■ (l + j- f ) - Hk n 

> v ' V ZE n -VlJ 



k n L 



where the term in parenthesis is the first terms of the Taylor expansion, which we can use 
since we are in a deep well and the ratio is very small. Now, if we let AE n = E — E n , we 
have: 

AE u ^ mk ' n{En ~ Vl) 
kn L Kk n 



4 v / 2m(E n -V 2 ){E n -V 1 )h 
L ■ 2m(E n - V x ) 



8 J^LZM (23) 
mL 



Remark 12: Eq (23) is valid only for deep potentials. 



4.5. SQUARE POTENTIAL BARRIER AND TUNNELING 

4.5 Square Potential Barrier and Tunneling 

Consider the following potential 
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V = V 



► 



V 







L 



The methods from §2.2-2.4 can be applied here, though we don't do them explicitly here. 
The new feature here is the tunneling effect, which can occur when < E < Vo- Classically, 
an incoming particle will be reflected at the barrier. Quantum Mechanicaly, a wave function can 
experience barrier penetration. In some cases, this penetration is far enough that the wave can 
emerge from the other side, as seen in Fig. 4.2. 



Thus, the quantum particle can tunnel through the potential. 

Remark 1: The size of the transmitted wave is determined by the height of the potential (which 
sets the decay rate for the exponential solution), and the width of the potential. 



4.6 General Properties of 1-D Schrodinger Equation 

Here, we summarize some of the results of §2.1-2.5. 

Remark 1: Here we are considering the time independent Schrodinger equation in 1-D with one 
restriction on the potential, V{x): 

• V{x) is a real function, bounded from below and piece-wise continuous over 
the real line. 

Remark 2: The following properties can be shown mathematically in Messiah. Here, we just state 
the results. These should be MEMORIZED! 




Figure 4.2: An incident wave tunneling through a step barrier. 
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4.6.1 Asymptotic Behavior of Solutions 

Here, let's state the results for x — > oo; similar results hold for x — > — oo. Choose xq 
so that either 

1. E-V(x) > Vx > x 




2. E - V(x) < \/x > x 



xo 




V{x) 



Case 1: Solutions to the Schrodinger equation oscillate and are bounded as x — > oo. 
If linxj^oo V(x) = V + < oo, (or V(x) < oo), and if it does so faster than -, 
then the asymptotic solution of the equation is: 

tp ~ x -+oc ^4+ sin(A;j; + <f>+), where hk = \/2m(E — V + ) 

Case 2: For E < V m i n < V(x), and for x > xq, one solution decays faster than 

tp ~ e~ KX , where Hk = \f2m{V m i n — E) 

All other solutions diverge at least as fast as e KX , and are eliminated on the 
basis of physical realizability. (§2.3 - Remark 4) 

4.6.2 Nature of Energy Eigenvalues 

Assume that 

V(x) V+, V(x) V+> V„ 

which will look something like: 



V- 




Then there are 3 energy ranges for which the following properties hold: 



4.6. GENERAL PROPERTIES OF 1-D SCHRODINGER EQUATION 



61 



1. E > V+: Energy spectrum is continuous and degenerate (there can be left 
or right travelling waves). The states themselves are unbounded. 

2. V- < E < V+: Spectrum is continuous, but non- degenerate (waves can only 
come in from the left and reflect back to the left). States are unbounded. 

3. E < V-\ If solutions exist, they have discrete energies and are bounded 
(ip — ► as x — > dboo). Spectrum is non-degenerate. 

Remark 1: Non-degeneracy is not a general property of bound states. In multi- 
dimensions, we can have bound, degenerate states. 
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Chapter 5 

Wave Packets 



In the previous chapter, we looked at a variety of 1-D problems, and particularly at solutions with 
waves extended from — oo to oo. Physically, such waves don't exist(though reflection, transmission, 
barrier penetration and tunneling do occur). Instead, we have wave packets. 

5.1 Plane Waves 

The simplest wave with frequency v and wavelength A is a plane wave: 

^(r,t) > 



-i(wt— kr) 



where u = 2nv and k = -y-. 

Remark 1: The wave propagates in the direction of k. 



Definition 1: The phase velocity of this wave is: 



v phase ~ ~ V ^ 



This gives the speed at which planes of equal phase move in the direction k. As a visual 
confirmation, consider a normal set of plane waves and a phase shifted set, as shown below. 



2tt 4jt 



Normal: 



Phase Shifted: 



2jt 
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We can see that the speed of the plane will be 

A uj 
— = Xu = — 
T k 

in both cases. 
Definition 2: The period, T, is T = -. 

Remark 2: We can not really associate plane waves with the motion of a particle (at least, not 
in the classical sense). For this, we need to use wave packets. 

5.2 Wave Packets 

First, consider the motion of a classical particle. Here, the velocity is given by 

p 2E p 2 

Vdass = — P = — ; Since E = for free particles 

m p 2m 

Remark 1: Quantum mechanical particles, like light, possess wave-particle duality, 
in which they can exhibit properties of both waves and particles, though the 
properties can not be simultaneously observed. 

For a quantum particle with energy E, we associate a wave-frequency by 

E = Hlo 



Remark 2: de Broglie developed the theory of matter waves based on the assumption 
that particles have wave properties. 

For a particle with momentum p, we assign a wave number k by 

p = ?ik = hkk 

This is called the de Broglie relation. Since we can identify a wave with a particle, we 
can also identify a phase speed: 

_u~_E~_E„ 

V phase — ~rk — K — p 

k p p 
which differs from the classical case by a factor of 2. 

To find a relation between de Broglie's matter waves and the speed of a classical 
particle, consider 

fc = £ = i^^ 

2m 2m 
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which gives 



And 



Or, in 1-D 



duj A <9u; „ <9u; „ 
kU; dfcj,^ Sfez^ 

= — {k x x + k y y + 
m 

_ _P 

m 



"777" — — — ^class 

ok m m 

Remark 3: A dispersion relation is given by the frequency as a function of k: uj{k). 



Remark 4' Here, we can obtain the classical velocity from the dispersion relation. 
To associate a wave with a particle, we need waves of limited extent (not plane 
waves, but wave packets). Wave packets can be constructed as a superposition 
of plane waves and can be expressed as 

/oo 
d 3 k' /(A/) e - i(a/ *- k '- r) (t) 
-oo 

Remark 5: f{k) gives the "weighting" of the superimposed plane waves. It is gen- 
erally complex, so it defines the amplitudes and the phases of the waves that 
are added together. 

Remark 6: Wave packets that are localized in space(like (t)), move as a whole at 
the group velocity. 

V group — Vk'k- 1 t/ = fc 

where k is the value of k' where f(k') is peaked. The group velocity corresponds 
to the classical velocity of a particle. 

5.2.1 Spacial Extent of Wave Packets 

For now, let us consider a 1-D wave packet. 

/oo 
dk' f(k')e- i ^ t - k ' x ^ (1) 
-oo 
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At time t=0: 

^(x,0) 

A generic form of \f(k')\ is: 



I fdJ\J,k'x 



dk' f(k')e 



(2) 



O- A* 




If we let k' = k + Ak', 

i/)(x, 0) = e 



kx 



dAk' f(k + Ak')e 



i\„iAk'x 



(3) 



The above is akin to shifting the peak to the origin in the Ak' integration. 
Now, consider (3) for different values of x: 
• x=0: 

/•oo 

^(0,0)=/ dAk' f(k + Ak') / 



If f(k') is real, this is just the area under the curve, and thus ^ 0. 

x = f(k') is only significantly different from for the range of Ak' 
given by — A/c < Ak' < Ak. This implies that the phase of the oscillary 
term inside (3) varies over the range —100 < Ak'x < 100. This range 
covers about 30 cycles of 2ir, so that the real part of f(k + Ak')e lAk x looks 
like: 




The imaginary part looks similar with oscillating phase shifted by ir. When 
the integral (3) is done over this function, it is very close to zero, which 
results in a limited spacial extent of the packet. 
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Remark 1: x = gives the largest contribution to the wave function (^;(a;,0)), and, 
for x S> we get almost no contribution. For ranges < x < ^E, there is 
a small, but negligible contribution to ^>(x,0). Take x = 4r: 




And integral (3) doesn't quite vanish because small values of x imply slower oscillations 
which give bigger contributions, while large values of x imply faster oscillations, 
which contribute nearly zero. This behavior gives a modulus wavefunction that 
looks like: 




Remark 2: This acts as an envelope function for e , multiplying the integral in (3). 
So the wave function is a rapidly oscillating function with peaks at x = and is 
only significantly different from zero in the region ~ —Ax < x < Ax ~ 

Remark 3: For instance, a classical particle: 

m = 1 kg, v = 10 m/s, Ax ~ 1CF 15 m 

where the value of Ax is approximately the size of an atomic nucleus. Now, 

mv = hk = 10 kg m /s 

k = 10 35 m _1 

A = 2 4- ~ 10- 24 m 

The range, A/c of k needed to achieve this localization (Ax) is 

Ak ~ - 10 15 m -1 
Ax 
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This implies that we need an accuracy to 21 decimal places if we want more 
than a single wave number to contribute to the wave packet. And thus a 
classical particle has a well defined position and momentum. 



Remark 4 : Now, for an electron: 



m 



x 10~ 31 kg, E = 1 eV = 1.6 x 10" 



-19 



Consider Ak ~ 10 m, which is about the size of an atom. 
hk = V2mE 4fc~2x 10 11 m" 1 



A 



2vr 
T 



3 x 10 -11 m 



Here, the range needed to achieve localization is Ak ~x:~3x 10 m 



10^-1 



Ak 

=^ 0.1 

k 

And so a wave packet living in Ax ~ 10~ 10 is made up of many wave numbers 
and does not have a well defined momentum. 




2 x 10" m" 1 



Figure 5.1: K Space representation of the wave packet. 
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5.2.2 Motion of the Wave Packet 



For a wave packet that represents a particle in the classical limit, the group velocity 
(as defined in §5.2, Remark 6) corresponds to the classical velocity of a particle. 
Consider (1) of the previous section. 



with k! —* k + Ah' and lo 1 



dk' f{k')e- i{u),t - k '^ 



+ AJ 



if)(x,t) 



-i(ujt—kx) 



/oo 
dAk'f(k + & k >y-i{^'t-Mx) 
-oo 



Following the same arguments, we can demonstrate that there is a localized wave 
packet, but, if t ^ 0, the wave packet has moved. 

Consider the phase <f> = Auj't — Ak'x. The position x is given as the position that 
causes <fi to change the slowest as Ak' varies over the relevant range. 







dAk 







Afc'=0 

dAu' 



dAk' 



t — x 



Ak'=0 

And this gives exactly the group velocity of the wave: 



J group 



X dAuj' 

1 = dAk' 



Afc'=0 



dJ 

dk 1 



k'=0 
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.2.3 Delay of Wave Packets at Resonance 

Consider reflection and transmission of a wave packet at a square potential for energies 
E > V : 



What is the delay in the appearance of the transmitted pulse? 
The incident wave packet is: 



^(s,t) = J dk'A+We-W*-^ (f) 



where A + (k') is centered at k' n resonance. 
The transmitted wave packet is: 



From Ch 4, Section 2, Eqn (13): 

C e ik ' L 



A + cos(fcL) - * (fc2 +g ,)2) sin(fcL) 



A + {k')e- i ^' t - k ' :i 



tp T (x,t)= I dk' — 77 — lk 2 +(kl) 2v . (t) 

(2kk') 2 



cos^fcL) + ^'P 2 sin 2 (fcL) 



where 



(f>(k') = -k'L + arctan ( - + ^ tan(fcL)^ 

\ 2kk J 

(Here we just have written the complex number C/A + in its polar form). For sharp 
resonances, k! jk <C 1, and in the vicinity of the resonance we have 



4>{k') « —k'L + arctan ( — \ tan ((k — k n )L)) 
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Resonance Width 



To figure out the time delay, we note that the phase in the integrals (f) and (X) must be 
minimized for the maximum contribution to the incident and transmitted waves. The 
incident wave simply comes in at the group velocity (§5.2.2), while the transmitted 
wave has a delay. 



_d 
dk 



-j(Jt - k'x - 4>) 



k'=k' 



dJ 
~dk> 



t — x 



V 



where 



d(f> 



group 



7T 



k' — k f n 
t — X 



dk' 







dk' 



= 



k'=k' 



dk 1 Res. Width 
by a rise over run argument. Then 



hk d(f) 
m dk' 



Hk 



k' =k' 



111 



(t + r) 



where r is the time delay at resonance. 



m dcf) 
hk dk' 



m 



7T 



hk' Res. Width 



Remark 1: In order for this picture to hold, we assume that the envelope function 
of the wave packet hold together so that there is a well defined peak across 
the barrier. To ensure this, the wave must be broad enough (in x), so that the 
delay is small compared to the motion of the packet: 

f hk Y\, 

t < — A A; 
\m J 
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Chapter 6 

Harmonic Oscillator and Second 
Quantization 



The harmonic oscillator (H.O.) is an important example in QM. A major application of H.O. is 
QED, where elementary excitations (modes) of the electromagnetic field are formally harmonic 
oscillations. This also leads to the generalization of QED, which is quantum field theory. Another 
major application covers all areas where there are harmonic oscillations about some equilibrium 
(like vibrations of molecules in crystals). 



6.1 Harmonic Oscillator in Position Representation 

The Hamiltonian for a single oscillator is 



f> 2 

H = T + V = — — h -muj 2 x 2 
2m 2 



(1) 



where u = yjk/m is the classical frequency. 



Remark 1 : A particle that is experiencing small oscillations about a stable equilibrium of potential 
V(x) may be approximated by (1). If xq is a stable equilibrium about V(x), then 



V(x) « V(x ) + 



dV 
dx 



(x - x ) + 



1 d 2 V 



■!■() 



2 dx 2 



(x - x ) 2 + • 



•I'D 



For small oscillators, 
V{x) 



1 d 2 V 



2 dx 2 



(X - Xq)' 



1 d 2 V 



X = XQ 



2 dx 2 



i 1 ii 
x = -mw x 

2 



with raw 



2 _ d 2 V I 



73 



74 



CHAPTER 6. HARMONIC OSCILLATOR AND SECOND QUANTIZATION 



Remark 2: For an n-oscillator system, the Hamiltonian may be diagonalized to a system of n 
decoupled oscillators, as was done in Problem 1.8.10 for 2 oscillators. In such a system, the 
Hamiltonian may be written as: 

N N PPS 1 N N 

i=l j=l i=l j=l 



with 

d 2 v 



v v = J?- 

dxidxj 



Vji 



The time dependent solution of the quantized HO can be captured by acting the propagator (§3.4.1) 
on the time-independent solution. So our goal here is to solve the time-independent Schrodinger 
equation in the coordinate basis: 



H \ E n ) — E n \ E n ) 



where H is given in (1). In position representation: 

(x\ H \E n ) - (-1- ±- 2 + l - m ^ ) ^(x) 
(x\E n \E n ) ^ Eii> E (x) 
So the time independent eigenvalue equation is: 

d 2 , 2m f „ 1 2 ; 



da . 2 -^ + -p-(^-^V)^ = (2) 



Following the procedure is Shankar: 

1. Non-dimensionalize (2). Let 



— —x and ^ = — (3) 
n nw 



Then 



d Imuj d 



dx \ h dy 
And thus, under the change of coordinates, (2) becomes 



^ E + (2£ - y 2 )^ E = (4) 
dy z 

Remark 3: The scaling defines a natural length scale in the problem ~ -Jjr- 
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2. Investigate the asymptotic behavior. 
Consider lim^o an d lim^oo. 



lim (4) => ^| - y 2 = 
The solution as y — ► oo then is: 

V> = Ay m e ±3/2/2 (6) 



cfjp 
y"-^b x ~' ' dy 

This has a solution of the form: 



lim(4) =► -j-J + 2£V = (7) 



^ = Acos( v / 2ey)+^sin(y2iy) (8) 

Remark 4 : In Eq (7), we have already neglected terms of order y 2 , so we should 
only keep terms with order < y 2 in (8). 

A + B^2~ty + 0{y 2 ) (9) 

3. Make an ansatz of a solution that satisfies the asymptotic behavior. Note that an 
ansatz is like an educated guess as to the solution. 

My) = u{y)e-y 2 ' 2 (10) 

Now, substitute (10) into (4) to derive a DE for u(y): 

dy 2 dy \dy J 

= %-v 2 / 2 - 2 d ^ye-y 2 ' 2 + uy 2 e-y 2 ' 2 - ne^' 2 
dy dy 

= -(2{ - y 2 )ue'« 2 ' 2 by (4) 

-£-*| + «-l>«-0 < U > 

4. Solve (11); Try a power series solution. Let's try: 

oo 

u{y) = Y,C n y 2 (12) 

n=l 
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Substituting (12) into (11) yields: 

= Cn [n(n - l)y n ~ 2 - 2ny n + (2£ - l)y n ] 

n=0 

oo 

= Y, C n Hn - l)y n ~ 2 + (-2n + 2£ - l)y n ] 

n=0 

oo oo 

= Cnn(n - l)y n - 2 + C„.(2£ - 1 - 2n)y n 



n=2 n=0 

oo 

= C n+ 2{n + 2) (n + l)y n + ^ C n (2£ - 1 - 2n)y n 

n=0 n=0 

oo 

= ^ [C n+2 (n + 2)(n + 1) + C n (2£ - 1 - 2n)] 2/ n 

n=0 

Since the y's are linearly independent, the coefficients must vanish and this results 
in a recursion relation for the coefficients: 

-C n (2£ - 1 - 2n) 
(n + 2)(n + 1) 

So given Co and Ci, (13) gives all the rest of the coefficients in (12), and thus the 
solution to (10). ..and thus the solution to our non-dimensionalized Schrodinger 
equation (4). 

Remark 5: The solution given by (3), (10), (12), and (13) doesn't seem to constrain the energy 
eigenvalues as expected. However, recall the asymptotic solution as y — > oo (See (6)). For a 
finite wavefunction, we take the negative in the exponential, and thus u(y) — — ^> y m . But 
the solution (12) does not terminate at finite n unless C n = for some n. The recursion 
relation (13) has C n +2 vanish for arbitrary C n when 

t = n +\ ( 14 ) 

So physically realizable solutions only occur at discretized energy eigenvalues. From (3) 
and (14): 

E = hw{n+ -), n = 0,1,2,... (15) 



Remark 6: Here we have only positive (or 0) values of n since we can easily show that the energy 
eigenvalues can not be negative. 

(h) = ^ (MP^+ l -mu 2 > 

positive positive 
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Remark 7: For n odd, let Co = 0, and for n even, let C\ = 0. Then the solution to (4) is 

My) = u(y)e-y 2 / 2 

' C + C 2 y 2 + C 4 y 4 + • • • + C n y n , n = even 



-y 2 /2 



Ciy + C 3 y 3 + --- + C n y n , 



n = odd 



These are Hermite Polynomials, H n (y) 



Remark 8: Normalizing (15) to get C$ or C\ gives the final solution 



1/2 

exp 



muj 



Trh2 2n (n\) 2 



-mux 
2h 



Hn 



moj \ !/ 2 



(17) 



6.2 Second Quantization: Creation and Annihilation Operators 
6.2.1 Harmonic Oscillator in Energy Eigenbasis 

Remark 1 : Here we solve the HO using an entirely different approach using creation 
and annihilation operators. This method is due to Dirac and is central to QED. 
One advantage is that we do not have to solve any differential equations-we 
can extract the energy eigenvalues from the underlying operator algebra. In 
this approach, we can also extract the physical content from the Hilbert space 
and the operators that act on it, rather than from the wave functions obtained 
from the solution of the Schrodinger equation. 

Definition 1; The Creation and Annihilation operators are given by the following: 

x jmuj / i A 
a = i^h{ q -^ P ) 

respectively, where q is the position operator. The reason for the names will 
become obvious shortly. 

Remark 2: The commutation relation for a and are: 

[a, a] = = [a\ a*] 
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mto 
~2h 



1 % 

H P, Q P 

moo muj 



muj ( 1 i r ^ i fA 

9i9 + — 2-2LP»P + b>9 9.P 

ma; / 2i 



Remark 3: we can obtain the operators p and g in terms of d and d^: 



(a + d f ) 



2mw 



P 



hmu> x 
! Y^— («- a 1 ) 



Therefore, we can write the Hamiltonian as: 
1 / hmu \ 



H 

2m 



— — j (dd — ad' — dd' + a^a^)+ 



H — muJ 2 ( | (aa + dd^ + a) a + d^) 

2 \2rna; / 

= -hu)(aa) + d^d) 

= hj0 ( N+ l \ (t ) 
Where we define the number operator N = d^d. 

Remark 4 : With the Hamiltonian given in the form (f), the problem is reduced to 
finding eigenvalues of N . 

To determine the eigenvalues of N, we can denote the eigenstates of N: 

N \v) = v \v) 

and establish a relationship between the eigenstates. 
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Remark 5: If \v) is an eigenstate of N, so is a\v) (an eigenstate). 

Na \v) = a)aa \v) 

= (aa^ — l)a \v) 
= a(a^a^) \ v) — a\v) 

N 

= au \v) — a \v) 
= (u — l)a \v) 

with eigenvalue (v — 1). Similarly, a) is an eigenstate of TV with eigenvalue 
Na j \v) = {v + l)a) \v) 

Remark 6: The eigenstates \u) , a, \u) , a' \u) of N are also eigenstates of H with 
energy eigenvalues: 



v H — ) huj for | u) 
2 



v — i hui 

2 



for a I v) 



v H — | huo 
2 



for a* 



Repeated application of a or a< generates a whole series of energy eigenstates. 
Eg: 

E-States E-Values 



a 2 \ u) 
a \v) 

W) 

at 1 1/) 



v — |) hu 



i) 



V 



Remark 7: Eigenstates a n will have negative energies if n > v. This would 
correspond to energies below the minimum HO potential, which is physically 
not allowed. To fix this, we require v to be an integer, say n. Then 



N(a n \v)) = {v-n){a n \v)) =0 



for v = n 
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and thus 



H(a n \u)) = -hw{a n \u)) 



where \ fuxi is the energy for the d n \ v) state. But that means that 



H{a n+l \u)) = --hu(d n+1 \v)) 



where we have negative energy-which we just argued was physically impossible. 
How is this possible? Let us calculate the norm of a n+l \ u). To see how this 
works, we'll go back to the wavefunction language. Thus, we have a \ip) — > dip. 



, , „ , , raw / K d . , 

(aipaip) = — ( (q + —)ip 

In \ muj Oq 



. h d , 



mu dq' 



raw 

-2K ' dq 



muj 
N 



Nip 



h d . 

+ — — 1 ^ 



dqip* [q 



muj dq 
h d 
mu dq 



h d 

+ 7T 

mu oq 

h d\ . 

+ 7T ^ 

mu; oq I 



by int. by parts 



Hence, the norm of the state aip is the expectation of N in the state ip. 

Therefore, the norm of a(a n \ v = n) equals the expectation value of N in the state 
a n \ v = n). But this equals zero, since a 11 \ v = n) is an eigenstate of N with eigenvalue 



equal to (as shown above). Thus we conclude that the state a 



n+l 



n) is identically 



0. Therefore, as long as v is an integer, the states with negative energy eigenvalues 
are not generated (So our z/s must by integer to be physical). 

Remark 8: To summarize: The eigenstates of N are represented by |0) , |1) , |2) , . . . , \ n) 

with eigenvalues 0, 1, 2, . . . , n, These states are also energy eigenstates with 

energy eigenvalues ^fouj, ^huj, . . . , (n + ^)hu>, . . .. 

Remark 9: We can see why N is the number operator: It's eigenvalue "counts" the 
number of hu energy quanta in a given eigenstate. Similarly, we can see why 
a and a) are the annihilation and creation operators. Prom Remark 7: 



a n) 



has norm = (n| d^d \n) = (n\ N \n) 
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and 



a 2 \n) has norm = (n\ (a^) 2 a 2 \n) = (n\ a^N a \n) 



E-state w/ 
E- value n— 1 



(n — 1) (n| a) a \n) 
n (n — 1| N \n — 1) 



Thus we deduce that 



a \n) = C_ \n — 1) 
C+\n+l) 



a) In 



with C- and C+ constants. We can see that the a removes (annihilates) one quanta of 
energy while a) adds (creates) one quanta of energy-and hence their respective names. 
To find C+ and C_, we take (n|n) = 1 and use: 



(n\ a' a |n) = (n\ N \n) = n (n\n) = 
= (n-l\ClC- \n- 1) 



C n 



|C_| 2 (n 



lln- 1) 



n 



IC- 



So that 
To find Ca 



a n) 



|n 

aa^ |n 
aa^ In 



{a) a + 1) |n 
(n + 1) \n 



n n 



C+|n+l) 
aC+ |n + 1) 
C+Vn + 1 |n) 
C+Vn + 1 |n) 
CWn + 1 In) 



So that 



Vn+ 1 
|n) = \/n + 1 |n + 1) 



Remark 10: we can construct any eigenstate of |n) from the ground state. 

(at)™ 



n > 



10) 
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2 Eigenfunctions in the Position Representation 

We can apply the above method to obtain the eigenfunctions in the position represen- 
tation without solving the Schrodinger equation. We have 

o|0) = 

which, in position representation, corresponds to 

( K d \ , 

aipo =>• [q-\ — = 

y muj dq J 



dipQ muj 

=>- —r- = ^q^o 

dq n 

dibo f mm 
—qt 



. . , . muj 2 „ 
hi{^)=-—q 2 + C 

tpo = Ae~^ q2 
where A is a normalization constant. 



i = J «*, =a*j « . . * = ^_ =► a = (_) 

Giving the ground state wave-function: 

The excited states can be obtained by repeatedly applying the creation operator: 

\n) = -±=(tf)»\0) 
Vn! 

which, in the position representation, corresponds to: 



/Mi;\ 1 / 4 1 fmu>\ n / 2 ( h d 
This is equivalent to: 



/mw\ 1 / 4 1 1 _mw„2„ / /mw 



where ff n is a Hermite polynomial. 
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6.3 Coherent States: Minimum Uncertainty Wave Packets 

6.3.1 Uncertainty Relations in Energy Eigenstates 

Let's calculate the Heisenberg uncertainty relations in the energy representation. From 
§6.2, Remark 3: 



(n| q In) = \ (n\ (d + d^) \n) 

V 2mu) 



2moj 



(n| (\/n \n — 1) + y/n+ 1 \n + 1)) 



Similarly, 



(n\ p \n) = 



Remark 1: This result is expected from the symmetry of the stationary HO wave- 
function ip n (q)'- 




Ml) 



The uncertainty in position is given by: 

(n| q 2 In) = — — (n| (d 2 + dd^ + d^d + (d^) 2 ) In) 
2mu! 

h 



2muj 

h 

2moj 



(n| (V™(™ - 1) |n - 1) + (n + 1) |n) + n|n) + yfjn + l)(n + 2) |n + 2) 



(2n + 1) 
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The uncertainty in momentum is given by: 

hmuj 



{n\p 2 \n) = '- :i:: —(2n + 1) 



Remark 2: We can rewrite these as 



(n\ q 2 \n) = — (n\ p 2 \n) = mE n 



AxAp = -J{n\ q 2 \n)J(n\p 2 \n) = \ — % = -(2n + 1) 



So the minimum in uncertainty according to the Heisenberg relations is h/2. 
This is satisfied by the ground state-all other states are not minimum uncer- 
tainty states. 



.3.2 Minimum Uncertainty States 

It is possible to construct non-stationary minimum uncertainty states. It turns out 
that these states are eigenstates of the annihilation operator. 

a \a) = a \a) 

where a is a complex number. 

Remark 1: We must make sure these states exist. Expand the state \a) in terms of 
the stationary states | n) . 



oo 

= ^ C n \r 

n=0 



oo oo oo 

o l a ) = CnVn \ n — 1) = ^2 Cn+w/n + 1 |n) = a C n \n) 

11=1 n=0 n=0 



C n+ Wn + 1 = aC n 



a 

C n = —p=C() 



which gives 



v'l 



OO 

a/=C V —j== \i 

n=0 

and Cq is obtained by normalizing: 



00 I I 2n 



1 = (a\a) = \Cn\ 2 = |0>| 2 £ = |C | 2 e H 



n! 

n=0 
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=HOj| 




C = e 



which has an arbitrary phase chosen to be 0. Thus we have 




n) 



which is known as the coherent state with complex amplitude a. 

Remark 2: These states were first constructed by Schrodinger as minimum un- 
certainty wavepackets. They are now used extensively in quantum optics- 
coherent states are important in describing optical coherence. 

Remark 3: We can show that \a) is indeed a minimum uncertainty state (See HW!) 

6.3.3 Physical Meaning of the Coherent State 

Let us relate the coherent state \a) to the ground state: 



Consider now |0) = e a a |0) where |0) is the only constant term in the exponential 
expansion that is not zero. Then 





e -M 2 /2 e «at ^ 




Then 



e A e B e [A,B]/2 



So let A = era' and B 



a*a. Thus we have: 



aa) ,—a*a] = \a\ 2 [a, af] 



= \a 



2 



e 
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and thus 

\a) =e aat - a * a |0) =V(a)\0) 

where T>(a) = e aa ~ a a is called the displacement operator. 

To see what this means, assume a = x (a real number) for simplicity. Then 



f)f x \ _ e x(at-a) = e x ^J hk^i-pi) 

pur d 

— g V tnu) dq 

-Q-f 



where Q = x\l ^ is a number with dimensions of length. (Which implies a and x 
are dimensionless. The displacement operator acting on an arbitrary function of q is: 

*^) = «^W = (i-«| + 5« , |i + -)^) 

= F(q - Q) for small Q 



Remark 1: So the action of D on an arbitrary function is to displace the function 
by an amount Q (and hence the name). 

Remark 2: Expressing \a) = D(a) |0) in the coordinate representation: 

Wavefunction of D(a) in coord. Wavefunction of 
coherent state representation HO in state |0) 

where \a)=\x) with a=x 

-OA. /niLOX 1 / 4 mw„2 

V irn J 



\ 7rh J 



So the coherent state has the same gaussian wavefunction as the HO ground state, it 
is just displaced a distance Q = x^J ^j. 
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► i 



Remark 3: The coherent state is not stationary (it evolves with time according to 
the Schrodinger equation). The formal solution is given by 



-m 



\a(t)} = e h 1 - \a) 

Init. State 



—frjjj (a) a + - ) 



cxp 



exp — iuj [a , a + -\ t\ exp 



exp 



\a\' 

~2~ 



]r ex p (-^ (n + J) 



n > 



n=0 



e 2 



|(^t+|a| 2 ) 



n=0 

°° / - --iojt\ n 



£ 

n=0 



(ae-*"*) 1 



-iut/2 



\e- iut a) 



overall arbitrary 
phase can be 
dropped 



Which gives the evolution of (q(t)} and (p(t)): 



(e~ iult a\q\e- iu;t a) 



2mw 



2muj 



2mu 



'e~ iu)t a\ a \e~ iult a) + (e'^al d) \e~ iujt a) 



{e~ iujt a + e iujt a*) 



(cos(a;t)(a + a*) — i s'm(u>t)(a — a*)) 



And so at t = 



<$(0)) 



2?nu; 



(a + a*) 
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And we can show that 

Thus we have that 

(?(*)> = (?(0) > cos(wt) + — (p(0)> sin(wi) 

m<jj 

And similarly 

= <p(0)} cos(u;i) - mw (g(0)) sin(a;t) 

which are the solutions to the Heisenberg equations of motion for the operators q and 
P- 

Remark 4 : The time dependence of the coherent state can be represented as a 
trajectory in the complex plane. 



Im(c<!)) {£(0)' ^Itt'mut 




And we see that the circular trajectory in the complex plane is related to the phase 
space trajectory of classical Hamiltonian mechanics. 

Remark 5: Consider the coherent state wave packet in coordinate representation. 
From Remark 2, we can see (q) = Q (the mean position), while all coherent 

states have Aq = (Remark 3 - HW). And the wave packet moves in a 

harmonic potential without changing it's width. The coherent state of the HO 
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is a non-dispersing wavepacket. This is significant since we can recall that free 
particle wave-packets always disperse (Ch. 4, §4.1). 



Remark 6: We can calculate the explicit time dependence of the wave packet. Re- 
mark 3 implies: 



\a(t)} = e^ Ht \a) 

= D(e- iu *a)\0) 



e~ iujt a) 



where 



D{e~ lu>t a) = exp 



e~ iwt aa) - e iujt a*a 



In the coordinate representation, the time dependent coherent state is represented by 
the wave function: 



irh 



1/4 



e 2ft 



Displacement 
Operator 

To find what the Displacement operator actually is: 



4 ^(^«-^ )+ ,(^_) 



exp 



ih 

= e imt)) e - ym)) e - <p(*)> <?(*)> [q,P] 
= e -4H (*)> <«(*)> e ^tf (*)> e - 

= e i(,(i)> <p(t)> e i9(p(t)> g- <?«> £ 

The last term is the displacement operator that we had when we let a = a* . If we 
have an arbitrary a, then the displacement operator generalizes and acquires time 
dependent phase factors. 

The coordinate representation of the time dependent coherent state becomes: 
^-**«(?)= e^ m)m \ (^) V4 e-W(<i-m))\ 

overall phase this phase carries Gaussian envelope 

changes with time info about mean carries info about 

and is indep of q momentum mean position and 

Aq,Ap 
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In the last comment, we can find Ap by knowing that \ot) is a minimum uncertainty 
and knowing Aq. The related picture is a gaussian wavepacket that oscillates back 
and forth along the classical HO trajectory with a fixed width. 




\ 1 ►<! 

Q 



If we want to represent this by the coherent state trajectory in the complex plane, the 
width of the gaussian can be related to the radius of a small circle. 

lm(o(t» 




Remark 7: There is another class of minimum uncertainty states which are closely 
related to coherent states. These are known as Squeeze States. 

As minimum uncertainty states, the squeeze states satisfy AqAp = h/2, but they 
differ from coherent states since they have either: 



Aq < 



Ap < 



2muj 



hmuj 



and 



or 



and 



Ap > 



Aq > 



h 

2mu 
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while preserving the minimum uncertainty. In this case, the uncertainty in p or the 

uncertainty in q is less than it is for a coherent state. 

To visualize the squeezed states, define quadrature phase operators: 

xi = 2 {a+a) = y^h q 

In a coherent state, these operators have uncertainties Ax\ = Ax2 = \. In a squeeze 
state, Axi < \ or Ax2 < \ such that AxiAx2 = \- Graphically, the minimum 
uncertainty states lie along the line AxiAx2 = \. 

A 




1 ► Ax " 

All other states along the line are squeeze states. 

A complex plot, in this case, can be represented by a plot on the x\ — X2 plane. In 
this picture, the squeeze states are represented as: 




ii • Quadrature Squeezed x 2 - Quadrature Squeezed 

The time evolution of a squeeze state is analogous to a rotation of frequency lv in the 

complex plane (in this case, the x\ — X2 plane). This is seen in Fig. 6.1. 

The wave packet in coordinate representation has a gaussian envelope whose width 

changes as it oscillates in the harmonic potential, as can be seen in Fig. 6.2 

More generally, the coordinate representation is squeezed twice during the cycle, 

though not necessarily at the end or middle of an oscillation. Eg: See Fig. 6.3 



92 



CHAPTER 6. HARMONIC OSCILLATOR AND SECOND QUANTIZATION 



x 2 

A 




Figure 6.1: The initial squeezing is in x±, which evolves in time to a squeezing in X2- 




Figure 6.2: The outer gaussians are squeezed in position, while the center gaussian is squeezed in 
momentum. 



A 




Figure 6.3: Example of squeeze cycles not fully squeezed on the main axes. 



Chapter 7 

Systems with N Degrees of Freedom 



7.1 Tensor Products (Direct Products) 

Definition 1: Let be a basis vector of the Hilbert space V\ so that the set of all spans 
V\. Similarly, let the set of all \u2) be basis vectors that span V2. Then \u\) 8> \uJ2) forms a 
basis vector in the combined Hilbert space V± ® V2. 

Remark 1: Physically, V\ and V2 may correspond to Hilbert spaces for 2 individual particles 
(1 and 2), while V\ ® V2 represents the Hilbert space for the combined 2 particle system. 
Alternatively, this may represent single particle system that spans 2 (or more) dimensions. 

Remark 2: Consider an operator Cl\ that operates on the space V\ with eigenstates and 
eigenvalues u\. Eg: 



where the superscripts denote the target Hilbert space (either V\ or V2) on which the 



In the combined system, we write: 



Cli 8) \co 2 ) 



U\' \Ul) (g) |W2) 

/ (2) 8> 
(8) /( 2 ) |w 2 ) 




operator acts, and where fi^ o>i 




\ui). Similarly, we can write 



nf) = n(W) = fd) g, ftp) 



So that 
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The bottom line is that Clf' only acts on the basis vectors \uJi) which span Vi- 



Remark 3: It is common to use a compact representation of the vectors which span the combined 
space. 

® \UJ2) = \u> 2 ) = \u>lU>2) 

Since this also represents the simultaneous eigenkets of operators fifp and Q 2 . 

Remark The inner product of the direct products 

% (uj' 2 \) (8) |^2)) = (u[\uji) (u/,|w 2 > 

Similarly, 

((^| ® (n^ ® Aj, 2) ) ® |w 2 )) = ^ |wi) <4| Af |w 2 ) 



Remark 5: Properties of the Tensor Product of Operators: 
Denote 



/( 2 )®A( 2) ^A 2 



acts only on |cui) 
acts only on \uj 2 ) 







Squares of sums: 



(1) /..\(2) 

)) <8)(rAj^ 



(fii + fia) = (^ (2) + S7 



(1)0(2) q(1)®(2) y 



(1) 



(2) 



Remark 6: The time evolution of the state vectors that are composed of direct products is still 
governed by the Hamiltonian: 

for a two particle system in 1 dimension. If the Hamiltonian is separable: 

H = H l + H 2 ^ V(x u x 2 ) = Vx{ Xl ) + V 2 (x 2 ) 
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then the two particles will evolve independently of one another, giving the time evolution 

as 

\iP{t)) = \Ei) e~ iElt / h <g) \E 2 ) e~ iE2t / h 

which is obtained by solving the individual Schrodinger equations (with E = E\ + E<i). 
Separation of variables is often employed in solving the Schrodinger equations in position 
representation. If the Hamiltonian is not separable (ie. V(x±,X2) ^ Vi(xi) + V 2 (x2)), 
then there is not a general method for solving the Schrodinger equations unless V(x\, £2) = 
V(\x± — X2I), in which case the problem can be separated into center of mass plus the motion 
about the center of mass (see Central Potential, coming soon!) 

7.2 Identical Particles 

Classically, two identical particles can be distinguished by fallowing their trajectories (ie. by looking 
at their non-identical histories), which can be done without disturbing the system. Quantum me- 
chanically, no such trajectory exists for particles, and thus there is no physical basis for distinguish- 
ing 2 or more particles. This implies that the same state vector must describe two configurations 
which differ only in particle exchange. 

7.2.1 Symmetric and Anti-symmetric States: Bosons and Fermions 

Consider the following state vector: 

\ip) = \xi = a, x 2 = b) = \ab) (f) 

which represents the result of a measurement that finds a particle at a, and one at 
b. If the particles are distinguishable, then the state (f) is distinct from the state 
\ip) = \x\ = b, X2 = a) = 1 6a). However, if the particles are indistinguishable, then we 
can not differentiate one state from another. Instead, we write the state vector as a 
superposition of these possible states. (These will define a 2-D degenerate eigenspace.) 

\tp(a,b)) = a\ab) + (5\ba) 

These two states are physically equivalent, so we can also write 

WM))=tWM)) 

where 7 is an arbitrary phase (complex number). 

=> a \ab) + (3 \ba) = 7 (a \ba) + (3 \ab)) 

=5> a = 7/3, (3 = 7a 
=4> a = a"/ 2 7 = ±1 
Taking the "+", we define the symmetric wavefunction: 



\tp(a,b)) s = \ab) + \ba) 
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and the "-" gives the anti-symmetric wavefunction: 

IV'O, b)) A = \ab) - \ba) 

Definition 1; Particles whose wavefunctions are symmetric are called Bosons , and 
those with anti-symmetric wavefunctions are called Fermions . 

Remark 1: One of the primary consequences of the anti-symmetric wavefunction is 

the Pauli Exclusion Principle: 
Suppose we have a 2 fermion system with the following state vector: 

so that one fermion is in state ui, and the other is in state w 2 . If u\ = ^2 =^ 
\i/j(ujx, u-2)) a = 0, which implies that no two fermions can be in the same 
quantum state. 

Consider the 2 particle Hilbert space Vi®2 made up of all vectors of the form \ui) ® 
|oj 2 ). Since each pair of vectors \w\ = a, w 2 = b) and \u>i = b, u>2 = a) can form one 
symmetric and one anti-symmetric state, if follows that Vi®2 can De decomposed into 
symmetric and anti-symmetric parts: 

Vi® 2 = V s g) Va, 

Bosonic Fcrmionic 
Hilbert Hilbert 
Space Space 

This leads to an ambiguity in determing whether the state of the system is a member 
of the symmetric Hilbert space (V5) or the anti-symmetric Hilbert space (Va). (We 
will see in §2.2 that the measurement made depends on which case applies.) To remove 
the ambiguity, we introduce the Symmetrization Postulate, which states: 



The state of a system containing N identical particles is either all symmetric 
or all anti-symmetric with respect to permutation (particle exchange) of the 
N particles 



The Corollary to this is that the symmetry (or anti-symmetry) will be independent 
of the number of particles. 

Remark 2: The normalized states for a system of n identical bosons is given by: 

|1, 2, . . . , n) s = —= (|1, 2, . . . ,n) + |2, 1, 3, . . . , n) -\ + n! permutations) 

vn! 

where 1,2,... , n are generic state labels. 



IDENTICAL PARTICLES 
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Remark 3: The normalized states for a system of n identical fermions is given by 
the Slate Determinant: 



i>i,2,...,n{x\,x 2 , ...,x n ) 



ipl(x 2 ) 



The determinant gives all the correct minus signs in the anti-symmetric wave- 
function, and vanishes if any two j's or x^s are the same. 



7.2.2 Distinguishing Fermions and Bosons 

The are two ways to identify whether a particle is a fermion or a boson: 

1. Spin (intrinsic angular momentum) 

2. Experiments 

We will discuss spin in much more detail later, though we note here that bosons have 
integral spin (0, h, 2h, . . .), and fermions have half-integral spin (|, 4^, . . .). 
To determine the nature of the particle experimentally, place two identical particles 
in a 1-D (for simplicity) box. Let 



ips{xi,x 2 ) 
ip A (xi,x 2 ) 



The 2 particle wave function for bosons and fermions 



In general, each particle could have different energy levels, so 

4>S,A = {tjj n (xi)tjj m (x 2 ) ± 1pn(x 2 )lp m {xi)} 

where n,m are energy levels and the "+" is for the symmetric cases and the "-" for 
the anti-symmetric cases. Also, n ^ m (they don't have the same energy). 
The probability distribution in x-space is determined by 

Ps,a( x 1' x 2) = 2 IV , 5,a(xi,x 2 )| 2 (t) 



Remark 1: The factor of 2 in (f) is a consequence of the normalization condition. 
We require: 

1 = jjdx 1 dx 2 \tps,A (xi,x 2 )\ 2 
and accounting for the double counting in the x% — x 2 integral: 

1 = JJ -Ps,a{x\,x 2 ) dxi dx 2 
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where 

Ps,A(xi,X2)= {xiX 2 \lpS,A)s,A 

= V2tp SiA (x 1 ,x 2 ) 

where the y/2 in for convenience. 
Taking (f): 



l(x 1 ,X 2 ) = 2 — (tp n (.Xl)lp m (x 2 ) ± tpn(x 2 )lpm(x 1 ] 

= \lpn{xi)lpm(x 2 )\ 2 + \i) m {xx)lp n {x 2 )\ 2 ± 

± ^(xi^^OV'm^V'n^) + Vn(xi)^ m (xi)^(x 2 )V'n(x2)| 



Remark 2: So the differenence between fermiens and bosons is detectable in the in- 
terferecence in the probability distribution. In the extreme cases, lets consider 
the limit that these particle have as x\ — > x 2 — > x. Then 

Ps{x\ — > x, x 2 — > x) — > 2{|^(x)| 2 |V m (x)| 2 + |^ n (x)| 2 ^ m (x)| 2 } 
P J 4(xi — > x, x 2 — > x) — > 

Where the antisymmetric case is consistant with the Pauli exclusion principle. 
Thus, two fermions will show probability of being in the same state, while 
two bosons double the probability density for two distinct particles to be in 
the same state. 



Chapter 8 



Classical Limit and WKB 
Approximation 

Correspondence Principle: In the limit that h — ► 0, laws of quantum mechanics must reduce 
to the laws of classical mechanics. 

Remark 1: h gives a measure of the energy spacing between discrete energy levels. Thus, in order 
to justisfy the limit as h — > 0, the system must have large quantum numbers. 

Remark 2: The limit as h — > is called the classical approximation, and the conditions of validity 
are the same as those of geometrical optics. 

Remark 3: Note that in the classical limit, we have: 

• [ihVi] = iftSij 

• AqiApi > | — > 0. This implies 

Aft = Vk 2 } - (di? 
which is the same as ignoring fluxuations about the mean. 

In order to ignore the effects of fluxuations, we require that: 

• The mean values approximately follow classic laws of motion, 

• Dimensions of the wave packets are small (and remain small as the system 
evolves) compared to the characteristic dimensions of the problem. 

Remark J^: Since wavepackets disperse with time, the classical approximation is, in general, only 
valid over a finite time interval. 

Remark 5: The correspondence principle, as stated above, is most useful when applied to com- 
mutators, uncertainty, and differences is energy levels, angular momentum, etc... A second 
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formulation of the classical approximation basically amounts to: 

lim Schrodinger Eqn. — > Eqns. of Classical Mechanics 

Geo. Optics 

In this view, the wavefunction represents a statistical mixture of classical systems. The 
density of this mixture, at some point in configuration space, equals the probability density 
of a quantum system at that point. 



8.1 Ehrenfest's Theorem 

Consider a system which evolves according to the Schrodinger equation: 

^H ^ 4 = ih 9 4- = -my 

dt dt x ' 

And the mean value of an observable A is: 

' A) 



Since tp in general evolves in time, \A) will also evolve in time: 

: 'A\ = imA\ 



dq n \ ^ r M + ^A-+^—^ 



dt \ I dt 
d 

~ ~dt 



dqi... dq n ifj*Aip 

j dqi...dq n (-±(HWAi/> + 4*A(J^H4) + 1>*^1> 
J dqi...dq n (-^f(HA)^ + ^r(AH)^ + r^A 



Thus, if A has no explicit time dependence, then 



Remark 1: An immediate consequence of (f) is that if [A, H] = 0, then A does not change with 
time and is a constant of motion (a conserved quantity). This is analogous to classical 
mechanics where we have {A,H} = (The Poisson bracket from §2.1.3). 
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To connect with classical mechanics, apply (f) to qi and pf 

Now consider the $j and p, commutator with 

A(q 1 ,q 2 ,...,q n ,Pi,P2,---,Pn) 
where A is only well defined when the order in each term is specified. We note that: 

qi,F(qi,q 2 ,...,q n ) =0 

Pi,F(pl,P2, ■ ■ ■ ,Pn) 

qi,G(pi,p 2 , ...,p r , 



(tt) 



ih 



dG 



Pi,F(qi,Q2, ■■■An) 

Generalizing the above, one can show that 



dpi 
dF 



-ih 



dqi 



Pi, A 



ih 



dA 
dpi 
dA 
dqi 



-ih 



Using this result in (ff ) gives: 



di 



(fit) 



d 

dt {Pl 



8H \ 

dpi J 

IdH 

\ dqi 



(©) 



which are the quantum averages of Hamilton's equations of motion. 

Remark 2: Equation (©) represents Ehrenfest's Theorem , which states that equations of motion 
for the quantum mean values coordinates and momenta are equivalent to those of classical 
mechanics. 



Remark 3: 
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WARNING! jCUIDADO! ACHTUNG! 

This is not exactly equivalent to saying that mean values of p and q follow classic 
trajectories in phase space. For this, we require that jM- and jM- to be at most 
linear functions of q and p since: 



(QiPi) + (<?;) (Pi) 




unless fluctuations about the averages can be neglected. 



8.2 Classical Limit and Wavepacket Spreading 

Though Ehrenfest's theorem provides a formal connection between quantum mechanics and classical 
mechanics, it doesn't allow us to think of the motion of a quantum wave function as we do the 
trajectory of a classical particle. 
To a good approximation, we can do this if 

1. Quantum means follow the classical trajectories 

2. Aq and Ap are small on a scale of dimensions of the trajectory. 
Using point 2 above to justify the approximate validity of point 1 (they are related), let 

q = (q) + 5q 
P = (p) + Sp 

Then 

(qP) = <«?) +6q)«p) + 6p)) 

= (q) (P) + (SqSP) 

Where the second term is approximately AqAp. If this term is small compared to q and p, we can 
approximately absorb it into (q) {p). For example, consider: 

H = ^ + V{q) 
2m 



Then 



Jt {4) = ^ and &™ — \-w/ 

V(q) = V((q) + Sq) 

^V((q))+V'((q))Sq+^V"((q))Sq 2 

^V(q)^V((q)) + \v"((q))(5q 2 ) 
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If the (8q 2 ) is small enough to neglect, then 

5»— ^««» 



Remark 4 : If (dq. 2 ) starts out small, does it remain small? To answer this, we use equation (f) 
from 88.1: 



1 

ih 
l 

ih 
1 

ih 



5q 2 ,H 



$q , 7T" 

2m 



(fi " <<?» 



2 r_ 

2m 



Similarly, 



((qp + pq) -2{q) (p» 



d 2 2 1 
di 2 ^ ^ ' i/im 



qp + pq -2 (q) p, H 



m 



(pf 
2m 



V((q))-V"((q))(l 



Let's let 



e = ( H 



(pT 

2m 



V((q)) 



Classical Energy 



dt 2 ^ ^ ^ 771 



One can use the above to determine the error introduced by replacing 



(ttt) 



— — \ by -rzVUq)) in Ehrenfest's Equation 
dq / dq 
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Remark 5: For a harmonic oscillator with potential V(q) = |mw 2 g 2 , 



(ttt) => 

while the equation for the mean is: 



dt 2 



4£ ■ .2/^2\ 



m 



4w 2 (5q 2 ) 



which we get from solving Ehrenfest's equations for the harmonic oscillator potential. The 
two above equations imply that (q) oscillates at frequency u>, while {Sq 2 ) oscillates at a 
frequency of 2u. 




But this looks just like our squeeze states! 



Remark 6: For a free particle, V{q) = 0: 

^ dt 2 



'(?) (pY 



2m 



2m 



const. 



:^ 2 > it) 



+ (5q 2 ) (0) 



t=o 



This implies that the wavepacket disperses in time, and thus the classical approximation 
remains valid only for a limited time. 



8.3 Classical Limit of Schrodinger Equation and the WKB Ap- 
proximation 

The alternative formulation of the classical limit is to make a connection between the flow of 
probability density by the Schrodinger equation and the flow of probability density in classical 
statistical mechanics. 
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8.3.1 Quantum Mechanical Probability Current 

The quantum mechanical probability density in quantum mechanics is given by 

P{r,t) =^*(f, #(r,f) 



J P(r,t)dr = 1 



Now imagine that the probability density represents a particle or fluid density. Then, 
the evolution of the probability density (according to the Schrodinger equation) may 
be thought of result of fluid flow: 

dP dP „ -> 

* - a +v ' J = 

where J is the probability current and the whole expression equals so that normal- 
ization is conserved. 

dP 



But 



dP d 

dt dt w w > 

= -i(£^)* + -U*(ify) 
in in 

= l -(im^-rm)) (t) 

And so J must be defined to agree with the right hand side of (t). 

Remark 1: There is some ambiguity in the definition of the probability current 
since satisfying (f) does not uniquely determine J. We choose a symmetrized 
definition: 

2 \ m \m J 



2 \ m m 

-i5.(^*v^-W) (t) 

2m 



Remark 2: Jfr- is chosen so that ($) looks like velocity times density as in classical 
current density. 
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So —V • J, with definition (J), gives: 
ik 

— V • J = - — (W* • Vt/? + v*v 2 v> - • V^* - */>V V) 

= Zl( ^*|%_^|?^/,* 
n \ 2m 2m 

dP 

= — — <— probability flow 
at 



8.3.2 Connect Probability Current in Classical Mechanics to Quantum 
Mechanics 

To formally connect the flow of probability in QM to that in CM, write the wave 
function in terms of amplitude and phase variable. 

ijj(r, t) = A(r, t) e iS ^' h 

amplitude P hase 

Substituting this into the Schrodinger equation to obtain 2-coupled equations: 

d - —h 2 

ifi— ib = H?P = V 2 V> + Vib 

ot 2m 



ih d ( Ae - s / h ) = ( —V 2 + v) Ae~^ h 
at \ J \ 2m J 



Equating real and imaginary parts: 



^A = Jl VA . VS -—av 2 s 

ot m 2m 
^- = -± iV S) 2 -V + ^V 2 A 



(t) 



dt 2m 2m 

From this we can obtain an expression for the equation of motion for the probability 
current. 

P(r) = = A(rf 
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at at m m 



0_ 

dt 



-A 2 = --V- (A 2 VS) 
m 



And we can identify the probability current density: 

1 BP 
J = —A 2 VS => — = -V • J 



Remark 1: The classical limit here is taking h — > in equation (f). This implies: 

dS 1 



2m 



(vsy - v 



(tt) 



This gives the flow of probability density in phase space. The S that satis- 
fies (ff) is the principle function of Hamilton as seen in the Hamilton-Jacobi 
formulation of classical mechanics. (See Phys 505!) 

Remark 2: For the special case of stationary states (Hip = Eip): 

dA 



dt 







OS 

and M 



-E 



and equation (f) becomes 

V • (^ 2 V,S) = 
(VS) 2 + 2m{V - E) 



h 2 V 2 A 



8.3.3 WKB Approximation 

Here, the idea is similar to the classical approximation, but less extreme, since we 
keep the h dependence to the lowest order. Taking the wavefunction: 

V'(r) = A(r)e iS{r)/h 

and write this as 

ip(r) = exp ( — W(r) 

\h 

with 

W{r) = S{r) + -M{A(r)) 

Remark 1: A(r) and S(r) still evolve accordingly to (f) in §8.3.2, but they are no 
longer restricted to be real functions. 
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The idea now is to expand our W(r) to order h 2 and drop higher level terms. Take, 
for example, the case of stationary states in 1-D: 



d 
dx 



dx ) 



(fV + 2^) 



E) 



h 2 d 2 A 
A dx 2 



(1) 
(2) 



(1) ^ 2A—— + A 2 — 

dx dx dx 2 

2 dAdS | A d 2 S 

dx dx dx 2 

1 dA , 
— —dx 

A dx 



ln{A) 



_d_ (dS\ 
dx \ dx ) 



dS 
dx 



dx 



A 



1, fdS 
In — 

2 \dx 



+ constant 



Substituting this into (2): 



dS\ 2 
dx J 



= 2m(E - V) + h 2 I 



fdS_V 12 ± 
\ dx J dx 



1 /dS 

2 V dx 



2m (E -V) + h z [ d ^ 
dx 



2m(E -V) + h 2 



1/2 



d 2 S' 

dx 2 

dS_ 

dx , 



dS\ 
dx J 



-5/2 



dx 2 
d 2 S^ 



dx 2 



dS_ 

dx 



^ 2 H 
dx 3 



d 3 S' 

dx d 
dS 
dx , 



(3) 



Remark 2: Solving (3) exactly is equivalent to solving the original Schrodinger 
equation. The WKB approximation enters here: we begin by expanding S in 
powers of h 2 . Let S = Sq + S\h 2 + . . ., and substitute into (3). Doing this 
will generate a system of equations for So, Si,... by equating powers of h 2 . 
To Zcroeth Order: 



(dSo\ 
\dx J 



2 

= 2m{E - V) 



To order ti 2 : 



2 



dS dSi _ 3 / \ 1 

dx dx A\ dSo 2 dSo 

\ dx / dx 
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Remark 3: In general, the expansion of S in powers of f? does not converge-it is 
an asymptotic expansion (initially decrease, and then increase in error), so a 
good approximation is obtained if a series is truncated at some finite values. 



The simplest approximation is truncation at zero order: 



±Jdx ^/2m{E - V(x)) + a±, E>V{ 



x 



±i / dx y/2m(V{x) — E) + /3±, E <V{x) 
where a p m and (5± are arbitrary constants. Then 

A = const • — — 



dx 



A (x) 



j±(2m(E -V(x)))~ 1/4 , E>V(x) 
5±(2m(V(x) - E))~ 1/4 , E<V{x) 



where j± and 8± are arbitrary constants. So the solution will be a linear combination 
of the plus and minus solutions: 

?/;(r) = ei w ^ with W(r) = S(r) + - In (A(r)) 

i 

gives 

• For E > V(x): 

ij){x) = (2m(E - V{x)Y 1/A (^ +e {Id*y/*riE-vm + C _ e fJ dx V^(E-V(x))^ 

(4) 

• For E < V(x): 

1/>(X) = (2m(V{x) - E)' 1 ^ ( D+e \Sd^2m{V(x)-E) + D _ g I / dx y/2m(V(x)-E) 



(5) 

Remark 4' To determine the accuracy of the solutions (4) and (5), consider the term 
0(h 2 ) in the expansion of S. Note that once we have So, we can, in principle, 
find Si. In order for this to be a good approximation, we require that 

|ft 2 Si| « |S | 

o \dV/dx\ 

=> mh 2 ! ' ! — - < 1 

\2m(E -V(x))\ 3/2 

(See Messiah Eqn VI. 47 for details) Thus, the potential must change slowly 
in space compared to the rate of change of zero approximation for the wave 
function. This does not happen when E = V (turning points in the classical 
trajectories). 
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Remark 5: The solutions (4) and (5) represent oscillatory solutions in classically 
allowed regions and exponential decay in classically forbidden regions. These 
solutions must still match at the boundaries; however, the solutions are only 
asymptotic and not necessarily well defined near the boundaries. If the form of 
the wavefunction can be determined unambiguously near the boundary, then 
it is possible to connect the two solutions. (See Messiah Ch. 6, §11 for HO for 
how to do this in this case) 



Chapter 9 

Symmetries 



In §2.1.4, we stated Nothers Theorem: symmetries in the system correspond to conservation laws. 
This holds for whether we are speaking of classical or quantum systems. Here we investigate from 
a quantum mechanical viewpoint. 

9.1 Translations, Translational Invariance, and Conservation of 
Momentum 

Classically, if the Hamiltonian is invariant with respect to x — > x + a, then p x is conserved. 

Remark 1: In quantum mechanics, we do not have a well defined position or momentum, so we 
replace these quantities by their quantum averages: 



x 



P (P) 



Remark 2: In analogy with classical mechanics, we expect translations to give: 




9.1.1 Active Translations 



Recall §1.2.4: There were 2 types of transformations: 

• the Active: transformed vectors in Hilbert space 
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• the Passive: transformed Hilbert space operators 
Here, we define the (infinitesimal) translation operator Tie) to have the following 



action: 



T(e)\iP) = \iP £ ) 



so that 



(x) —> (x) + e 



or that 



(lp e \ X \lp E ) = (lp\ x\lj))+e 

Consider the action of the translation operator on the position eigenstate: 

f (e) \x) = \x + e) 

which shifts the position to the right by amount e. Then, if (x\ip) = ip{x), we have 
(x\A) = (x\f (e)\i>) 



Remark 1: T (e) |x) = \x + e) is not the most general result. In fact, we could write: 




= S xy (y - 
= (x-e\i/j) 
= ip(x — e) 




x + e) 



space dependent 
phase factor 



This result still gives: 



(x) — > (x) + £ 



but gives 



(p) -> {x\f\e)pf{e) \x) 

= (x + e\ e - i£9{ - x)/h {-ih)^e i£9{ - x)/h \x + e) 

ax 

= (x + e\ e ~ ie ^l h ((-ih)^g'(x)e i£9( - x)/h \x + e) + {-ih)e i£g{x)l 




) 




TRANSLATIONS, TRANSLATION AL INVARIANCE, AND CONS. OF MOM. 113 



For translations, we require that 

(p) — > (p) =>■ <? = constant 
and we take this constant to be zero (for simplicity). 

Definition 1: Translational Invariance is defined by 

(i>\H\ijj) = (i, E \H\i> E ) 

i.e., the expectation of the Hamiltonian is invariant under infinitesimal trans- 
lation. 

Remark 2: According to Nother's Theorem, translational symmetry implies a con- 
servation law. (Here, it will be conservation of momentum.) 

Proposition 1: Momentum, p, is the generator of infinitesimal translations. 

• Proof: Since e is small, we expand T (e) in a Taylor series about 
e to ©(e): 

f(e) = I + eG 
where T (0) 44> no translations. Then: 

(x\ T (e) \x) = (x\ip e ) = ip(x — e) 
Expanding both sides to ©(e): 

lx\i + e6\il))=il>(x)-e^ 

ax 

tp(x) + e (x\ G \tb) = ip(x) — e^- 

dx 

in 

^f(e)=I+^-p 
in 

Proposition 2: Momentum is conserved in a translationally invariant system. 

• Proof: We begin with the definition of translational invariance: 

{ijj\H\^) = ^ e \H\^ e ) 

= ^\f\e)Hf{e)\^) 

in J \ in 
H\^)-^r pH - Hp\ifj) + 0(e 2 
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Since we are translationally invariant, 0(e 2 ) = (ip\ \p,H] \ip) = 
0. And by Ehrenfest's theorem: 

(\p,H]\ = (^p^ = =>■ momentum is conserved 



Remark 3: So far, we have considered only infinitesimal translations. If we need finite 
translations, say (x) — > (x) + a, where o > £, we divide a into N intervals of 
a/N and let N — > oo. Then we consider infinitesimal translations: 



vat; k/\r 



and define our finite translation operator as 



T(a) 

In coordinate basis: 



lim 

N^oo 



N 



N 



-iap/H 



So 



(x\f(a) 



T(a) 



{x\e~ a -c 



dip 1 2 d 2 ijj 



dx 



dx 2 



+ 



ijj(x — a) 



Remark 4 : For consecutive translations 

f (a) f (b) = e -™P/n e -ibp/K = e i(a+b)p/h = f (a + b) 



9.1.2 Passive Transformations 

Passive translations leave quantum states invariant, whil shifting the coordinate sys- 
tem to the left by the amount e (assuming position is shifted to the right in the active 
case). See Shankar for the complete derivation of the translational operator in passive 
version. Briefly, the requirement that the position and momentum operators must 
obey is: 

f ] (e)x f (e) ^x + el 
rt(e)#f(e)-p 
Expanding T (e) in powers of e results in: 



in 
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which we note is the same expression as in the active case. In the passive transforma- 
tion picture, translational invariance requires 

f ] {e)HT{e) = H 

which again requires that the Hamiltonian is invariant under translation. Now we 
want to show that this is equivalent to 



H (x + ei,p) = H(x,p) 



Proposition: For any Q(x,p) that can be expanded in a power series, and for and 
unitary operator U: 



Proof: Expand in a power series and consider a typical term 
(i.e. pxp). Then, for this term, we have 

U ] pxpU = U ] pUU ] xUU ] plJ 

Collapsing the expansion, we get: 

n (u ] xu,u ] pu y 



Corollary: 

h{x, p) -> f tftf = n (f ] xf, f ] pT] =n( x + ei,p 



• Proof: T is unitary. The rest follows from above. 
Thus we see that for £1 = H, translational invariance implies that 

H(x,p) = H(x + el,p) 



9.1.3 Translations For a Many Particle System 

For a system of N particles, the translated many body wave-function becomes: 

(xi, . . . , x N \ f (e) = ip(xt - e, x 2 - e, . . . , x N - e) 
and, to 0(e), 

(xi,...,x N \yi- = ^(xi, ...,x N )-^2 e J~r 
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. N 

i=l 

where V is the total momentum operator. 
Remark 1: For many particles, translations generalize to 

f\e)xif(e) =Xi + eI 
f\e)p l f{e)=p l 
and translational invariance becomes: 
f ] (e)H(xi,x 2 , . . .,x N ;p 1 ,p 2 , . . . ,p N )f (e) = H(x\ + el,x 2 + el, . . . , x N + el;p 1 ,p 2 , 

= H(xi,x 2 , ■ ■ .,x N ;p!,p2, . . . ,p N ) 



,Pn) 



Remark 2: For a single particle system, translational invariance implies free particle 
which implies that V = 0. For a system of particles, translational invari- 



ance implies than V = V(£i 



-so that the Hamiltonian depends only on 



interactions between particles and not on "external" interactions. 



Remark 3: In a translationally invariant system: 

f\a)Hf{a)=H and f f 

^ H = f{a)HT ] {a) 



T 



T(a),H 



T(a),U(t) 







This tells us that a system that starts translationally invariant remains trans- 
lationally invariant. 



9.2 Time Translational Invariance and Energy Conservation 

We expect that homogeneity of time will give the same result if an experiment is repeated at 
different times. 

A system is prepared in an initial state, \ipo) at time t±. The Hamiltonian at this time is H\, and 
the system evolves for a short time r. 

m 1 +T)) = U(r)\^ Q ) 

= (i- '-Hit) |^o) 



9.3. DISCRETE SYMMETRIES: PARITY AND TIME REVERSAL 
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where Hi = H{t\). The experiment is repeated, only system is prepared at time t 2 : 



m 2 + r))= [I--H 2 r\ |^o) 

where H 2 = H(t 2 ). Since the outcome does not depend on when the experiment is done, we require 
that 

hK*i+r)) = |^2 + t)> 
=> H(h) = H{t 2 ) 

This holds for arbitrary t\ and t 2 , so we can conclude that 

8H 



dt 







which implies that H has no explicit time dependence. 

Remark 1: Time translational invariance mean H is time-independent. 



Remark 2: Remember Eherenfest's Theorem: 



d 
It 


(A) = 




dA 

~bl 


Setting A = H and using that 


(#> = 


0, we have that 










= 



which implies that we have conservation of energy. 



9.3 Discrete Symmetries: Parity and Time Reversal 
9.3.1 Parity Invariance 

Definition 1; The Parity operator, II, is defined by its action on the position eigen- 
state: 

n \x) = \—x) 

And it has the following properties: 

1. n = n- 1 



Proof: 

n 2 = i ipfr 1 = in- 1 n = u- 1 



n 2 \x) = u\-x) = \x) 
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2. A has e- values dbl, which follows from IF = I and II acting on an 
eigenstate. 

Eigenvalue =1 Even parity 

Eigenvalue = — 1 =>• Odd parity 

3. f[t = fi (Hermitian), and fttft = I (Unitary) 
The action of the parity operator on an arbitrary ket is 

/oo 
\y) (vW) dy 
-oo 

/oo 
\-y) (y\tp) dy 
-oo 

POO 

= / \y) {-y\ip) dy 



where the last line is by a simple change of variables. We know that {x\ip) = 
i/)(x), and thus from the above we find that 

\x) n = ip(—x) 

and even or odd parity corresponds to the even or odd nature of the function 
ip(x). 

Remark 1: fi \p) = \p) (which follows from the definition of p in position represen- 
tation.) 

Remark 2: The parity operator is the mirror image of a function about the origin. 

Remark 3: For passive transformations 

fix = —x 
lip = —p 

and parity invar iance happens when 

&H(x,p)fl = H(-x, -p) = H(x,p) 
which means that HU(t) = U (t)TL and thus parity is preserved over time. 



Remark 4 : If weak interactions are present, II, H ^ and thus parity is not 
preserved over time. 
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Remark 5: Parity transformations are different from space and time translations 
(and rotations) since this is a discrete transformation . It may be thought of 
as a mirror image of some process. 

9.3.2 Time-Reversal 

Definition 1: The time-reversed state is defined by 

x R (t) = x(t) and p R (t) = -p(t) 
(The system is moving in reverse.) 

Time- Reversal Invariance (TRI): From an initial state of a system x(0) and p(0), 
the system evolves and at time T is at x(T) and p(T). Then we reverse the 
system and run for time T. We have TRI when 

x(2T) = x(0) and p{2T) = -p(0) 



Remark 1 : TRI exists when you can time evolve the system forward or backward in 
time and there is no violation of physical laws. 

Remark 2: Quantum mechanically, time reversal is consistent with our ijj — > ip* 
since, in position representation, 

x — > x and p — > —p 

under complex conjugation. For example, if we have an initial state ip(x,0) 
and we time evolve that state for time T: 

i/>(x,T) = e- lflT/h ifj(x,0) 
If we then time-reverse the system for time T: 

1>(x,T) 1>*(x,T) = e^* T/, V*(x,0) 
We then evolve the system forward again (for time T) 
e -iHT/n e iH*T/n r{X}0) 

Now, for TRI, we want vp(x,2T) = ip(x,0) and H = H* , the latter of which 
happen automatically whenever the Hamiltonian in real (for ever powers of p). 



Remark 3: This is another example of discrete symmetry. 
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9.4 Rotations, Rotational Invariance, and Conservation of Angu- 
lar Momentum 

9.4.1 Rotations 

9.4.1.1 Rotations in 2D 

Recall in classical mechanics, we define a counterclockwise rotation of a vector 
about the z-axis by: 




X 



X 



y 



xcosOo) - ysin(<fo 
x sin(</>o) + y cos((/>o 



And we identify a rotation matrix 





which also rotates momentum vectors. 

In QM, define the operator U[R((j)z)] as the operator which rotates a vector in 
Hilbert space: 




(x) R = {iJjr\x\i/} R ) = (x) cosOo) - (y) sin(^> ) 
(v)r = {4>R\y\i>R) = (x) sin(^o) + (y) cos(0 o ) 



where 



(x) = (tp\x\ip) 

(y) = (ip\y\il>) 
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and similar expression hold for the x and y components of momentum. 
Remark 1: Operating U[R] on position eigenkets: 

U[R] \x,y) = |xcos(V>o) - y sin(^ ), x sin(0 o ) + ycos(0o)) 

Now consider infinitesimal rotations about the z-axis: 



U[R(ez)} = 1 + 



where L z is the generator of infinitesimal rotations. Consider now this action on 
a position eigenket: 

U[R] \x, y) = \x - ye z ,y + xe z ) 

From this, we can show that 

£zT z 



x,y 



1 + 



iH 



V> ) = tp(x + ye z ,y - xe z ) 



Expanding this to 0{e z ): 



x,y 



L- 



x,y 



L, 



dip dip 



So in position representation 
Lz - 

or, more generally, 



x ( ~ in ir I ~ y (~ in ir 

ay J \ ox 



L z = xp y - yp x 



(t) 



Remark 2: Putting (f) into momentum representation and acting on ifj(p x ,p y ) 
rotates the momentum space wavefunctionso that momentum expecta- 
tions are consistant with classical rotations. 



Remark 3: The passive version of a rotational transformation has: 

U\R]xU\R] =x-ye z ' 

U^[R}yU[R}=x£ z + y 
tf[R]p x U[R]=p x -p y e z 
tf[R]p y U[R] =p y +p x e z ) 



(tt) 
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Remark 4 : Substitute the infinitesimal version of U[R] (I+zfa) into (ft) and 
get, to 0(e z ): 



x, L z 

y,L z 

Pyi L z 



-ihy 

ihx 

—ihp~y 

ihPy 



> -> L z = xp y - yp x 



Remark 5: For finite rotations, take N infinitesimal rotations of size (fro/N as 
N oo 



N 



-i(f>oL z /h 



Writing L z — > x ( — iftg^J — V {~^^) m polar coordinates yields 



L 7 — > —ih 







So 

And 

—6— 
e 8 <i>i{;(p,(j>) 



U[R(<j>oz)\ = e 



1 " 4>0^7 + 



il>(p,<f>- 4> ) 



d 1 ,,3 2 



2r°d<p 



+ ■■■) Mp, 



Remark 6: Two consecutive rotations 

U[R{^z)]U[R^z)] = U[R((<& + <Mz)] 



Remark 7: Physically, L z is the angular momentum operator, which is analo- 
gous to the classical definition, and is the generator of infinitesimal rota- 
tions about the z axis. 



Remark 8: If the system is invariant under rotations about the z axis, then 
U*[R]H(x,p x ,y,p y )U[R] = H(x,p x ,y,p y ) 
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Expanding this to order e z for infinitesimal rotations gives 

l z ,h] =0 

which implies 

1. (Lz) equals a constant, which then implies that angular mo- 
mentum is conserved. 

2. Outcomes of experiments on rotationally invariant systems will 
be the same for differing orientations of the system. 

3. L z and H can be simultaneously diagonalized, and therefore a 
common eigenbasis exists 

Remark 9: A transformation which consists of a product of translations and ro- 
tations will, in general, depend on the order of individual transformations- 
e.g. translations and rotations do not commute. 

9.4.1.2 Rotations in 3D 

The results of the previous section may be generalized to include rotations about 
all 3 coordinate axes. The corresponding rotation matrices are: 

(cos(</>) — sin((/)) 0^ 
sin(>) cos(0) 
' T 

/l 

R(cf)x) = cos(4>) — sin(</>) 
\0 sm(4>) cos(0) 

(cos((/)) sin(</>) 
1 ' 
— sm((j)) cos((/>) / 

The generalization gives the components of angular momentum: 



(T) 




Remark 1: Equatios (f) are related by cyclic permutations of indices, i.e. 

xyz — > yzx —> zxy 
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Remark 2: One can consider products of infinitesimal translations to derive 
commutation relations between components of angular momentum. We 
know that 

( x \ ( x \ ( x \ mR(s y m ( x + £ y z \ ( x \ umeM ( x - £zV \ 

y > y — £ x% K y > y i y > y + 

\zj \z + £ x y) \zj \z-Eyx) \zj \ z J 

So now if we consider the following sequence of rotations: 
R{e x x)R{e y y)R{-e x x)R{-e y y) 




R(e x x) 



R(£yi 



R(—e x x) 



R(-e y y) 




X + E y Z + E y £ x y 
y -E X Z + £ x (z + £ x y - EyX) 
K Z + E x y - EyX - E x (y- e x z) / 

x I s y z | ^ x y ^y ^y *^ I ^ X 
y + E 2 x y- E x E y x 

- E y X + £ 2 X Z + E y (x + EyZ + Ey£ x y) j 



X -\- £y£ x y ~t~ £y% £y£-xZ 1 

y — *- 



£ x £yX -\- E x y 
-2 
'V 



z+ (4 + e\)z + E x E 2 y y 



x -\- SyS x y 



y 



Which is equivalent to 

y £x&y% 



From this, we require that the quantum rotations must satisfy 
U[R(-Eyy))U[R{-E x x))U[R(E y y))U[R(E x x)] = U[R(-E x E y z)} (ft) 
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Writing 



U[R{e x x)]=i + 



ih 

£uLi, 



U[R(e y y)]=I+^ 
U[R(e z z)] = I + 



ih 

£zL z 

ih 



and expanding and matching coefficients gives 



L x , Ly 
Similarly, one can find 

Ly, L z 

L z , L x 



ihL z 

ihL x 
ihLy 



Remark 3: The commutation relations are also related by cyclic permutations 
of the indices. 

Remark 4 : Alternative ways to express these commutation relations are 

L x L = ihL 

which is legitimate since the components don't commute. 

Remark 5: Also, introducing the Levi-Civita fully antisymmetric tensor e^-fe 
where 

if i = j, j = k, k = i 
tijk = \ 1 for even permutations 
-1 for odd permutations 

and ei23 = 1, we can write 



Li , Lj 



ihf-ijkLk 



where summation over repeated indices is implied. 
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Remark 6: For finite rotations in 3D-let 6 = ^ be an arbitrary axis of rotation. 
Then Lq is the generator of infinitesimal rotations about that axis and 

N 

U[R(6)) 

JV— » 

_ e -i8-L/h 




9.4.2 Eigenvalues of Angular Momentum 
9.4.2.1 Eigenvalues of L z 

For a problem that is invariant with respect to rotations about the z-axis, 



L Z ,H 







there exists a common eigenbasis for L z and H 



To find that basis, start with the eigenvalue problem for L z : 

L z \l z )=l z \h) 
In the coordinate basis, L z — > —ih-g^: 

d 

where R(p) is an arbitrary function of p. 

Remark 1: For tpi z (p,(p) to be normalized, we require: 

! /•2-k 

pdp d(j) tf{J)i z = 1 



(t) 



Remark 2: For the normalization requirement, l z seems arbitrary. However, it 
is not since we require our angular momentum operator to be Hermitian: 



ipl 



i>2 ) = ( 4>: 



L 7 



In the coordinate basis: 



tp* \~ ih 'gl ) 4>2pdpd(j) 



00 f 27T ( d \ 

Jo WP** 

Jo ^{ ih dj)^P d P d t 

2?r ^oo p2tt 



ih / pdp i>2^\ 

10 



o Jo 







ipl I —ih— ip2pdpd<p 
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Thus, to be Hermitian we want the first term to vanish, which happens 
when if)(p,0) = ijj(p,2ir). Combining this with (f) implies that 



l z = mh where m = 0, ±1, ±2, . . . 



Remark 3: Eigenvalues of L z are discrete and m is called the magnetic quantum number. 



Remark 4 : At this point, R(p) is arbitrary. However, if we are looking for a 
simultaneous eigenbasis for L z and H , then we use the eigenfunctions of 
L z and (f) in the eigenvalue problem for H. Thus the energy eigenvalues 
and eigenfunctions determine R(p). 



Remark 5: It is convenient to introduce the following functions: 

V Z7T 



This is the normalized angular part of (f): 

f27T 

$ m (0) d(j) = 5 nm 

o 



9.4.2.2 Solutions to Problems with Azimuthal Symmetry 



Once we have determined the eigenfunctions of angular momentum, we can pro- 
ceed to solve the eigenvalue equation for H. 



L Z ,H 







H is not an explicit function of 



=>vM) = v(p) 

In cylindrical coordinates, the Schrodinger equation is 



2p \d P 2 + P dp + p 2 d(p ' " 



Remark 1: p here is equal to the mass, and is used so as to not be confused 
with the magnetic quantum number. 

Trying a solution of the form: 
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which is an eigenfunction of L z and <3? r 
Schrodinger equation: 



_1 im< 

2^ 



e* m ^. Substituting this into the 



-h 2 ( d 2 Id 
+ 



m 



2 fx \dp 2 pdp p 2 



+ V(p) 



REm{p) = ER Em (p) 



Remark 2: The azimuthally symmetric problem is reduced to solving a 1-D ra- 
dial problem. The angular contribution is an effective repulsive potential 
that corresponds to the centrifugal force. 



9.4.2.3 Eigenvalues of L 2 and L z 



Define L 2 = L 2 + L 2 + L 2 . One can show that 



L 2 , Lj 



for all i = x,y, z 
Remark 1: We already have shown that if the Hamiltonian is invariant under 



rotations about the z-axis, then then 



H,L Z 



and L z is thus conserved. It 



follows that if H is invariant under arbitrary rotations, then 



H, Li 



and 



for all i = x,y, z 



H, Lf 



implies that L 2 , L x , L y , and L z are conserved. 



Remark 2: Since the components of the angular momentum vector do not 
commute, we can not construct a common eigenbasis between H,L 2 ,Li 
for all i = x,y, z. So we pick a single component, usually L z , and form a 
common basis with H and L 2 . 

Now to find a common eigenbasis between L z and L 2 . Let \af3) be a common 
eigenvector so that 

L 2 | a/3) = a \ ac0) 
L z \a(3) =P\af3) 



Remark 3: Recall the HO problem where we defined the annihilation and cre- 
ation operators in terms of x,p. Here we define L+,L_ in terms of L x 
and L v . 
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Definition 1: Let 



L± = L x i iLy 
be the raising and lowering operators. We note that 



±hL± 



L 2 ,L± =0 



Remark 4 : L± raise/lower the eigenvalue of L z while leaving the eigenvalues of 
L 2 fixed. 



and 



L z [L ± \af3)) = [L±L Z ± HL±j |a/3) 
= L± (L z ± ft) | a/3) 

= L± {15 ±h) | a/3) 
= (/3±ft) L± |q/3) 



L 2 (L ± |a/3)J = L± [L 2 |a/3) 
= a ( L± | a/3) 



Remark 5: Since L± raises/lowers the eigenvalue of the z-component of the 
angular momentum, we deduce that 

L± | a/3) = C±(a,/3) \a,P±h) 

However, given a state |a/3), we can not raise or lower (3 an arbitrary 
amount as \l z \ < V^ 2 . This implies 



a/3 



L 2 -L 2 



a/3 ) = a — /3 2 
a/3 



a/9) > 



a > /3 2 



Remark 6: Since /3 2 is bounded by a, there must exist a state \a(3 max ) that can 
not be raised any further. 



L+ \a/3„ 
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Similarly, there must be a \aj3 m i n ) that can not be lowered. 

L_ \a(3 min ) = 



Starting with the "top state" and considering the above, we have that 

L_L + \a(3 max ) = 

= (t, x - iLy^j [L x + iLy^J \a(3 max ) 



L 2 x + L 2 + i \ L x L y - L y L x ) ) \a/3 v 
L 2 - L 2 Z - hL^j \af3max) 

(oi @max foPmax^) l^Pmax / 
=^ Ot = ftmax (ftrnax ~\~ K) 

Similarly, starting with the "bottom state" : 
L + L_ \af3 min ) = 

L 2 - L 2 + hL z ) \a(5 m 



OL — f3 m in \Pmin h) 

Comparing values of a implies that P m in = ~Pmax- 

Remark 7: One can also start at the top state and use L_ to get to the bottom 
(or use L + to go from \a(5 m i n ) to \a(3 max )). In general, this will take k 
steps, where k = 0, 1, 2, . . ., so that 

Pmax Pmin = "^Pmax = hk 

hk 

=> Pmax = ~y where k = 0, 1, 2, . . . 



Remark 8: Let ^ m ^ x = -| denote the angular momentum of the state. 

Remark 9: Note that in §9.4.2.1, we found the eigenvalues of L z to be rah (where 
m is an integer), while here we find that the eigenvalues can be half- 
integers. The reason is that here we used only commutation relations and 
not the specific form of L z (=> —ih-jj^). So, in fact, the result here is valid 
for any angular momentum with the same commutation relations without 
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specifying the type of angular momentum. We will study two types of 

angular momentum: orbital angular momentum (denoted by L), which 
has only integer eigenvalues, and spin angular momentum (denoted by 

S), which can have integer or half-integer eigenvalues. The total angular 

momentum, J = L + S , obeys the same commutation relations that were 

derived for L, and hence the eigenvalues derived above hold for J. For 
total angular momentum then: 



j 2 \jm) = j(j + l)h 2 \jm) , where j = -, 1, -, . . . 



Jz \jm) 



mh \jm) , where m = j, j — 1, . . . , — j + 1, j 



For J = L: 

L 2 \lm) = 1(1 + l)h 2 \lm) , where / = 0, 1, 2, . . . 
L z \lm) = mh\lm) , where m = I, I — 1, . . . , — I + 1, — I 

Remark 10: We will discuss rules for addition of angular momentum shortly 
- for now, we will just use the notion of total angular momentum and 
determine general eigenfunctions. 

9.4.3 Matrix Representation of Angular Momentum 

Finding the matirx representation of J 2 , J x , J y , J z is sufficient for determining the 
eigenvectors of angular momentum. We will look at the coordinate representation of 
these eigenvectors when we undertake central potential problems. 

9.4.3.1 Matrix Elements of J 2 and J z 



These are both diagonal and therefore straight forward to compute. For j 2 
have: 



we 




where j 






2h 2 



2h 2 



2h 2 
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And for j, we have: 



jm 



jm) = mh 



J, 



\0 



'■■/ 



.4.3.2 Matrix Elements of J x and jj, 

These are non-diagonal. We will use 



Jx = 2 ( 



1 



and J y = — [ J+ — J. 



and find the action of j± on the state \jm). We know: 

J ± \j m ) = C±(j, m) \j, mil) 
where we need to determine C±(j,m). In dual space, we have 

(jm\ J| = (jm| J T = C± (j, m ± 1| 
Consider just C+ for now. We know that 



J_J 4 



J7TI 



|C+| 2 0',m+ l\j,m+ 1) = |C 4 



and 



jm 



J 2 - Jl - HJ Z 



jm 



\C 4 



\C + \ 2 = j(j + l)h 2 - m 2 h 2 - mh 2 
= H 2 (j — m)(j + m+l) 



C+ = hy/{j-m)(j + m + \) 



Remark 1: In general, C+ can have an overall phase, which we choose to be 
unity. 



Similarly, considering C_ 



jm 



|C_ 
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We find that 



C_ = hy/(J +m)(J -m+ 1) 



=> J± \ jm) = %\J (j T m)(j ± m + 1) | j, mil) 
Using this to determine the matrix elements of J x and J y : 



jm 



, ■> i 
2 U m 



J+ + J_ 



- ^j'W ft\/ (j - m)(j + m + 1) j, m + 1^ + - ^j'W ft a/ (j + m)(j - m + 1) 



-ft 
2 



Thus 



V U - m)(j + m + T)5jj'6 m +i,m> + V (j + - m + l)^/^-!,, 



/...P. 



- 

u 2 

- 

2 u 



V o . 

Similarly, for J y we have 






v/2 





h 





/< 









V2 






.P..\ 



•3 (i 



J+ - J_ 



jm 



- 1 
~ 2i 

Thus 



V (i -m)(j + m + 1)6 jt - V + m)0 -m + l)<Jj )3 -'5 m -i,m' 



1/77 



V o I 



i n - 

; u 2i 

i =h n 

; 2i u 






773 
— h q ?i 

173 



.Q..\ 



Remark 2: Though J x and J y are not diagonal, they are block diagonal (all 
matrix elements have 5jj> ) . This implies that the blocks do not mix upon 
multiplication, and we can conclude that: 



fU) tU) 



ihj[ j) for 



o,ii,.. 
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.4.4 Finite Rotations 

Recall from §9.4.1.2 Remark 6 that finite rotations may be expressed as 

U[R(0)] = e ~ i8J/h 

Remark 1: We have all components of J; however, J is block diagonal with only 
blocks of the same j multiplying. So, schematically, we represent U[R] by the 
following: 



U[R]- 



\jm) basis 



D {0) 






\ 




£)(l/2) 










D (l) 











'•• ) 



where D^' is the (2j + 1) dimensional block for a given j. 

Remark 2: Since a vector \ipj) need only be spanned by 2j+l vectors \ jj) , . . . , \ j, —j) 
we only need D^[R] to rotate (which stays within the subspace Vj). 

Consider 



Y- 

n=0 



(T) 



Remark 3: This series expansion can be truncated at n = 2j. 



2j 



D^[R]=Y,fn{0) (9- 



■(j) 



n=0 



The terms (j) ■ j^J for n > 2j can be written as linear combinations of the 
first 2j terms, and the f n {9) represents that combination. 



Remark 4 : Consider the subspace Yj spanned by basis vectors \jm) where m = 
j, . . . , — j. The subspace is identified by the eigenvalue of J 2 : j(j + l)h 2 . Since 



J 2 ,U[R] 







then a rotation of any state spanned in this basis will not change the eigenvalues 
and therefore this is called an invariant subspace (invariant under arbitrary 
rotations). 
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Remark 5: Furthermore, the invariant subspaces are irreducible, i.e., it does not 
contain any other invariant subspace. 

Remark 6: The block diagonal matrices, which make up U[R] are an irreducible 
representation of rotation in the \jm) basis. 



Remark 7: Also, since 



H,J± 



energy eigenvalue Ej . T 



= 0, each block represented by will have a single 
lis implies that all states with a given j are degenerate 



in the rotationally invariant problem. 



Remark 8: Classically, this degeneracy is because one state can be rotated to an- 
other without changing the energy. Quantum mechanically, a rotation com- 
bined with J± is required to take one state to another, but these operators all 
commute with H, so the energy is unchanged under U[R] and J±. 
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Chapter 10 

Central Potential 



10.1 Hamiltonian in Spherical Coordinates 

For a particle with mass \x and momentum p, the Hamiltonian in a spherically symmetric potential, 
V(r), is 

H=?- + V(r) 

which corresponds to a time independent Schrodinger equation: 



Hift(r) 



■0(r) = Eift{r) 



(1) 



In spherical coordinates, this is: 



1 



sm(0)^ + 



t) 2 



2fi \r 2 dr 8r r 2 sin(#) 86 89 r 2 sin 2 (9) 8 '(ft 2 



+ V(r) 



4> E (r,9,<j)) = E4>(r,9,, 



(2) 



Remark 1: Though we can tackle (2) as is, it is useful at this point to note the spherical symmetry 
and thus that we have conservation of angular momentum. Consequently, we can express 
(1) in terms of a conserved quantity. 

Remark 2: Recall that L = fxp = ih(rxV) = angular momentum, and that all three components 
of angular momentum commute with the Hamiltonian. 



Remark 3: For the radial momentum, we will not use —ifr-§p 



p since it is not Hermitian. 



Instead, we will use the symmetrized form: 

1 



Pr 



- -P + P-- 

2 \r r 



Remark 4 : Pr commutes with any function of 9 and eft, and also with l x ,l y ,l z , but 

[f,p r ] = ifl 
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Proposition: 

L 2 



p =p r H 2 for r / 



Proof: 



L 2 = (f x p) • (f x p) 

= r 2 p 2 — r 2 p 2 by vector identities 



So an alternate form of the Hamiltonian for a spherically symmetric system is: 

Vl , L 2 



H 



2fi + 



+ V(r) 



Remark 5: Comparing the Schrodinger equation obtained by (3) 



« 2 ^ 2 



with (2) gives us 



-h 2 



sin 2 (9) 



d ( d\ d 2 



(3) 



(4) 



(5) 



10.2 Solution of Spherically Symmetric Schrodinger Equation 

Recall from Chapter 9 that spherical symmetry implies 
• H is invariant under transformation 



H, Li 
H, L 2 



for i 
■ 



x,y,z 



• There exists an eigenbasis common to H,L 2 and one Li, usually L z . 
The steps to finding a common eigenbasis are: 

1. Separate the angular and radial parts of the Schrodinger equation. 

2. Solve the eigenvalue problem of L z and L 2 for the angular part. 

3. Solve the full problem to get the radial part. 



Remark 1: The particular form of V(r) enters only in Step 3! 
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10.2.1 Solution of the Angular Part 

We've actually done this using the bra-ket notation in previous chapters using: 

H |n, I, m) = E\ m \n, I, m) 

and starting by solving the eigenvalue equation for I? and L z . Again, we will solve 
the eigenvalue equation for L 2 and L z , but this time we'll do it in coordinate repre- 
sentation. Following the same strategy as before, we'll start with 

L+ \l,l) = 

In coordinate representation: 

i ± = ± ^(| ±i cot W |) 
We let \l, I) — ► ip\(r, 9, <ft) so that we have: 

(|^<crt(0)A)v*(r,M)=O (t) 



Remark 1: ipl is an eigenfunction of L z with eigenvalue lh: 

L 7 = —ih—— 



So let 

4^,9,0) = Ul(r,e)e il <t> 
Substituting this into (f) gives: 



(j^+cot(0)zWr,0) = O (t) 



Remark 2: We can ignore the r dependence in (J), and we can write 

Uj sm.{9) 8in(0) 
which is satisfied if 

U}(r,6) = R{r) {sin{9)) 1 



Remark 3: R{r) is normalizable, but otherwise arbitrary-in actuality, it is deter- 
mined by the radial part of the solution. 
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Remark 4 : Requiring R(r) to be normalizable with respect to r implies that the 
angular part must be normalizable with respect to the angles. A function that 
satisfies this requirement is: 



Y 



i \ u -,< 



21 + 1 
2vr 



1/2 



1 



(sin(#)) l e 



l il<j> 



(ft) 



Remark 5: Applying the lowering operator to this top state 

L-\l,i) = h[{i + i)(i)] 1/2 \i,i-i) 
= hV2i\i,i- 1) 



Y 



i-i, 



— H (~ 



icom^Yi 



Repeating this gives the spherical harmonics. For m > they are 



"(2/ + 1)!" 


1/2 x - 


4tt 


2H\ 



(l + m)\ 1 1/2 



(2l)\(l-m)\ 



s- 



<i(cos 9) 



(sinS) 

l—m v ' 



2/ 



And for m < 0: 



Y~ m = {-l) m (Y^)* 
These function satisfy the orthonormality condition: 



(9,<t>) dn = 6 l>l ,6 ri 



with d£l = sin 9 



Remark 6: The Spherical-Harmonics are closely related to the Associated Legendre 
Polynomials. For P™ with < m < I, they are: 



Y t m (0A) 



(2l + l)(l-m)\ 



n i/2 



4tt(1 + m)\ 



(-l) m e im< ^P z m (cosl 



10.2.2 Solution to the Radial Part 

Since the solutions will be simultaneous eigenfunctions of 1? and H , we assume a 
solution of the form 



4>nUr,9,cj>) = R El (r)Y l m (9, l 



1 0. 2. SOL UTION OF SPHERICALLY SYMMETRIC SCHRODINGER EQ UATION 



141 



Substituting this into the spherically symmetric Schrodinger equation: 

S 2 12 



Pr 



+ 



2/i 2\ir 2 

f r l(l + i)h 2 



+ 



2/z 2fir 2 
p 2 r 1(1 + i)h 2 



+ V{r) 
+ V(r) 



Y l m (9,<f>)R(r)=E lm Y l m (e,cl>)R(r) 
Yr{e,cj>)R{r)=E lm Yr{e,4>)R{r) 



Since 



2fi 2/iH 



Pr 



+ V(r) - E 



R E l{r)=0 



, t i9 . t / a i 

-ih-—r = -ih — + - 
r or \orr 



( 



Pr 



d 2 id_ 



_a i 

dr r 



+ - 



V 



1 ld_ 

1*2 r dr 



Or 



2d_ 

r dr 



h 2 ( d 2 d d 
r— + — + — 



di 



n 2 



r 

r dr 2 



d_ f d_ 

dr \ dr 



dr dr 
d 



+ 



dr 



Thus we have that 

-h 2 d 2 



1(1 + l)H 2 x y. , 

r + — tt— + V(r) 



2/i dr 2 2/ir 2 
To solve this, let Uei(t) = rREi(r): 



h 2 d 2 1(1 + 1)H 2 V(r) 



E 



R E i(r) = 



E 



2jir dr 2 



2^ 



h 2 d 2 i(i + i)h 2 . 

2fidr 2+ 2/ir 2 + {r) 



E 



U(r) 
U(r) 



(t) 



Remark 1: (f) is called the radial equation and closely resembles the 1-D Schrodinger 
equation. The differences are: 

• < r < oo instead of — oo < x < oo 
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• For I 7^ 0, there exists a repulsive centrifugal barrier ■ So the 



total effective potential is 

1(1 + i)h 2 



Veff = V(r) + 



2 fir 



2 



• Boundary conditions on U(r) are different. 

Remark 2: As in the 1-D Schrodinger equation, the asymptotic solutions (ie, the 
solutions as r — > oo) are exponentials which decay to (E < 0) for fixed 
discrete values of E (bound state solutions), or they are oscillatory solutions 
(E > 0) which are the continuum or unbound states. 



Remark 3: The solutions are generally degenerate with respect to the energy eigen- 
value: 

• Energies with discrete spectrum are degenerate with respect to quan- 
tum number m, and may be, but not necessarily, degenerate with re- 
spect to the angular momentum quantum number /. 

• Continuous energy spectrum are always infinitely degenerate, since for 
each E > there are eigenfunctions for all values of / = 0,1,2,... 
and m = 0, ±1, ±2, . . . , ±Z (assuming that V(r) — > monotonically as 
r — ► oo). 



Remark 4 : We must determine boundary conditions at r = so that the solutions 
are physically acceptable. First of all, we require that 



. . h 2 d 2 l(l + i)h 2 , . 



to be Hermitian with respect to functions Uei- Here, 

D l (r)U E i = EU El 



is the eigenvalue equation representing the radial equation, (f). 
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Remark 5: Requiring that D{{r) is Hermitian is equivalent to requiring that p r be 
Hermitian (Convince yourself of this with Exercise 12.6.3 in Shankar). 



\Pr 



f°° I 1 d \ 

dQ. I r 2 drip* I —ih——r I ij) 
Jo \ rdr ) 

f°° ( d \ 

dft / rdr ib* I —ih—r ib 
Jo \ dr J 

™drrYr*R* El (-ih^rYTR 
dQYrrr J°° dr (-ih)rR* El ^- (rR El ) 



dn 



-in \tRei\ Iq + z " 



dr 



Tr {rRh) 



rR 



El 



/roo 

-ih \rR E i\ 2 \^ + (ip\p r \ip)* 



ITr {rRh) 



rREl 



So p r is Hermitian if 



-ih \rR El \ 



= \Up 







io - i^l lo 

Since we require ib to be normalizable, this does vanish for sure at the upper 
limit. For physical solutions then, we have the requirement that 

U m {r = 0) = 



Our goal now is to understand how solutions to the radial equation behave at the 
origin. To investigate this, let's try an expansion of the form 

U E i(r -> 0) ~ r s (1 + air + a 2 r 2 + . . .) 

where s, ai, 02, • • • are constants and s > 0. 

We will assume that if V(r) — > 00 as r — ► 0, then it does so no faster than - (the 
Coulomb potential). Let's substitute this into the radial equation: 

'-h 2 d 2 i{i + \)h 2 



2fi dr 2 2/j,r 



+ V(r) - E 

-2 



* [ s(s _ iy-z + air .-i (a + i )s ] + 

Keeping only the dominant terms leaves 



r s + ai r s+1 + a 2 r s+2 + ...)= 
+ air° _i ) + y(r)r s - Er s = 



'2/i 



s(s - 1) - i(Z + 1) 



radial 
KE 



Centrifugal 
KE 



„s-2 
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^s(s-l) =1(1 + 1) 
s = I + 1 or s = —I 

So, near the origin: 

U El r ~*°> r l+1 which is the regular solution. 

or 

Uei ► r which is the irregular solution 

Remark 6: The irregular solution is rejected since it is not normalizable, and it does 
not satisfy the boundary conditions Uei(0) = 0. 

Remark 7: It is clear that the irregular solution must be discarded for I ^ 0, but 
it also must be discarded for 1 = 0. If I = 0, then the corresponding wave 
function: 

V>o = Y°R E0 = Y^Ueo ~ \ 
which doesn't satisfy the Schrodinger equation: 

(H -E)^ = 2 -^5{r) + 
since V 2 ^ = &it5(r). 



10.2.3 Examples: Free Particles and Central Square Well Potentials 

There are several "classic" problems that correspond to central potential problems: 

• Free particle in spherical coordinates: V(r) = Vq = constant 

• Isotropic H.O: F(r)/i / uu; 2 r 2 

2 

• Hydrogen atom: V(r) = 

• 3D versions of square wells and step potentials 

10.2.3.1 Free Particles and Spherical Bessel Functions 

Consider a free particle such that V(r) = Vq for all r > 0. Then the solution to 
the spherically symmetric Schrodinger equation is: 

1>Eim(r,e,<l>)=REi(r)Yl n (0,<l>) 

And the radial equation is: 



ft 2 d 2 | Z(Z + l)ft 2 | Vo E 
2/jl dr 2 2fir 2 



U E i{i 
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with Uei{t) = rRzi{r) (See §10.2.2). Consider now the case where E > V$. Let 
k 2 = ME-v ) 



Or, in terms of Rei- 



d 2 . , 2 1(1 + 1) 



dr 2 



+ V 



U El = 



d 2 | k 2 1(1+3 
dr 2 r 2 



tRei = 



£_ 2d_ k<2 
dr 2 r dr 



d 2 



+ 



d 



d(kr) 2 (kr) d(kr) 



+ 1 



1(1 + 1) 

1(1 + 1) 

(kr) 2 



R 



El 



Rei = 



(t) 



Remark 1: (f) is a differential equation whose solutions are spherical Bessel 
functions. 

ji(kr) — > The proper spherical Bessel functions 
ni(kr) — > The Neumann functions 

h^ + \kr) — > The Hankel functions of the 1st kind 
hy\kr) — > The Hankel functions of the 2nd kind 



Remark 2: The general solution is a linear combination of ji(kr) and ni(kr) 
which are defined in terms of ordinary Bessel functions: 

/ 7T \l/2 

ni(kr) = (-!)< (^) 1/2 J_,_ 1/2 (*r) 



Remark 3: Near the origin: 



n,(A;r) (Ax)-( z+1 ) 



regular 
=>• irregular 



Remark 4 : The Hankel functions are given by 



h 



(±) 
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which are both irregular at the origin. 

Remark 5: For physical solutions, only regular solutions are allowed in regions 
containing the origin. However, the other functions may be important 
for matching discontinuities in the potential (Eg. a central square well 
problem) . 

Remark 6: Asymptotically, 

. , fcr^oo 1 . f, In 

7/ (kr) ► — sin kr 

J v ' kr \ 2 

/ , \ kr—>oo 1 / , ^ 

Tii(kr) ► — — cos kr 

kr V 2 



and 



h^fkr) J_ e ±i(*r-J W /2) 

kr 



Remark 7: Asymptotically, the factor of ^- in ji(kr), ni(kr), and (kr) gives 
the inverse square law behavior expected in a central potential problem. 

Remark 8: To lowest order, the spherical Bessel functions are: 

,, . s'm(kr) ,, . sin(fer) cos(kr) 

Jo{kr) = — : , ji{kr) - 



kr (kr) 2 kr 

. cos(kr) , cos(kr) sin(kr) 



Remark 9: The solution to the spherical Schrodinger equation with V = Vq and 
E>V is 

4>Elm{r,e,<f ) )=ji{kr)Y l m {eA) (t) 
where k 2 = 2 ^(^~^° and the solutions satisfy 

where 

jl(kr)ji(k'r)r 2 dr = ^p-<K fc ~ 
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Remark 10: For V(r) = Vq and E < Vq, we simply let 



h 2 V h 2 



bolic functions. In this case, it turns out that hi is the only solution 



so that the analogs in trigonometric functions become analogs of hyper- 
bolic functions. In this case, it 
that doesn't diverge at infinity. 



-Kr 



Kr 

whereas all other functions diverge as 



e Kr 
Kr 



Remark 11: In Cartesian coordinates, the solution for the free particle is given 
in terms of plane waves: 



which are eigenfunctions of the energy ^ and the components of linear 
momentum p x ,Py,Pz- These plane waves form a complete set of functions 
which may be used to expand any wavepacket. The spherical standing 
wave (f) in Remark 9 is an alternate complete set. 

We can expand the plane waves in terms of the spherical standing waves (J): 

oo I 



e ifeV = 47r£ J- i l Mkr)Yr (£) *T (?) 

1=0 m=-l \ ' 



where the two terms in parenthesis denote angles 9 and <j> which define the unit 
vectors in the k and r directions. If k is along the z-axis, this simplifies to 



oo 



e ik-f = e ikr cos e Y^(2l + l)i l j,(Jfer)Jl(cose) 



1=0 

where Pi(cos6) are the Legendre polynomials. In this case, a particle moving along 
the z-axis has no angular momentum in this direction. Thus, Y™ — ► YJ° and the m in 
the sum in completed explicitly. 

10.2.3.2 Central Square Well 

Consider now the central potential defined by 



V(r) 



V , < r < a 
0, r > a 
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V(r) 



For E < (Bound case):: 

Aji(kr), r < a regular at origin 
r > a finite at infinity 
Solutions must match at r = a. 

(+)/ 



Rei{t) 



Aji(ka) = Btif~\iaK.) 
Akj[(ka) = Bh\ + ^ (ian)in 

Recursion relations relate derivatives of spherical Bessel functions to other 
spherical Bessel functions. Let us consider the I = case. Then: 



As'm(ka) 
ka 



B 



B- 



cos(ina) s'm(iKa) 
— r +^ — r 



i(iKa) 



B- 



(The B absorbed the i) 



(1) 



and 



Ak 



cos(ka) sin(£;a) 



ka 

A ( cos(fca) 



(ka) 2 
sin(A:a) 



ka 



Bin 



B 



ic. 



KG 



i(Ka) 5 



HO 



(2) 



Adding equations (1) and (2) yields: 
Acos(A;a) = -Be~ Ka 

= —Asin(ka) 

tan(/ca) = 



by (i) 



Solving graphically: We found the bottom line by finding k in terms of k 
and plugging into — ^: 




10.2.3.3 The Coulomb Potential: Hydrogen and H-like Atoms (Ions) 



Consider the 2 particle coulomb interaction between an electron and a proton, as 
in a Hydrogen atom. In the center of mass frame, this becomes a central potential 
problem, with a Schrodinger equation: 



Hip 



1L 
2m 



V 2 + 



ip = Eip 



where m = ™ i 1 ^£ 2 is the reduced mass and z\ = +l(protons) and z^ = —l(electrons). 
If one neglects nuclear structure, then this equation also describes other 1 electron 
atoms such as He + , Li 2+ (z\ = +2, +3, . . .); it also approximates atoms with a 
single electron outside one or more complete shells. 

Recall the Bohr model, in which we toake the electron to be distributed uniformly 
over a spherical surface of radius tq. In this case, the mean potential energy is 
given by 

V(r) 



>•() 



With this model, we can assume the particle is confined on the scale of ro, and 
by the Heisenberg uncertainty principle, the uncertainty in momentum is at least 



h 

ro 



Letting this be the mean value of momentum, the total energy in the ground 
state is: 

1 ShV e 2 



E„ 

2m \tq / ro 
The stable configuration is where this is at a minimum: 



dr 



2h 2 e 2 
+ -. 



2mrj] 
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h 2 

^r = 

And plugging this into the ground state energy gives: 

h 2 fme 2 \ 2 2 ( me 2 \ 1 me 4 
a = 2^\~¥~) ~ 6 \ W) = ~2l^ 

Remark 1: ro is called the Bohr Radius . 



Remark 2: The simple arguments actually give the exact results calculated 
from the Schrodinger equation! However, in order to get the excited state 
energies, we must actually solve the Schrodinger equation. 

Thus, we have a classic central potential problem, so the angular and radial parts 
separate and the main problem is solving the radial equation. The angular part 
is the same as always, for the radial we have: 



d 2 2m / | e 2 
dr 2 h 2 V r 



1(1 + 1) 



U E i(r) = 



We are solving for the bound states, and hence are looking in the region where 




E < 0. We begin by making the change of variables p = r 
from the free particle solution). 



2mE 



(i.e., p = kr 



dp 2 



1 + 



Let 



-E 



-E 



Now, we know that 



Uei 
Uei 



2mE 



1(1 + 1) 



U E i(i 



2mE 



n 2 



—2m 

We 



as p — > 
as p — > oo 
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Hence, we will let 



U E i = p l+l e- p v El ( P ) 



to explicitly capture the asymptotic behavior. Substituting this into the radial 
equation: 



d 2 ( ,^ _„ , ,\ d r 



dp* 



[ P l+1 e-"v El {p) 



dp 



(I + l)p l e-Pv El - p l+1 e~ p V El + p l+1 e-Pv' El 

(I + \)lp l - x e~ p v El - (I + l)p l e- p v El + p\l + l)e- p v' El 
- p l (l + l)e-Pv El + p l+1 e~ p v El - p l+1 e-Pv' El 
+ p l (l + l)e- p v' El - p l+l e~ p v' El + p l+1 e- p v' El 



1(1 + l)^- 1 - 2(1 + l)p l + p l+1 ) v EX + {(I + 1)P 1 - P l+1 ) 2v' El + p l+1 v' El 
And the radial equation becomes: 

p 1+1 ^ + ((/ + Dp 1 - p 1+1 ) 2^ + (V +1 - 2(1 + 1) + i(i + i)p 1 -) + (-i + ^ - ^) p 1 ^} 

{p^ + 2 (Z + 1 - p) ^ + (-2(1 + 1) + Ae 2 ) }^ = 



(t) 



A;=0 



We will try solutions of the form: 

oo 

VEl(p) =Y,°kp k 
k=0 

Substituting this into (f): 
{Hk ~ l)C k pk - 1 + 2 (I + 1 - p) kC k p k ~ x + (-2(1 + 1) + Ae 2 ) C k p k \ = 

oo 

{[-2(1 + 1) + Ae 2 - 2k\C k p k + [(k(k - 1) + 2k(l + l)]^" 1 } = 

k=0 

oo oo 

Y^[-2(l + 1) + Ae 2 - 2k]C k p k + ^T[k(k ~ 1) + 2fc(Z + l)]C k p k ^ = 

fc=0 fe=l 
oo oo 

Y^[-2(l + 1) + Ae 2 - 2k\C k p k + ^[(fc + l)k + 2(* + + l)]C k+lP k = 

:0 k = 

oo 

{ [-2(1 + 1) + Ae 2 - 2k]C k + [(k + l)k + 2(fc + l)(l + l)]C k+1 ) p k = 



(t) 



fc=0 
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a 



fc+1 



-2(1 + 1) + Ae 2 - 2k 



-Xe 2 + 2(l + l + k) 



C k (k + l)k + 2(k + l)(l + l) (k + l)(2(l + l) + k) 

And hence we see that we get a recursion relation between the successive coef- 
ficients in the polynomial expansion (so long as we are assuming a power series 
solution for R(r)). Consider now again the series solution (f). To see how this 
behaves for large p, consider the form of the radial equation given earlier in this 
section: 

" d 2 Xe 2 1(1 + iy 



dp 2 



1 



P 



P' 



u El = o 



For large p, try a solution of the form: 

U E i ~ P a e bp 



2 Uei ~ — ( ap a - l e bp + bp a e bp 
dp V 



dp 2JEl 



~ a(a - l)p a - 2 e bp + 2abp a - 1 e bp + b 2 p a e bp 
~ [a(a - l)p a " 2 + 2a&p a - 1 + 6 2 p a ] e 6p 

Keeping the terms up to order p a ~ l in (ft): 



Xe 2 



2abp a - 1 + b 2 p a - p a + Ae 2 p a ~ 1 



6 = ±1 a = =F- 
So, Uei ^ p~ Xe2 ^ 2 e p dominates in the p — > oo limit. Now, 
U E i = f} +1 e- p v E i(p) 

with 

oo 



fc=0 



But, we also have: 



(tt) 



fc=0 



P ^ 00 , /9 -Z-l-Ae 2 /2 e 2p 



fc=0 



In order to avoid a divergent solution, we can force the series to terminate at 
finite k by requiring the coefficient to vanish at some k. This can happen if we 
set the numerator in the recursion relation for successive coefficients to 0. 



=> e 2 X + 2(k + I + 1) = 
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Since we have that A = w -^gr (where m = mass), then we can conclude: 

_ 2m me 4 
^ E = -WV = - 2V(k + l + iy for fc, i = 0, 1, 2, . . . 

Remark 3: We define n = k + I + 1 to be the principle quantum number, and 

the allowed energies. 

Remark 4-' Given n = 1, 2, 3, . . .: 

I = n — k — 1 = n — l,n — 2,. ..,1,0 (The ang. mom. Q #'s) 

are allowed values for /, and there in degeneracy for different I at each n: 

J^Z = n_1 (2Z + 1) = 2 • -(n- l)(n) + n = n 2 (# degeneracy) 



Remark 5: A Rydberg (.Ry) is a unit of energy that measures the energy levels 
of Hydrogen: 

_ me^ _ 1 



Remark 6: Equation (f), 

^ + 2{l + l ~ p) T P + (~ 2(/ + + Xe ^ } VEl = 

is the Laplace Equation, which has one solution that is regular at the 
origin. That solution can be expressed as a confluent hypergeometric 
series: 

+ 9 x 2 97 ^ 9 ^ ^ r(2Z + 2 + fc-Ae 2 ) (2Z + 1)! p k 
F 1 (2/ + 2-Ae,2Z + 2;p) = ^ ^ _ (2f + 1 + fc) , ¥ 

fc=0 

where 



/•OO 

r(z) = / t z ~ x e-Ht 

Jo 

r(n + 1) = n! for integer, positive n 
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Remark 7: The final wave functions are obtained from the recursion relation 
for the coefficients. Note that the series: 

n-l-l 
VEl(p)= Yl Ck P k 



k=0 

terminates at n — I — 1. The recursion relation can be solved: 

-.2 



C k+1 2(1 + 1 + A) - \e 2 



C k (k + l)(2l + 2 + k) 
2(l + l + k- n)C k 



Ck+l (k + l)(2l + 2 + k) 

where we used that fact that Ae 2 = 2n since n = 1 + 1 + k. Letting Cq = 1 
implies: 

2(1 + 1 - n 



21 + 2 



Thus, 



_ 2(l + 2-n) 
° 2 ~ 2(21 + 3) Cl 

_ 2 k (I + 1 - n)(l + 2 - n) ■ ■ ■ (I + k - n) 
k ~ ~k\ (2l + 2)(2l + 3)---(2l + k + l) 
- 2 1( 1)k (n-l-k)\ (21 + 1)1 
k\ ' (n-1)! (2l + k + l)\ 



1 (-l) k (n - I - k)\ (21 + 1)1 , n . h 



VEI{P) ~ ^ —\ („-!)! (21 + k + 1)! {2P) 

which in the Laguerre Polynomial : L 2 '+ 1 _ 1 (2p). 

Remark 8: Finally, we write the normalized wavefunction for the Hydrogen 
atom: 



_2_ l(n-l-l)\^ 3/2 (2r 

'rnlm 9 



n 2 \ 3(n + l)! 
with 



£-) (2 P ) l e-^U2p)Y l 



-2mE r /2m 2 e 4 rme 1 r 
K 2 H\ 2h 2 n 2 nti 2 n ro 

where 



h 2 

ro = 9 

is the Bohr radius 



Chapter 11 

Angular Momentum: Spin 



Spin in another concept that completes the angular momentum picture. Loosely, spin is imagined as 
angular momentum associated with a particle rotating about its axis (like Earth). However, there 
really is no such mechanical analogy. Nevertheless, spin represents an intrinsic angular momentum 
of a particle that is not associated with the angular momentum operator L. 

Principle Evidence of Spin: The principle evidence of spin is in the behavior of electrons in 
a magnetic field, where the degeneracy is lifted (Zeemann effect). For example, for 1 = 1: 
where H is the external magnetic field. Spin is an angular momentum quantity and thus 



it's corresponding operator behaves much the same as L (See §9.4). 



11.1 Brief Review of Properties of Angular Momentum Operators 

From §9.4.2.3, we know that eigenvalues of angular momentum as determined strictly by commu- 




m--l 



m = 1 



tator algebra are integer or ^ integer quantum numbers: 




j(j + l)^ 2 \jm) 
rah \jm) 
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where j = 0, 1, 2, . . . or j = i, |, |, . . . and m = ±j, . . . , and J is the total angular momentum. 
The commutation relations hold for all types of angular momentum: 



J J 2 


= 






J±, Jz 


= tJ±k 


J+J- 


= 2hJ z 



where j± = J x ± iJ y are the raising and lowering operators. 

J± \ jm) = hy/(j^fm)(j ±m+l) | j, mil) 

11.2 Spin: Evidence in Atomic Spectroscopy 

Let us briefly review the physical context in which electron spin was proposed. Consider a multi- 
electron atom in a magnetic field. If we did not know about spin, we would write the Hamiltonian 
as: 

H = H - H ■ M 

where 7i is the externally applied field and M = -^.L = ^ Sj=i( r i x Pj) 1S the magnetic moment 
obtained by summing over the z magnetic moments of the individual electrons. 

Remark 1: Let's recall that the concept of electron magnetic moment comes from the classical 
definition for the magnetic moment in a current loop. 



magnetic moment /J 




Change e, mass m 



11.2. SPIN: EVIDENCE IN ATOMIC SPECTROSCOPY 



157 



-irr 



1 

c 2~Kr/v 
1 e 1 
c 27rr m 



r x p 
\r x p\ 



x p 



7rr 



r x p 
\f x p\ 



2mc 

e 

2mc 



| pj- Xif| =jj sin(90)=« 
f X p 

L 



Remark 2: Hq is invariant under rotations about the nucleus (coulomb interaction depends only 
on distance between particles), and therefore commutes with all components of angular 
momentum Hq, L 2 , L z form a mutually commuting set of operators. We choose TL to be 
in the z-direction and form a representation in terms of \nlm) where: 

H \nlm) = E% l \nlm) 

L 2 \nlm) = h 2 l(l + 1) \nlm) 

L z \ nlm) = hm \ nlm) 

and Etf energies are 21 + 1 -fold degenerate since they are independent of m. The Hamil- 
tonian is written 

H = H — - — HL Z 

ZfAC 

Since H commutes with Hq,L 2 , and L z , the \nlm) states are also eigenstates of H, with 
eigenvalues: 

-—hUm 
2/ic 

where \xb is the Bohr magneton. Thus, in the presence of a magnetic field, the degeneracy 
is lifted: 

-m = -1 

Ef: 1 = 1 / n 
x -m = 



m = l 



158 



CHAPTER 11. ANGULAR MOMENTUM: SPIN 



Remark 3: The above energy splitting was not always observed in atomic spectroscopic exper- 
iments. If the number of electrons was odd, the number of splitting was even. Also, the 
spacing in energy levels was not always given by figH. Given the even number of levels, 
it seems there must be | integer spin values of angular momentum (2j + 1). Since orbital 
momentum can not give \ integer angular momentum quantum numbers, it was hypothe- 
sized that electrons contain some intrinsic angular momentum that is independent of orbital 
angular momentum. This implies electron spin. 

Remark 4 : The magnetic moment associated with spin is: 

Ms = g S 7, — * 



where 



f2 s = spin magnetic moment 
S = spin angular momentum 



and g s is the factor that produces the correct energy level splitting in spectroscopic exper- 
iments; g s ~ 2. 

Remark 5: Since S is an angular momentum operator, we assume § 2 has eigenvalues h 2 (|) + l) 
for S = i, and that the S z eigenvalues are ±ifo. 

Remark 6: Spectroscopy gave physical evidence of spin. Dirac demonstrated mathematically that 
electron spin was required to make the wave function covariant. 



11.3 Spin Kinematics 

11.3.1 Spin Operators, Spinors, and Eigenvalues 

Consider an infinitesimal rotation of an n-component wave function about the z-axis: 

\ 



IE 



<(n) 



-ih-Sz 



-ih- 



IE 



S, 



(2) 



(1) 



where 



1. Corresponds to a physical rotation which assigns points in space to rotated 
points 
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2. Corresponds to taking components of wave functions and transforming 
them into linear combinations of other components 
Equivalently: 



W) 



I--J 

n 



where J z = L z + S z is the generator of rotations about the z-axis. More generally, 

J = L + S 
is the total angular momentum. 

Remark 1: J and L obey angular momentum commutation relations (§1). Since 
L, J = (L acts on the spacial part of wave functions x,y,z, and S acts on 

components 1, . . . ,n), it follows that S also obeys angular momentum com- 
mutation relations. 



Si , Sj 



Remark 2: In Chapter 9, we found that rotation operators are block diagonal matri- 
ces and each block represents a rotation matrix for a specific value of angular 
momentum. Since we are assuming that the electron has spin eigenvalues 
the relevant block in Z^ 1 / 2 ) , which implies that the electron spin operators are 
given by: 

1 " 

-1 



5;, 



1 

1 



Q - H 

ty- 2 



-i 

1 



2 



Remark 3: Given the spin operators, the electron can be described as a 2-component 
wave function: 



ip-(x,y,z) 



+ V- 



where tp is called a spinor. 



Remark 4 : This can be represented as 



+ = 



1 1 

2' 2 
1 

2'~2 
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In an S z basis, |+) points up the z-axis and | — ) points down. Here, S z in 
matrix form is: 



(+ 



s z (\+) 



■)) 



(+ 
(- 



h , , h 

2 l+) "2 



h (\ o 

2 VO -1 



which is what we said above. 



In the spirit of raising and lowering operators: 

1 



S+ + S. 



1 



Sy ~2i [ 



s. 



where 



s+\+) 


= 


s+H 


= K\+) 


S-\+) 


= h\-) 


S-\-) 


= 



These give S x and S y matrices that are consistant with Remark 2. 
Remark 5: Note that 



S 2 \±) 



S z \±) 



2 



1 + 



1 



2 1 ' 



Remark 6: One important difference between spin and orbital angular momentum is 
that the magnitude of spin is fixed whereas the magnitude of orbital angular 
momentum can change (as in an external field). 



Remark 1: So far, we are only considering the case where 



L,S 



0. In some 



instances, this is spin-orbit coupling, where 



^ 0, but we'll worry about 



that later (perturbation theory). For now, we assume: 

H = H Q + H s 

Since we are taking the spin to evolve separately, the wavefunction separates: 

\m) = \Mt))®\xs{t)) 
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Remark 8: Generalizing the spin components to an arbitrary direction n, rather 
than in the z-direction: 

(n,±\S\n,±) 

are the spin states that point up or down the n axis. So if n points in the 
direction of (6, 6): 



n x 


= sin t 


I COS 


n y 


= sin( 


) sin 


n z 


= cos 


9 



And 



n ■ S = n T 



n z 
n x + in 



n x — in 
— n. 



using Remark 2 



cos v sin oe 



sin »e r —cost 
Solving the eigenvalue problem gives the following eigenvectors: 



cos (f) e 
sin (|) e 

11.3.2 Pauli Spin Matrices 



i<t>/2 
id>/2 



Sill 



COS 



t\ e -«f>/z 

i<f>/2 



(I) 
(I) 



It is conventional to drop the ^ factor in the spin matrices S x ,S y , S z to get the Pauli 
Spin Matrices: 



1 

1 



b = -rr 

2 



-i 

1 



1 
- 1 



Remark 1: The Pauli matrices anti-commute: 

[<rh°j] + = + = ° (* 7^ j) 



Remark 2: a x a y = ia z and other cyclic permutations (by the anti-commuting na- 
ture). 



Remark 3: Tr(<7; = 
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Remark 4- <5f = Moreover, (n • &) = I. 
Remark 5: = 25ijl 

Remark 6: \a x , a y ] = 2ia z and cyclic permutations 

Remark 7: Tr(<7j<5j) = 2£y. If we choose <ro = (J ?)> then 

Tr (a a frp) = 25 a ^; a,/3 = 0,x,y,z 

This implies that the <7 a 's are linearly independent, and thus any matrix M 
can be written 

M = rn a a a where m Q = — Tr (M<7 a ) 

a 



11.3.3 Rotations 

The rotation operator for spin ^ is in a simple form. Consider rotations about the 
z-axis: 



U[R(z)\ = e 



(2fc!) l2 2 / ^ (2fc + l)! \2 2 J 2 

fc=0 v V y V 7 , fr=f K V 7 



Even terms Odd terms 



Since <r? = I: 



= cos [ — ] I — i sin — cr 2 
For rotations about a general axis, this becomes: 



U[R(p)] = co*[-\-isml-]0-& 



11.4 Spin Dynamics 

Consider a stationary electron in a magnetic field B. The Hamiltonian is given by 
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The time evolution for a state is given by 

\m) = H) i#>)> 

where the propagator 



UH) = e ™/* = exp f-JL.S-Bt) 
\ men J 



Remark 1: We see that the propagator closely resembles the rotation operator e t6 S . This 
implies that the propagator rotates the state by an angle 6 = —-^-Bt. 



Consider now the case where B = Bz. Then 

Uit) = exp (^ c S z t) = e 

where uj = —. Let's let 

|V>(0)> = \n+ 



iuj a z t/2 — ( e U 

q e —iujQt/2 



cos (f) e-^/ 2 
sin(|)e^/ 2 



where \n+) indicates a spin up along an arbitrary axis. Then: 

\m) = u(t) i-0(o)> 

'e^ 2 \ /cos (|) e"^ 2 



e'^/y V sin(|)e^/ 2 
cos (|) e -i(*-«ot)/2\ 



sin 



'I) e i(^-wo*)/2 J 
implying that the rotation angle changes at a rate loq. 

11.4.1 Spin \ and the Bloch Sphere 

Let us consider a geometrical representation for the state of a spin-^ system. This 
representation is the work of Felix Bloch and is used extensively in optical physics (a 
spin-^ system is isomorphic to a 2-state system and this makes a nice way to study 
transitions in atoms and molecules). The most general state of a spin-i system is a 
superposition of up and down states: 

M = C 1 \+) + C 2 \-) 

where C\ and Ci are complex numbers. 



Remark 1: C\ and C2 are complex, and thus there exist 4 parameters that specify 
the state. We can reduce this to 2 by noting: 
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• (ip\ip) = 1 =► |Ci| 2 + |C 2 | 2 = 1 

• Overall phase of a state does not matter, we only need the relative 
phase of C\ and Ci- 

Remark 2: Let the 2 parameters be the angles 6, <j) so that we can write 
|V>=cos(£) e-^ 2 |+> + sin(£) e^ 2 |-) 

Let the state \ip) be represented by a point of the surface of the unit sphere. This is 
called the Bloch sphere. 




Figure 11.1: The Bloch Sphere 



Remark 3: The axes on the Bloch sphere have physical significance. Consider: 

= O^> = |+> = |I>3> 

e = -K^\f) = \-) = \\-\) 

Thus the points represent eigenstates of S z with eigenvalues ±|ft. Expressing 
these in terms of Pauli matrices, the north and south pole of the Bloch sphere 
represent: 



(a z ) = l (9 = 0) 
(a z ) = -l (6 = tt) 
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Remark 4 : The expectation values for a x ,a y at these locations are: 



9 = : 



7r : 



<*.) 

<*.) 



+ 



+ 



(_ t -)( 5+ _ 5 _) 



(S+ + S-) 



(-i)(5 + -5_) 



+ =0 



which suggests that the projection of the state onto each of the Cartesian axes 
of the Bloch sphere corresponds to (a x ) , (a y ) and (a z ): 




represents |— ) state with ((<r x ) , (a y ) , (a z )) = (0,0, —1) 



Remark 5: Consider the general case 



cos(|)e-^ 2 |+)+sin(|)e 



i<t>/2 
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<*.) 



cos 



) e ^ 2 (+|+ S in(f)e-^ 2 ( 



e\ -i<t>/2 



S x + S. 



cos e 



-i<j>/2 



+) + sin (J) e 



0\ i0/2 



cos (I) + sin (§) e i<? 
sin(6>) cos(</>) 
sin(0) sin(</>) 



cos 



) e 1 ^ 2 (+| + sin (|) e-^' 2 ( 



i<t>/2 



cos(|)e-^ 2 |+) + sin(|)e 



i0/2 



cos 2 (|) -sin 2 (|) 
cos(fl) 

These correspond to the Cartesian coordinates for a point located at 9, <j> on 
the unit sphere. 



11.4.2 Rotations on the Bloch Sphere 

Consider a rotation about the z-axis. From §11.3.3: 



U[R(z) 



cos 



1 sin 



(!) 



cr 2 



where x 1S the angle we have rotated through. Applying this to the general state: 



U[R(z)} |V> 



cos ($J J — i sin 



[cos (§) 



f)] 



(7- 



cos(|)e-^ 2 |+) + sin(|) e 



„i0/2 



ism{*)\ cos J e 



-i0/2 



|+)+[cos(f)+isin(f)]sin(|)e^ 2 |-) 



cos 



(|) e -^+x)/2| +)+sin( |) e ^+x)/2 



Which shows that a point on the Bloch sphere represented by state 8, cj) gets rotated 
about the z-axis by an angle x- 




U[R(y)} \i>) 



Now, consider a rotation about the y-axis: 
= cos (I) I - i&y sin [cos (|) e"^ /2 |+) + sin (f) e^ /2 |- 



s(f) J-iff y sin (|) 

;$) cos (|) e"^ 2 - sin (f ) sin (|) e^ 2 ] |+) + [cos (§) sin (§) e^ 2 + sin (§) cos (|) , 



cos I 
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This is more complicated, but looking at the case where (f> = 0, (in the (a x ) — (d~ z 
plane), we have: 

U[R(y)] |V# = 0)) = cos (»±2£) |+) + sin |-) 




Remark 1: Since the evolution of a state can be considered a rotation (§11.4), the 
evolution can be described as a point moving with time on the surface of the 
Bloch sphere. 



11.5 Particles with Spin 1 

Electrons gave the first clues to the existence of an intrinsic angular momentum. However, this 
properly applies to a wide assortment of particles. 

Remark 1: Recall that particles with ^ integer spin are called fermions. Particles with integer 
spin are bosons. 

Consider the rotated state of §11.3.1: 



is 



I--(L Z + S Z 



The rotation matrix for this rotation is: 



''cos e — sin s s 
Rz(z) = I sine cose 
1, 
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11.5.1 Rotations of a Vector Field 

(Disclaimer: Sometimes H = 1 in this section) 

Now, if we consider an arbitrary vector field A(r), we can imagine this as having 
a vector associated with every point in space. 




When we rotate the vector field, we are rotating the vectors A as well as the positions 
r. we want to relate this rotation to angular momentum. 

A A', r -» Rr 

So the transformation: 

A'(Rr) = RA(r) 
-or- 

A'(r) = RA(RT l r) (f) 

where the two different forms correspond to active or passive transformations. The 
rotation matrix given for rotations about the z-axis has 

R z \e) = R z (e) 

Thus, the vector field A'(r), according to (f), is given by 

(cose —sine 0\ / A x (R^ 1 r) 
sine cose A y {R~ l r) 
l) yA^R^r) 

=> A' x = A x (R~ 1 r) cose - A y (R^ 1 r) sine 
A' y = A x (R~ l r) sine + A y (R~ 1 r) cose 
A' z = A z {R~ l r) 



where 



' x cos e + y sin e N 
R~ 1 r = J y' J = J y cos e — x sin e 
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Since e is small we can write this as: 



A' = A x \ x cos e + y sine, y cos e — x sin e, z ) — eA v (x + ey,y — ex, z) 

y . " . ' ] 

= A x (x + ey, y - ex, z) - eA y (x, y, z) to 0(e) 



A x (x,y,z) +e 



dA„ 



dx' 



dx' dA x 
+ 



=x , de dy' 



y=y 



i de 



eA y (r) 



d d 

A x {x,y,z)+e{ y —- x -^) A x (r) - eA y (r) 



it 



te 



I--L z )A x {r)-eA y {r) 



Similarly, 



A' y (r)=eA x (r)+(i-^L z )A y (r 



I-jL z j A x (r) 
- |4) A(r) - j (-iA y (r)x + %A x {r)y) 



Usual Rotation 



Spin part 



Thus we have 



A'(r) 



is 



I--[L Z + S l 



A(r) 



where we define 



S z A(r) = -^A y (r)x + l -A x {r)y 



since 

'A x {rf 
.Mr), 

Therefore, we can represent S z as the matrix 



-iA y ( r y 

iA x (r) 




S z 



'0 -i 0> 
i | 
,0 Oj 
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Similarly, rotations about the x and y axis give 

/0 \ / i\ 

s x = n o o -i , s y = h o oo 

\0 i / \-i 0/ 



Remark 1: S x ,S y ,S z are related by row-column interchanges. 

Remark 2: if S x ,S y ,S z are components of angular momentum, they must obey 
commutation relations. This can be checked directly. For example: 

/0 0\/0 OA /0 OA/0 \ 

&4 - S„Se = & -i - h 2 -i ) 
\0 i / \-i 0/ \-i 0/ \0 i / 

/ 0\ /0 -1 0\ 

= ft 2 -l o o - h 2 o o o 

\0 00/ \0 0/ 

/0 -i 0\ 
= /i 2 i i 
\0 0/ 



Remark 3: The magnitude of angular momentum is given by 



'0 x 



,0 i 

A N 
2 10 
V0 1, 












A) 


— i 




(• 








H 














oy 




VO 





0/ 



= s (s + i)i 

s( s + 1) = 2 =^ s = 1 => Spin 1 



Remark 4 : Working in the standard representation for a s = 1 system, |l,m), with 
S z |1, m) = mh |1, m) ; m = 0, ±1 
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We find that, in this representation, 

S z = (1,0| S 2 (|l,l) |1,0) |1,-1)) 




To find the other components, we'll use the relations from §11.1: 
5+ 1 1,0) = ftV(l-0)(l + + l)|l,l> = hV2\l,l) 



5+|i,-i) = nV(i + i)(i-i + i)|i,o) = W2|l,0) 
s~ |i,i) = hy/(i + i)(i -i + i) 1 1, o> = W2|i,o) 

5+|l,0) = ^v / ( 1 -°)( 1 + ° + 1 )l 1 ) - 1 ) = ^\/2 |1, — 1) 



So that 

5* x i 





S P y^ [ (1,0| |4(|1,1) |1,0) |l,-l))=-L[l 






1 (A 






1 


1 





i oy 







-1 















1 






Remark 5: Notice that the matrices in Remark 4 are different than those given at 
the beginning of this section. This is because the first set are identified with 
vector cartesian components, and the matrices in Remark 4 form a basis that 
is not identified as such. To understand, recall that we defined the "spin" part 
of the rotation by 

S z A(r) = -iA y (r)x + iA x (r)y 

Let 

u± = — (x± iy) 
with «o = z. Then any vector can be expressed in terms of these as 



A — A x x + Ayij + A z z 

= ~^ A x (U+ + U-) 

= A + u+ + A_u— + vlo^o 



1 1 

—=A X (u+ + u-) H -=A y (u+ - u-) + u A z 

y2 iy 2 
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with 



A± = -j= (A x t iA y ) , A = A Z 



So 



S z A(r) = -1 (u+ + u_) -J= (A+(r) - A_(r)) + («+ - «_) (A + (r) + A_(r)) 
= u+A+(r) — M_^4_(r) 
And thus 



S z 



Matrix from Remark 4 

Similar arguments work for S x ,S y . Physically, this change of basis can be 
related to the polarization of a photon, which is a spin 1 particle. These spin 
states of a photon are associated with circular polarization in the x-y plane 
where m = dbl, and linear polarization in the z direction (m = 0). Hence, it 
is not associated with the 3 linear polarization directions. Therefore, u+,U— 
and A + , A_ describe circular polarization of an EM field. 






Chapter 12 

Addition of Angular Momentum 



Consider the situation in which we have many particles that make up a system, and this system 
has rotational invariance, and therefore angular momentum commutes with the Hamiltonian. The 
total angular momentum of the system is the sum of the spin and orbital angular momentum for 
all the particles. Given that we know the eigenstates for the individual pieces, we want to find the 
eigenstates of the total angular momentum. 

12.1 Addition of 2 Angular Momenta 

Consider 2 angular momenta, j\ , j 2 , and let 

J = 31+ h 

be the total angular momentum. The eigenstates for the individual pieces we already know: 

3i b'l, mi) = hji(ji + 1) |j'i,mi) 
jiz |ji,™i> = mi^|ii,mi) 

We have similar results for j 2 . Let 

\ji,j 2 , mi,m 2 ) = |ii, mi) ® \j2,m 2 ) = \ji,mi) b'2,^2) 
be a representation based on the common eigenstates of jf , 3\ z ^3iz- 
Remark 1 : Since J z = j± z + j 2z , 

J z \ji,j2,m 1 ,m 2 ) = H(mi +m 2 ) j 2 , m 1 , m 2 ) 
So the state \j±, j 2 , mi, m 2 ) is also an eigenstate of J z . 

Remark 2: \ji, j 2 ,mi,m 2 ) is not an eigenstate of J 2 , so we need to find the eigenstates of J 2 . 
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Remark 3: Since j\ commutes with every component of j\ and j 2 (and similarly for we have 
that 





3iJ 



31 J 



So we take j 2 , j 2 i J 2 , Jz as our complete set of commuting observables, and instead of 
|ii> J2> mii ^2); w e use J2, J, M) (so j\ z and j 2z are replaced by J and M). In this 
representation, we have: 

fx J2, J, M) = ^ 2 ii(ii + 1) |ji, j 2 , J, M) 
32 b'1^'2, J, M) = ^ 2 i 2 (i 2 + 1) |ii, j 2 , J, M) 
J 2 \h ,j 2 ,J,M) = hJ(J+l) |ji , j 2 , J, M) 
£ lii , J2 , J, M) = HM \jt , 32 , J, M) 



With this representation, we need to find 

1. The allowed values of J and M for a given ji,j 2 , and 

2. How to express |ji,j 2 , J, M) as a linear combination of | j±, j 2 , mi, m 2 ) 

12.1.1 Finding Allowed Values of J and M 

For individual particles, we have allowed m values given by 

mi = -ji, -ji + 1, . . . ,+ji 
m 2 = -j 2 , -j 2 + 1, . . . ,+2i 

and 



hz \h,j2,mi,m 2 ) 
hz \hJ2,m 1 ,m 2 ) 



hmi \ji,j2,m 1 ,m 2 ) 
frrn 2 |ii,j2,mi,m 2 ) 



and 



■/* |ii,j2,mi,m 2 ) = ft(roi + m 2 ) | Ji, j 2 , mi, m 2 ) 
So we expect the M eigenvalue to be given by 

M = -(ji + j 2 ), -(ji + j 2 ) + 1, . . . , +(ji + j 2 ) 



Remark 1: In general, there will be more than 1 state with the same M eigenvalue, 
but these states will have different J eigenvalues. 
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Remark 2: In order for us to get the allowed M values as \M\ < j\ + j 2 , we need a 
J = ji + ]2 eigenvalue. This would give: 

J ■ h+h^M = -(ji + j 2 ), + j 2 ) 

ji+32-l->M= -Ui+h) + !,•••, (ii +h) - 1 
ji + h ~ 2 -> M = -(ji + j 2 ) + 2, . . . , (ja + j 2 ) - 2 



But where does the sequence terminate? We can determine this with vectors. 
If ji and 32 can be represented as vectors (and they can), then the maximum 
angular momentum should be obtained when they are parallel: 




Max Magnitude of J = ji + j 2 



The minimum occurs when they are anti-parallel: 




Min Magnitude of J = \ji — J2I 



Remark 3: We can now find the total number of \jx,j2,J,M) states. Without loss 
of generality, let j\ > j 2 - Then 

# states = [2(j! + j 2 ) + 1] + [2{ji + j 2 - 1) + 1] + . . . + [2{h - j 2 + 1) + l] + [2(ii -32) + A 

V v ' 

Allowed M if 
J=ji+h 

Many of these terms cancel, leaving only (2ji + 1) in each term. Since there 
are (2j 2 + 1) terms, the number of states is 

# states = (2ji + l)(2j 2 + 1) 

which is the same number of states as we had in the \j±, j 2 , mi, m 2 ) basis. So 
we have the correct number of J and M eigenvalues. 



Remark 4' To summarize: for a fixed ji and j 2 , the (2ji + l)(2j 2 + l)-dimensional 
space spanned by \ji,j 2 ,mi,m 2 ) has possible values of J given by ji +j 2 ,ji + 
j 2 — l,...,\ji — j 2 \, and for each J, there are (2 J + 1) states \J, M) with 
different M values. 
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.1.2 Example of Addition of Two S=| States 

Consider 2 particles with spin ^. The individual eigenstates are: 

51 = \hh)i and li'-i>i and 

5 2 = |f, \) 2 and |3,-5) 2 =► |+) 2 and |-) 2 

We define 

S = S 1 + S 2 

and we denote the eigenstates of 5 2 and S z by |5M). 

Remark 1: For the 4 |si, S21 mi, 7712) states: | — ) 2 ; |+)i \— ! I - )i l+)2 i l — )i I - )2> 
there are 4 ISM) states: 

s = i=Hi,-i>,|i,o),|i,i> 

s = ^ |0,0) 

The state with S 1 = 1 is called a triplet and the state with S = is called the 
singlet. 

We now want to express the \SM) states in terms of the \s\, S2,mi,m,2) states. Note 
that there is only 1 state in each set with M = mi + m 2 = 1 and one state with 
M = mi + mi = —1. So we can write: 



I1.D = l+)il+) 

->l|->2 



l\ IT-/2 

11,-D = 



We can find the other 5=1 states by using the lowering operator: 

|1,0) = -^5_|1,1) 

= -L(5i_+5 2 _)|+) 1 |+) 2 

= ^(|-)il+) 2 + l+)i|-) 2 ) 

For the state corresponding to 5 = 0, we must have a linear combination of |— ) 2 , 
which is orthogonal to the other states: 



= 4= (i 
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Remark 2: In the \s\, S2, mi, m.2} basis, we have a tensor product of i ® i, which is 
equivalent to 1 © in the \SM) basis. In other words, 

\®\^ (2si + l)(2s 2 + 1) = (| + 1)(1 + 1) = 4 
1 

1©0^^(25 + 1) = 1 + 3 = 4 

s=o 

which must be equal. 



12.2 Two Nucleon System 

Let's consider a system of 2 nucleons each with spin „-> but with the addition of orbital angular 
momentum. Only the angular momentum of the relative motion is relevant, so the total angular 
momentum is given by: 

J = L + S 1 + S 2 = L + S 

where S = Si + S2 is the total spin of the system and L is the total angular momentum of relative 
motion. We then have 

\L, S, Ml, Ms) = \L, M L ) ® \S, M s ) = \L, M L ) \S, M s ) 

with 

L = 0,1,2,... M L = -L,...,L 
5 = 1 M s = -1,0,1 

5 = M s = 



L + 1,L,L-1 L>1 
1 L = 



Giving the allowed J values of 

For 5 = 1 J = I 
5=0^ J=L 

The states of total angular momentum are given by \L, 5, J, M). 

Remark 1: The quantum numbers L, 5, J are related to spectroscopic language: 
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Total Spin 

2s+l 



L D, _ 

jf J I p=1 

i Gives L e-values according, to < D=2 

Total Ang. Mom. * I F — 3 

Total L+S \G = 4 

(a) General Form (b) Specific Instance 
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12.3 Clebsch-Gordan Coefficients 

In §12.1, we noted that for a given ji,j 2 we needed to find the allowed values of J,M (done in 
§12.1.1) and an expression for |ji,j2, J, M) in terms of |ji, j 2 , mi, 7722) (done as an example of 2 
spin ^ particles). We now wish to come up with a general scheme for expressing the eigenstates of 
ji,j 2 ,J,M. Formally, we can write: 

\ji , j 2 ,J,M) = y~] (ji , 32 ,mi,m 2 \ji , j 2 , J, M) \ ji , j 2 ,m 1 ,m 2 ) 

m~m V v ' 

11 2 Clebsch-Gordan Coefficients 

Remark 1: The Clebsch-Gordan coefficients depend only on the angular momentum quantum 
numbers, not on the specific details of a physical system. Thus they can be calculated 
from the addition of angular momentum in general (without repeating for each physical 
application). 

Remark 2: In practive, extensive tables of these coefficients are available. 

Remark 3: To simplify notation, let 

Oi , 32 , mi , m 2 \ji , 32 , J, M) = (ji , j 2 ,mi,m 2 \J, M) 



Remark 4- We already have some of the Clebsch-Gordan coefficients from previous examples: 



From §1.2 ^ |0, 0) = |-) 2 - -L 



Similarly, 



+ 



T2 mi 



V2 



mi 



2,m 2 = 
-\,m2 



otherwise 



(iimim 2 |ll) 



1 mi = 7TT-2 : 

otherwise 



■11 

,22 



mirri2|l, — l) 



1 mi = 772,2 : 

otherwise 



(ii772im 2 |l0) 



772i 
772i 



1 

V2 
1 

V2 

otherwise 



2-,m 2 = 
-\,m 2 
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Remark 5: Note that there can be an ambiguity in the overall phase (minus sign) of the states 
| JM). This is fixed by convention. 

Remark 6: General properties of Clebsch-Gordan coefficients: 

1- {jih m i m 2\J M) = unless M = mi + m 2 and j\ + j 2 < J < \ji — j 2 \- 

2. They are real. 

3. The convention for sign is that 

0'iiaii, J-ji\JM) > 

where mi = j\ and m 2 = J — j±. 

4. The equation 

(jij 2 m im2 \JM) = (-l)-'' 1+ ^- J 1 j 2 ,-7n 1 ,m 2 |J,-M) 

gives the coefficients of — M in terms of M. 

5. Starting at the bottom (top) state and applying J + ( J_) gives recurrence 
relations between successive coefficients: 

\J J(J + 1) - M{M + 1) (jij 2 m im2 \J, M + l) = 

Vhih + 1) - mi (mi - 1) (jij 2 , mi - l,m 2 \JM) 
+ y/h(32 + 1) - m 2 (m 2 - 1) (j'ij2mi,m 2 - 1| JM) 

and 

y/J(J + 1) - M(M - 1) (jihm im2 \J, M - 1) = 

Vhtii + 1) _ rn\{m\ + 1) (jii 2 ,mi + l,m 2 \JM) 
+ Vi2(j2 + 1) - m 2 (m 2 + 1) (jij2'mi,m 2 + 1\JM) 



Remark 7: The transformation from the |jij 2 mim 2 ) basis to the \j\j 2 JM) basis is an example 
of a unitary transformation between 2 orthonormal basis. Thus, for jij 2 : 



such that 



and 



matrix elem. \ co i. i n d ex 



JM 

row index' 



uu ] = U ] U = I 
^2 {hhmim 2 \J M) (JM|jii 2 mim 2 ) = <S mi , m ><5, 



m,2,m 9 



JM 



^ {J M\iij 2 mim 2 ) (jij 2 mim 2 \ J' M') = 8jj>8. 



MM' 



m\m2 
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12.4 Irreducable Tensor Operators 

12.4.1 Irreducable Invariant Subspaces 

Consider a state \X,j,m) such that 

J 2 \\j,j) =hj(j + l)\X,j,j) 
Jz I A, j, j) = hj \X,j,j) 

and the (2j + 1) ortho-normal states obtained by §11.1: 

• v / (j' + m)! M-m u • .v 
|A ' J ' m) = V(2^-m)! J - ' A ' J ' J) 

which span the subspace of Hilbert space having eigenvalues of J 2 and degeneracy 
index A. Let £^ denote this subspace. Rotations involve only operators J x ,Jy,J z 
(or, alternatively, J + , J_,J Z ), so that a state that is initially in £^ is rotated into 
another state in We say that is invariant under rotations. 

Remark 1: Consider a state in £^ J '. All possible rotations of will span the 
entire Because of this property, we say that £^ is irreducible under 

rotations. 



Remark 2: These rotations can be represented in standard representation by a ma- 
trix: 



&( Q .A7) = (A, J,M 7) 



A, J, M' 



(2J+l)x(2J+l) 
matrix 



The K^-j^/ are called rotation matrices and represent rotations in Hilbert 
space. The rotated basis vectors R\X,J,M) can be expressed in terms of 
the non-rotated basis using these rotation matrices: 



R | A, J, M) = ^2\X,J, M') /A, J, M' 



AT 



R 



X,J,M 



AT 
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12.4.2 Irreducable Tensors 

Scalars, vectors, and tensors are distinguished by their transformation properties. For 
instance, consider transformations under rotations: 

Scalar: S -» S' = U ] SU = sU^U = S 

So scalars are invariant under rotations 
Vector: v ; — > v\ = U^viU = RijV^ 

and, equivalently, 



Vi ->• v'i 



Uv^ = R 



i.i 1 



since U[R] -► [/[iT 1 ] = U ] [R\ 

Definition 1: A rank k tensor operator is a collection of (2k+l) operators denoted 
Tg k \q = —k, —k + 1, . . . , k that transform according to: 

f(k) _^ ftjWff = f q/ R q ' q (f) 



Remark 1: These tensor operators are irreducible under rotations. 

Remark 2: Equation (f) represents a transformation of a group ("basis") of opera- 
tors. The equivalent and more familiar case is a transformation of states: 

R\k,q) = J2\k,q')Rf g 

{kq'\R\kq} 



12.4.3 Wigner-Eckart Theorem 

This theorem is widely used in atomi and nuclear physics and is used to calculate 
matrix elements for evaluating transition rates between states. It says that the matrix 
elements in a standard representation of the q th component in a k th order irreducible 
tensor operator is given by: 



XJM 



T (k) 



X'j'M' 



1 



V2J + \ 

unimportant 
conversion 
factor 



(AJ| 



T (k) 



| A' J') (j'kMJjM) 



reduced matrix clement 



CG- coefficient 



Remark 1: The Wigner-Eckart Theorem does not tell us what the reduced matrix 
element is - this depends on the particular tensor operator under consideration. 
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The value lies in the way the q, M, M' quantum numbers are extracted in 
the Clebsch-Gordan coefficients. For example, in atomic physics, the M, M' 
quantum numbers for the initial and final states in a radiative transition are 
connected with the polarization of the emitted photon. In this example, the 



Wigner-Eckart theorem tells us the ratios of radiative decay rates for different 
polarizations. These ratios are given by ratios of the squares of the Clebsch- 
Gordan coefficients. 

Remark 2: The Clebsch-Gordan coefficients {J'kM'q\JM) give the angular momen- 
tum selection rules: 



1 = 1: M=-l M = M = l 

•* >v Linear Polarization ^ 




Linear Polarization 



=0: 



q = M -M' 
J-j'\<k<J + J' 
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Chapter 13 

Perturbation Theory 



Most QM problem can not be solved exactly, thus approximations are an important part of the 
actual applications. Perturbation theory is one of the main approximation techniques used. The 
idea is to approximate eigenvalues and eigenstates of a Hamiltonian which differs slightly from a 
Hamiltonian whose eigenvalues and eigenstates are known exactly. 



13.1 Perturbation Theory for a Non-Degenerate, Discrete, Energy 
Level 



Consider the Hamiltonian: 

H = H + XV 

where Hq is a Hamiltonian whose eigenvalues and eigenstates E n °\ |n^) are known. So 



Ho 



n 



(o) 



(0) 



For now we consider only the case where the states are discrete and non-degenerate. XV is a 
perturbation with A the size of the perturbation and V has the same "size" as Hq. The eigenvalues 
and eigenstates of the full Hamiltonian are given by: 



H n In) = E n In) 



where we write 



So that 



E n = E^ + XEP+X 2 E^ + ... 



n 



lim E n 



\n) + X 



and 



(2) 



+ 



(t) 
(t) 



lim In) 



(0) 



Substituting (|) into (f) gives: 



(h 



XV 



n (0)^ +A | n (D^ A 2 



n 



(2) 



+ 



4°) + XEW + A 2 M 2) + • • • 



n 



(0) 



+ A 



n«)A 5 
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Now, equating power of A on both sides: 



Hn 



n (fc-l)^ = E (0) \(k)^ + E {1) l(fe-l)^ + ... +E (k) L(0)^ 



or 



' (h - E^) = - (V- - |n(°)> , fc = 1 



where terms with superscripts (k) refer to the fcth order term in the perturbative expansion. 

13.1.1 First Order Approximation 

Consider k = 1. Then we can apply (n^^l to 



H n -E^ 



n 



(i) 



V — E, 



(i) 



(tt) 



n 



(0) 



H n -EM 



n 



(i) 



,.(0) 



n 



(o) 



= - (n 



,(0) 



r 



n<P>)+E<P 



(0) 



— 'M 



71 



(0) 



V 



n 



(o) 



Thus the first order correction to the energy is just the expectation of the perturbation 
taken with respect to the unperturbed state. To get the 1st order correction to the 
eigenstate, we must expand [n^ 1 )) in the basis of the unperturbed state: 



77 



(1) 



??? 



(0) 



77 



(1) 



77? 



(0) 



where m is the unperturbed state. To find the coefficients, we will apply ^m^l to 
(tt): 



77? 



(0) 



H n -Ei°A n W 



77? 



(0) 



v - eP 



,(») 



77 



(0) 



777 



(0) 



7? 



(1) 



7)7 



(0) 



V 



ir 



(o) 



p(0) p (0) 



where the term goes to zero since we are considering when m ^ n and thus 

( m (0)| n (°)> = 0. 
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Remark 1: At this point, it becomes important that we are considering a non- 
degenerate state since this trick will not work if Em = En for m 7^ n. 



Remark 2: In the expansion of Ira^) in terms of |m^ ^, there will be a coefficient 
corresponding to n = m which appears to yield a diverging coefficient. To fix 
this, we set 



n 



(i) 



n 



(0) 







We can justify this by considering 

| n ) = | n (0)^ +A | n (D^ +A 2| n (2)^ + . 

If |n^) has an |n^ ^) term in its expansion, then 



n W\ = / n (o) 



n 



(k) 



n(°) ) + 



n 



where the last term is an expansion of |n^ fe ^) in terms of |m^ ^ with m 7^ n. 
Thus we have 



n 



n 



a 



(o) 



00 \ / 

l + ^A fc (n(°)| n W) +A| n «) +A 

k=i ) v 



all n=m terms here 

(o)^ + A | n (i)y + A 2| n ( 2) y 



expansions with n ^ m 



+ 



where 



n 



(0) 



fe=i 



is a scaling factor that can be absorbed into the length of |n). Dividing by a 
and letting ^ |n) — ► |n) and |n^ fc ^) = — \n^ k ^ gives: 



|n) = |n(°))+A|n( 1 )) // + A 2 |n( 2 )) ,/ + 



which is the expansion of |n) with states having no projection onto the 

state |n(°)). This is equivalent to writing 



|n) = |n(°))+A|n( 1 ))+A 2 |n( 2 )) 



+ 



and setting 



n 



(0) 



n ( fe ))=0, Vfc = 1,2,3,... 
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Remark 3: The procedure used to obtain the first order correction in eigenvalues 
and eigenstates can be applied successively to get higher order corrections in 
terms of previous determined corrections of lower orders. Taking terms \ k 
with k > 2 implies: 



(Ho ~ 4 0) 



(/v) 



V - J5« 



n 



(fc-i) 



n 



(fc-i) 



i=2 



on which we apply (riA°)| to get 



n 



(o) 



(o) 



7i 



(fc-i) 



i=2 



n 



(0) 



n^ 1 ) ) + 



k-2 



E^+J^Et 



V 



/ n (o) 



n 



(i; 



t=i 



V 



(fc-i) 



Now, applying (m^l to terms 0(X k ) implies 



(0) 



Ho - En 



n 



(fc) 



m 



(0) 



V - EW 
( 



i=2 



n 



??? 



(0) 



(fc) 



- £ -'m - C/ n 



fe-1 



(0) 



n (fc-i)U ^ #oo / m (o: 

i=2 



(fc-i) 



V 



to fc- 1 since (m<°) \nW)=0 J 



which are coefficients for the expansion 



n 



(fc) 



■m 



(0) 



(fc) 



■/7J 



(0) 



13.1.2 Example: Harmonic Oscillator in an Electric Field 

Consider a particle of charge q and mass m in a harmonic oscillator potential and 
corresponding to an unperturbed Hamiltonian: 

{r f , 1 2-2 

iin = — mw x 

u 2m 2 
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Applying an electric field with magnitude / along the +x direction corresponds to an 
electrostatic potential = —fx and potential energy 



V = -qfx 



The full Hamiltonian then is 



H = H + XV 



where A is the strength of the perturbation. The first order correction to the energy 
is given by 



(o) 



V 



n 



(0) 



- q f 



n 



(o) 



-qf 



2mu> 



1/2 



n 



(0) 



a) + a 



n 



(o) 



Remark 1 : Physically, the average interaction of the particle in the electric field is 
since the particle's unperturbed wavefunction is a H.O. 



4 1} = -<?/ / (V£ 0) )*z (Vi 0) ) dx = -qf 



OO 9 

(0) 



x dx = 

odd 



The first order correction to the eigenstate is: 



??? 



(0) 



n 



(i) 



(mW\(-qfx)\nW) 



qf 



ft \ x / 2 <m(°)|a + at|n(°); 



p(o) _ p (o) 



P (0) P (0) 



(\/n(5 m ,n-i + V™ + 15 min+ l) 



n 



(i) 



^ Z F (0) 1 F (0) (v^m.n-1 + V^TT^^+i) 
m c/n 



7?? 



(0) 



9/ 



2muj 



(n-l)(°)) V^+|("- + l) (0) > 



+ 



P (0) p(0) 



-hoj 



hw 



9/ vS^(^ TT l ( " +1,(01 )-^|("- 1 » < ° ) )) 
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So, to first order, our eigenstate is: 
\n) = \n^) + x\n^) 

= n<°>) + (V^+l\(n + 1)( )) - ^ |(n - if')) 

Since the first order energy correction was 0, let's calculate the 2nd order correction: 



n 



(0) 



V 



(i) 



qf(n 



,(°) 



11 



(1) 



9/ 



2mw 



(o) 



a + a 



(qf? 

2mu 2 

Ml! 

2mw 2 

Ml! 

2mu 2 



2mw' 



(o) 



a + a 



(i) 



,n + 



l)(o) 



(o) 



(n + 1) 



(o) 



vWn- 1 (n - 2) (0) \ + Vn + Wn + 1 (n + 2) 



\(0) 



n 



(o) 



(n + 1 — n) 



Remark 2: This is a problem that can be solved exactly, which makes it a nice test 
for perturbation theory. Consider: 



H = 1 — mu> 2 x 2 — Xqfx 

2m 2 



p2 

H — muj x 



1 



2m 2 



Mf 



raw 



1 q 2 f\ 2 



2 muo 2 



shifted H.O. 



const, potential that 
we don't care about 



resulting in energy eigenvalues given by 

E n = 4°) - £££ = 4 0) + +X 2 E^ 

energy shift 

and eigenstates of the H.O. with oscillator displaced by: 

qf 



muj z 



qf [rruj 



raw 



muj 2 V 27i \j2mhm>J 2 



\n) + D(Xx) 



n 



(0) 
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where D(Xx) is the displacement operator of §6.3.3. 

,Ax(at+a) ^i + Xx(a j + a 



n ) 



11 



(0) 



D(\x) = e 
Mf 



Vn+l 



(n + 1 




1 



,(0) 



\j2mfrufi 

Here we note that computing higher order terms in the perturbation expansion 
is equivalent to taking higher order terms in the D{\x) expansion. Also, note 
that the energy eigenvalue is exact at 2nd order-computing higher order terms 
in E n will give contribution. 



13.1.3 Example: Stark Effect 

See Homework! 

The picture is: Here we have that L \\ u and the rotation occurs in the plane defined 




by the angles 9 and (f) relative to the Cartesian axis. 

: 2 , f 2 , f 2 



E = -Ioj 2 
2 



2 I 



L x -\- Ly -f- L*. 

T 



with 



Hint: 



cos BY, 



1 



mD 2 



(l + l + m)(l + l-m] 
(2Z + l)(2Z + 3) 



2 

nl/2 



ym i 

Y l+1 + 



(/ + m)(l — m) 
(2Z + 1)(2Z-1) 



1/2 



Y l-1 



We will see that 2nd order perturbation removes the degeneracy in the unperturbed 
eigenvalue. 

Note: a degeneracy remains with respect to the sign of m which physically corresponds 
to the configuration of a rigid rotor by invariance with respect to 9 — > tt — 6. 
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13.2 Degenerate Perturbation Theory 

Now, consider the case where the unperturbed eigenstates are degenerate. Start with the simplest 
case of 2-fold degeneracy (higher order degeneracy is straight forward). 



Ho 



n(°),l 
nC°> 2 




(1) 



E ni K 1) 
E n2 \n,2) 



(2) 



where the l's and 2's and the degeneracy index. So eigenstates |n^°',l) and |n(°\2) (which are 

assumed to be orthonormal) , give the same eigenvalue E$ ■ We expect the full solution to look 
like: 

H\n,l) 
H\n,2) 

to give E n -> E^ as A -> 0. Also, for now, assume that one eigenstate approaches |n^°\l) and 
the other approaches |n(°)). To simplify notation, lets let: 

n( k \j) = \k) j 
with j = 1,2 and k = the order of perturbation. Then, 

E ni =E^ + XE^ + ... + X k E^ + - ^ 
E n2 = 4°) + \E$ + ■■■ + \ k Ei k J + • 

and 



(3) 



"2 



112 



(4) 



|n,l) = |0) 1 + A|l) 1 + --- + A fc |fc) 1 + ---| 

|n,2) = |0) 2 + A|l) 2 + --- + A fc |A:) 2 + ---J 

Now, substitute (3) and (4) into (2) to get the expansion in terms of powers of A. As outlined in 
§13.1, the terms to 1st order in A give: 



^o-^jliK 

'h -E(°A\D 
Multiplying on the left by (0^ gives: 



V-E^IO), 
V-E^)\0) 2 








(5) 

(6) 
(7) 
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While multiplying on the left by 2 (0| gives: 

= - r 
= - 



(8) 
(9) 



Remark 1: This analysis is sufficient if equations (7) and (8) hold. In reality, however, they may 
not. In fact, if they do hold, this situation is identical to the non-degenerate results (as in the 
rigid rotor problem) and we've gained nothing. If (7) and (8) do not hold, then we have to 
do something new. In this case, what we've done must be incorrect since (7) and (8) would 
be contradictory. What went wrong above was the assumption that \E n ,l) and \E n ,2) 
approached specific |0) 1 and |0) 2 respectively as A — > 0. In fact, the full eigenstates can 
hypothetically approach any mutually orthogonal linear combination of those states (since 
any linear combination of these states is also an eigenstate of Hq). We must determine the 
linear combinations for the degenerate perturbed states. 

Let us proceed without making any assumptions on the actual form of the unperturbed degenerate 
eigenstates, (let those be mutually orthogonal, normalized, linear combinations of the degenerate 
states), such that the expansions given by (4) still hold. 



\n,l) = |0) 1 + A|l) 1 + --- + A fc |fc) 1 + 
|n,2) = |0) 2 + A|l) 2 + --- + A fc |fc) 2 + 



(4) 



where |0) 1 and |0) 2 are now unknown states. Then the 1st order in the perturbation expansion still 



gives 



(5) 



(h -E^)\1) 2 = -(v-E^)\0) 2 
where we multiply (5) by any (n^\ for which (1) holds to get on the left hand side. Formally 



we can express this as: 



p(y-4 1 1 ))|o) 1 = o 
p(y-4 1 2 ))|o) 2 = o 



(6') 
(7') 



where P = , l) (n^, l| + |n( ), 2> (n^, 2\ is the projection onto the subspace of degenerate 
states. |n(°),l) and |n(°),2) for m an orthonormal basis over this subspace. In practice, this 
calculation reduces to the diagonalization of a 2 x 2 matrix. To see this, let's let 



|0>!=< 

where a, (3 are unknown. Then (6') gives 



n 



(0) 



(0) 



n 



(0) 



l)(n(°),l 



+ 



(0) 



V - 



a 



n 



(0) 



1 +/? 



n 



(0) 
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n 



(0)1 



n 



(0)1 



1' 



+ 



(0) o 



77 



(0) o 



V 



E^)a + {n^,l\v\n^,2)(3 



(0)i 



+ 



a 



V 



s£>W(n< \l 



?i (0) ,2)/3 = 
a = 




((nW, 2|y|n(°), 2) - E$) (3 + (n®, 2\vU°\ l) 





1 

1 



(9') 



which is the standard 2x2 matrix eigenvalue equation. Solving the eigenvalues, eigenvectors, and 
normalizing gives first order energy corrections, and the right linear combinations of the unper- 
turbed degenerate ground states. 

Remark 2: Solving the characteristic equation given by (9') can result in 1st order energy correc- 
tions (En^ 2 ) that are the same or distinct, depending on the roots of the quadratic equation. 
So, to first order, the degeneracy may or may not be lifted. 



13.2.1 Example: Coupled Harmonic Oscillator 

Consider 2 identical H.O.'s with frequency lo that are coupled by an interaction en- 
ergy proportional to the product of oscillator displacements from equilibrium. The 
Hamiltonian is: 

H = H + V 

where 

H = fwj (a) a + tfb) 
V = hg(a + a f ) (b + 6 f ) 

and where a, a) and b, ¥ are annihilation and creation operators, and we drop 2 factors 
r^hui in Hq. 

The unperturbed energy is 

En = E n + E m = Nfouo = (n + m)huj 

And the degeneracy arises because there is more than 1 way to break up the energy. 
All we require is that N = n + m and the degeneracy increases with increasing energy. 
n and m are numbers of quanta in independent, uncoupled oscillators. 

States Degeneracy 

1 0,0) none 

|1,0),|0, 1) 2-fold 

|0,2),|1,1),|2,0) 3-fold 
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Consider 1st order perturbation theory for the 1st excited state (2-fold degeneracy). 
In the limit g — > 0, the exact eigenstate goes to 



E 



(0) 

iVl,2 



a|0,l) +/?|1,0) 



The matrix elements required in (9') of the previous section are 



0,1 
1,0 
0,1 
1,0 



hg [a + a+ J (b + 
kg [a + a f ) (b + 6 f 
hg (a + a f ) (b + b ] 
kg (a + aA (b + 6 f 



0,1 
1,0 
1,0 
0,1 






hg 
hg 



And (9') becomes 



-E 



(i) 

Ni 



hg 



-E 



(i) 



V h 9 ~^N 2 . 
Solving the characteristic equation gives 



(i) 

(1) 



EV T ! = hg 



E% = -hg^ 

And the ground state, unperturbed eigenstates follow: 



E 



E 



4= (I 
1 



10,1) 



|i,o» 



Now we need 1st order corrections to the eigenstates 

E$y 



N ^ 



-E 



jr(0) -I 

N ' X 



E 



V 



N > X 



(0) 
M 



(0) 

V! 



(10) 



(11) 



(12) 



Remark 1: This is the same procedure used in non-degenerate perturbation theory, 
but the summation is more restrictive since E^ and E^ have no projection 
on the degenerate subspace. 



Remark 2: An expression similar to (12) gives 



EP,2). 
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To get the 1st order correction to the eigenstates, calculate the matrix elements 



E 



V 



EP,(1 or 2 



and specifically for the H.O.: 



n,m) hg(a + a ] ) (b + tf) (|0, 1) ± 1 1 , 0> ) = (n,m\ -^L y/2 |1, 2) ± \/2|2, 1 



all unperturbed 
states other than 
11,0} & |0,1> 



V2 



<1,2|-^(|0,1>±|1,0» = % 
<2,l|-^=(|0,l)±|l,0)) = -% 

Remark 3: The states |1,2) and |2, 1) both have energy E^ = 3Hlo. 



Plugging (13) into (12) gives: 
,(l)\ % 



E 



iVi 



3hui — huj 



1,2 



3?ia> — ^iw 



2cj 



(|1,2) + |2,1)) 



-^(|1,2)-|2,1)) 



(13) 



(14) 
(15) 



Remark 4 : The interacting Hamiltonian in this example is often approximated by 

Vrwa = (a'b + atf 



This gives the 1st order corrections that are exact (no higher order approxima- 
tions from a) b, ab>). This interaction is referred to as the rotating wave approximation. 



Remark 5: In this case, the 1st order corrections lifted the degeneracy, and higher 
order corrections may proceed with non-degenerate perturbation theory. If 
degeneracy was not lifted at first order, we must diagonalize another matrix 
(use non-degenerate perturbation theory) again. 

The second order corrections in energy for coupled H.O. are given by 
p(2) _ / „(o) , 

= -L (|0, 1) + (1, 0|) % (a + at) (S + St) (Z£J (|l, 2) + |2, 1)) 



7 
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and 



p(2) _ / p (o) 9 



7 



ER2 



So, to second order, we have 



E Nl =h\uj + g 



E N2 = h { u - g - — 



Remark 6: Note that it was the non-rotating wave approximation part of the inter- 
action that gives the second order corrections (ab term). This is an example 
of a non-energy conserving, or virtual transition . 



or 



Energy-Conserving Transition 
ot6|0,l) = |1,0) 



"b" 



Non-Energy Conserving Transition 

a t S^(|o,i>±|i,o» = ^(|o,i>±|i ) o» 



"b" 



Only the intermediate step is non-conserving, overall, energy is conserved. 



13.3 Time Dependent Perturbation Theory 

For this analysis, we are going to consider the unitary propagator (aka time-evolution operator): 

U(t,t ) = e- i ^*-*°)/ R (1) 
which gives the solution to the Schrodinger equation: 

m)) = u(t,t )m )) (2) 

The propagator itself obeys the equation of motion. 

^-U(t,t ) = ^HU(t,t ) (3) 
at in 
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with the initial condition that U(to,to) = 1. Formally integrating (3) gives 



1 



U(t,t ) = l + - / HU(r,t )dr (4) 

Time dependent perturbation theory is based on (4). Our 1st goal is to rewrite (4) so that the 
dynamics generated by the unperturbed interacting part of the Hamiltonian are separate. 

H{t) = H (t) + V(t) (5) 

Here we are allowing Hq and V to explicitly depend on time. 

Remark 1: We will use the interaction picture to separate the perturbed and unperturbed dy- 
namics. Let Uo(t,to) represent the propagator which governs the unperturbed motion such 
that Uo(t,to) satisfies 

U(t,t ) = ^-H (t)U (t,t ) (6) 
in 

with initial condition Uo(to,to) = 1. Now, Uo(t,to) can be thought of as a "rotation" in 
Hilbert space. If we go into a reference frame that rotates at a rate determined by Uo(t, to), 
then the state only evolves in that frame if V(t) ^ 0. In the case that V(t) ^ 0, we are in 
the interaction picture. 

Remark 2: Alternatively, a non-rotating frame in which the state evolves according to both Hq 
and Vq gives the usual Schrodinger picture, while a frame that rotates at a rote determined 
by both Hq and Vq (the state is stationary in this frame) gives the Heisenberg picture. 

Remark 3: The states in the Schrodinger and interaction pictures are related by 

\Mt)) =u (tM\Mt)) (7) 

Schrodinger Interaction 
picture picture 

so that if V(t) = 0, \ipi(t)) is stationary, and \ipi(t)) = \ipi(t )) = |V>s(*o))- If V{t) / 0, 
then the state in the interaction picture evolves according to 

\4>i(t)) = u I (t,t )\Mto)) (8) 

where Ui{t,to) is the interaction picture evolution operator. 

Remark 4: The evolution equation for Ui(t,to) can be obtained via (7) and (8): 

\Mt)) = ^o(Mo)#/(Mo) |^r(io)> 

= U (t,to)Ui(t,to)\Mto)) (9) 

Comparing to (2): 

=>Us(t,t ) = U (t,to)U I (t,to) (10) 
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Substituting(lO) and (5) into (3): 



at 



at j 



\h Uo] U: + 4 VUoU! 
in I in 



dt ih 

since Uq is unitary (and thus UqUq = I). 
Defining 



TT" U VU t 



o) Uj 



we have that 



Formally integrating gives 



V I (t) = U J (t,t )V(t)U (t,t ) 



-^/(Mo) = ^Wi(*,<o) 



1 /■* . . 
E7> (t,«o) = 1 + - / Vr?7/(r,to)dT 



Time dependent perturbation corrections are obtained by iterating (14): 

1 



1 r* - 

Ui(t,to) = 1 + / Vi{t)cLt 



1 + 



ih 



y/(T , )^/(r , ,to)dr / 



= 1 + k, f ^ /(T) dT ~¥ I dT f dT> ^(r)^(r / )c>/(T , ,to) 

J to J tQ J t() 

Higher order corrections can be obtained via an infinite number of iterations. In powers of 
this becomes: 



n=l 



Ujit, to) = l + J2[j%) j dTn [J dr "-i • • • P drxV^Viirn-i) ■ ■ ■ Vi(n) 



to 



where the order of integration is such that 

t < n < t 2 < ■ ■ ■ < T n _i < T n < t 



Remark 5: We now transform (16) back into the Schrodinger picture using 

U s (t,t ) = U {t,t )U I {t,t ) 

and 



Also, note that 



Uo(Tj,t )U^(Tj-l,t ) = U (Tj,Tj-l) 
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Thus we have that 



U S (t, to) = U (t, t ) + ^ i^ — j J dT n J dTn-l ■■■J dn 



■ Uo{t,T n )V{T n )Uo{T n ,T n ^) V(T n _ 1 )---Uo(T 2 ,T 1 )V(T 1 )Uo(T 1 ,t ) (17) 

V v ' 

u (t,t ) Ul(r m t )V(T n )U (T n , to) ul{rn-i,t ) 

V v ' 

Remark 6: Starting at the earliest times (right) and moving to interacting (left) in the integrand 
represents 

1. Free evolution from to to n 

2. Interaction at t\ 

3. Free evolution from t\ to t 2 

4. Interaction at r n 

5. Free evolution from r n , to t 

So that the full evolution represents a sum over all elementary time evolutions, taking into account 
all possible numbers of interactions between to and t. 

13.3.1 Transition Properties Between Energy States 

For this analysis, consider the case where the free Hamiltonian Hq is time independent. 
Then 

U (Tj,Tj-i) = e-^fo-^-OM (18) 

and we consider a system with discrete energies Ei, E2, We develop the perturba- 
tion theory in terms of jumps between the levels. 



Ek+i 



E k uJ k ,l = 



E k - 



The transition probability between some initial state \E a ) and some final state \Eb) is 



defined as 

2 



W, 



a— *b 



E, 



U s (t,t )E a ) (19) 
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Remark 1: The transition probability depends on the interval to —> t during which 
the transition takes place. We expect non-zero transition probabilities if there 
is an interaction V(t) which couples the initial and final states \E a ) , \E^). The 
interaction is specified in terms of matrix elements in the unperturbed energy 
basis: 



V{r)=yVki{r)\E k ) {E l 



k,l 



where Vki(r) is: 



V u (t) 



E r 



V(t] 



Ei 



(20) 



(21) 



Substituting (17), (18) and (20) into (19) gives: 



1 



Vk n ,i n (T n ) (Ei n 



to J to 



-iHo(t-T n )/h 



E kr. 



-iHo(T„—T n -i) /h 



E k n -i 



-iHo(i-i—to)/H 



■^-.^.(Vi) (Ei^ | • • • e~ lH ^-^ h \E kl ) V^n) (Et 

E(i) ^n-i-/ *i-EE-E^ (Mw)/ 'w 



/ dr n 

n=l v 7 17 *0 



£n 'n— 1 ^2 



t )/h 



(22) 



Remark 2: 

• ^^Lx — ► sum over the number of intermediate state interactions between 
to and t 

• — > integration over all possible times for intermediate state interactions 

• Vi ji i,_ 1 — > strength of interaction connecting intermediate states 
and \E tj ) 

• (• • • ) — y transition amplitude is sum over amplitudes for transitions tak- 
ing place via a set of intermediate or virtual states \Ei 2 ) , \E[. 3 ) , . . . , \Ei n ). 
There are n — 1 intermediate states (stepping stones) so that there are n 
steps between \E a ) and \Eb). 

Remark 3: Diagrammatically, the transitions can be represented as: 
where the vertices are where the state of the system changes. 
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Tl 




Remark 4- To first order in perturbation theory, (22) looks like 

2 



Wn 



where oj ba 



1 

l 

~ ¥ 

Eh—E a 



-iE b (t-T)/h 



to 



V ba (r)e 

2 



-iE a {r-t )/h 



dr 



e iuJbaT V ba (r) dr 



to 



(23) 



13.3.2 Time-Independent Interactions 

The time integral in (23) is easy if the interaction is independent of time: 



W a ^ b = ^\V h 



ft 2 ' lOJ ba 



iu) ba (t-t ) _ e -iui ba (t-t ) 



ba 

2 2 1 

TolVfeal — (1 - COS [UJ ba {t - t ))) 



Let At = t — to and define 



2(1 - cos(u;At)) 



(24) 



(25) 



Remark 1: f(u, At) contains information about energy conservation during the tran- 
sition. It is peaked at lo = 0, with a maximum /(0, At) = (At) 2 . This is seen 
in Fig 13.1. 

In the final limit At — ► big, this function becomes highly peaked about lo = with 
half-width Aw ^ and 



At) ► A5(lo) 



(26) 
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Figure 13.1: Plot of /(w, At). 



where A is found by integrating (25) over u. 



f°° l-cos(wAt) 
A = 2 / £ - dco 

= 2nAt (27) 

where the integral is done by contour integration. (Let f(z) = 1 ~ e z 2 — and the 
integral is the real part of this along the x-axis.) So, in the limit of infinite interaction 
time, the transition probability is unless 



Eb — Eg 

h 







(28) 



This is a statement of conservation of energy: A mismatch of energies is tolerated for 
finite time intervals as long as 



AEAt « ttH, 



with Ah = TvjJ ~ - — 
At, 



(29) 



Transitions occurring in a finite time At involve an uncertainty about the energies. 
These transitions can occur if E a and are close enough together so that the total 
system energy uncertainty can bring them together during a finite interaction. Note 
that AE in (29) is the difference in unperturbed energies, not the energy shift due to 
the perturbation. 



Remark 2: For transitions between states of same energy: 

W a ^ b = ^\V ba \ 2 (At) 2 



(30) 



which has limited validity since W a ^b must not exceed 1, even as At — > oo. 
Hence (30) is valid only if 



At < 



h 

VbA 
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13.3.3 Transition to a Continuum of States 

Often, \Eb) is not a discrete state, but is part of a continuous spectrum: 

\E b ) \E b ±6E 



\E a ) 



In these cases, we must find the transition probability to the state whose energy falls 
in the range of energies centered at E b . Let B denote the subspace of final states with 
energies in the range Eb — 5E to Eb + 5E and define 



W, 



a—*b 



(31) 



all states 
in B 



where W a —,b is given by (24). For a continuum of states, we replace the sum in (31) 
with an integral which includes the density of states: 



pEb ■ dEb = Number of states with energy in the range Eb to Eb ± 5E 
So (31) becomes 



rE +8E 9 

W a ^ B = / dE p(E)~2 

JEo-SE h 



E 



V 



E, 



1 - cos At) 



(32) 



(33) 



Remark 1: If the range of final energies is sufficiently small, then 

p(E) « p(E b ) 



E 



V 



in which case 



W a ^ B = 2p{E)\VbaV I dE 
'Eo-SE 



E +8E 



V 



E a 



= \VbaV 

1 - cos At) 
(E-E a f 



(34) 



Remark 2: From §13.3.2.R1, we know that the function in (34) has a max value 
°f (if) a ^ E = E a , then secondary maxima get smaller as deviations from 
the resonant energies (E = E a ) decrease. The width of the central maximum 
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is ~ tt, and we assume that the integration range is much larger than this 
width. 



SE > — 
At 



(35) 



With this approximation, we consider (34) for 2 distinct cases: 

1. | E b — E a \ > 5E implies that the final states are non-conserving by an 
amount larger than the width of the transition probability maximum 



\E b ) y Integration Range 



l«„) 



l-cos(g^At) 
(E—E a ) 2 



In this case, 



So that 



And so 



{E-E a y E b -E a y 

Eq+SE 



dE 



E Q -SE 



1 — cos 



E-E n 



-At 



28E^>-hAt 



it ; 



Ap(E b )5E\V ba \< 
{E b - E a f 



(36) 



2. \E b — E a \ < 5E implies that the transition probability peak falls within 
the range of final energies. This means that the transition is energy 
conserving. Now, 




Eo+SE 
'Eo-SE 

and (34) in given by (27). Thus 



in (34) 



W, 



2irAt 

h 



p(E b ) \V ba \' 



(37) 
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Remark 3: Eqn (37) is an important result since it represents the dynamical behavior 
in many quantum systems. I.e., imagine N a systems in state \E a ) at time t. 
A short time At later, the number of systems in state \E a ) is 



So 



N a (t + At) = N a (t) 
N a (t + At) - N a (t) 



At 



l-^p(E b )\V ba \ 2 At 



—p(E b )\V ba \ 2 N a (t) 
n 



N a = -TN a 



where 



W, 



a~>B 



At 



YP( E b)\ V ba\ 2 



(38) 



(39) 



and thus T is a transition rate (transition probability per unit time), and (38) 
describes an exponential decay at this rate. 



N a (t) = N a (0)e 



-rt 



(40) 



Examples of where this holds: /3-decay, spontaneous emission from an excited 
state of an atom or molecule. 



13.3.4 Periodic Perturbations 

Consider a harmonic interaction of the form: 

V{t) = A^e~ iujt + Ae iujt 
with A an unspecified time-independent operator. Then (23) becomes 



W, 



a—>b 



1 

¥ 



t +At 



dT 



to 



E b 



it 



E a )e 



J(w ba -Ul)T 



+ (E b 



A 



E a )e 



(41) 



And the integral is largest when uj ba — to = or uj ba + uj = 0. These are resonance 
conditions corresponding to energy conservation. 



E b = E a + huj 
E b = E a - hio 



(42) 



Physically, these represent an energy quanta absorbed or emitted (respectively) from 
the system providing the periodic potential, as seen in Fig. 13.2. 



Remark 1 : We require that w>^so that the frequency is much greater than the 



width of the transition resonance peak. This establishes that the system has enough 
time to "feel" the periodicity of the perturbation. 
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m 

A 

IK) 

(a) Resonance and 
energy conservation 
with absorption of a 
quanta 



m 

— \Ea) 

(b) Resonance and 
energy conservation 
with emission of a 
quanta 



Figure 13.2: Resonances of Periodic Perturbations 



Following through with the integration in (41) gives us: 



W a . 



2 

¥ 



it 



E a 



2 1 - COS ((Ufa - w)At) 2 



(uJ ba - to) 2 



+ 



h 2 



E b 



A 



E a 



2 1 — COS ((UJba — Uj)At) 
(U ba - Lo) 2 



Considering transitions to a continuous band of final states gives energy conserving 
transition rates: 

1. For the case where E b + 5E > E a + hw > E b — 5E: 



W, 



a—*b 



2vr 



p(E a + hw) 



E a + flLU 



At h 

2. For the case where E b + 5E > E a — hw > E b — 5E: 
W, 



it 



'a^b 27T 

AT = H p{Ea ~ ^ 



E a — hw 



it 



E a 



E n 



(43) 



(44) 



Remark 2: Equations (43) and (44) give transition rates that are often referred to 
as Fermi's Golden Rules . The hw supplied by the perturbation provides the 
energy to make the transition energy conserving. 



13.3.5 Sudden and Adiabatic Perturbations 

Two situations in which the response of the quantum systems changes is easy to 
calculate is: 

1. Sudden Approximation: Here we assume the system can not change it's 
state arbitrarily fast, so a sudden change in external field will not change 
the state of the quantum system during the field change. 

2. Adiabatic Perturbation : When perturbation is slow enough that the state 
of the system changes adiabatically with applied perturbations. 
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Chapter 14 

Scattering 



14.1 Scattering Cross Sections 

Consider a general scattering experiment: 



Detector has solid angle df2 



Beam of 

Monchromatic 

Particles 




Scattering Sample 
containing N Scatterers 

Where the particle flux J is the number of particles crossing unit area per unit time. Also, the 
beam is weak enough so that the particles scatter independently (no coherence). 

Remark 1: We neglect coherence of scattered waves between scatterers-each scatterer acts as if 
it were alone. This assumption fails for x-ray diffraction. 

We define: 

N = # of particles recorded at the detector per unit time 

^N=JNa(Q)dQ (1) 

where a(Q) is the differential scattering cross section, and is a parameter that characterizes the 
collision between the incident particles and scatterers. The "differential" refers to the detector only 
monitoring a small solid angle d£l. 

Remark 2: The total scattering cross section is 

a tot = fa(Q) dQ (2) 
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Remark 3: Though the target can be quite complicated (an atom made up of electrons and 
protons, etc.), it is often sufficient to describe the scatterer by a central potential V(r). 
(if we consider elastic collions only!) This is what we will do in this section. 

Remark 4 : The goal is to relate what is measured in experiments (the differential cross section), 
to a solution to the Schrodinger equation. 

Assume incoming particles have energy E, and momentum hk. Then the stationary solution must 
satisfy 

r fi 2 

-V 2 + V(r) 



2m 



ipk = Eip k (3) 



where V(r) defines the interaction between particles and scatterers, and E is the energy of the 
incident and outgoing particles. We want a solution that describes the scattering process, so we 
look for a solution of the form: 

Mr,t) = e~ lEt / h Mr) (4) 



with asymptotic form 



p ikr 

fc(r)~ + f{n) — (5) 

incoming 



particles outgoing spherical 
wave of scattered 
particles 



as r — > oo. This means that ijjk(r,t) gives an expanding spherical wave. 

Remark 5: The e ik v represents a beam of particles with speed ^ and a density of 1 particle per 
unit volume. The corresponding particle flux is 



j = i.«* 

m 



The flux of particles entering the detector at angle £1(0, (j>) is: 

(number of outgoing \ fik 
particles per unit ( De a t r e e c * or ) 
volume / 

For large distances r, the first term is and the detector area is r 2 dCl. 

hk 

N = \f(n)\ 2 —dn (6) 



??? 



Comparing (6) with (1) gives: 

a(0) = \m)\- 



measured calculated 



Remark 6: Now the main task is to calculate what is called the scattering amplitude f(£l). 
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Remark 7: Often we are considering the scattering of individual particles off one another. In 
this case, an incoming singe particle would be described by a time independent wave-packet 
rather than a plane wave, and the picture would look like: where the big particle is being 



2d 




scattered by the smaller. The incident particle (wavepacket) has velocity v = ^ = 
The longitudial and transverse dimensions of the wavepacket are determined by the spread 
in momentum. 

In the above picture we've made the approximations that: 

• l,d 3> a - otherwise the scattering depends on the details of the shape of the incoming 
wavepacket. 

• I,d> A - otherwise the wavepacket description of the incoming particle will not give 
well defined energy and momentum to be small compared to the scale on which the 
scattering properties change with energy and momentum. 

• d > b - for non-trivial results (if d <C b, then there is no significant perturbation, 
assuming short range potentials that fall off faster than ~). 

• D S> a, A - so that the influence of the potential on the detector has dropped to 
(detector is asymptotic regime). 

• D sin 9 S> d - so as to separate the transmitted and scattered parts of the incident 
wave packet: 




Mathematically, for the description we would multiply the incoming plane wave by a real envelope 
function which depends on the distance from the scatterer (more precisely, on the impact parame- 
ter). We would eventually derive a description of the scattered wave packet that depends on the 
solution to the Schrodinger equation in a central potential. Given this description and the above 
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constraints, we can show that the differential cross section (which depends on the probability for a 
wavepacket to scatter within d£l of f2 (location of the detector)) is still given by 

a(u) = \f k (Q)\ 2 (7) 

where the asymptotic solution to the Schrodinger equation is: 

V>fc(r)~e <fc - r + /*(fi)^ (5) 

14.2 Waves and Phase Shifts 

Our goal now is to find f(£l) for a given central potential. Recall that in §10.2.2 and §8.2.3 we 
found that the solution to the spherically symmetric Schrodinger equation could be written 



^ r ,9,cj>) = Y J A l J^-Y l 



lm\ 

r 

l,m 

where Uei is a solution to the radial equation (§10.2.2.E2) and Yi m (9, 0) are the spherical harmonics. 

Remark 1: The scattering problem has azimuthal symmetry. For initial propagation direction, 
aligned with the z-axis, the solution is (^—independent, and only m = terms are involved. 
Thus we have 

oo 

V(r,0)=5>-^PKcos0) (8) 



where Pi{cos9) is the Legendre polynomial and ai = J ^^-A^q. 
Remark 2: The goal is to match (8) to the asymptotic form (5) 

pikr 

VM) >e ik - r + f(6) 

r— >oo f 

where we are using arrows to represent asymptotically similar. 
If we take k in the z-direction, the incident plane wave can also be expressed as: 

oo 

e lk r = + !)(*)' ■ ji(fcr)fl(cos0) (9) 

(See §10.2.3.R11). Expanding f(6) in Legendre polynomials: 

/(0) = X>iK«>s0) (10) 

/=0 
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We require Uei(t) to be a regular solution to the radial equation, and asymptotically, we can write 
it as 

U E i(r) ► siii(Ax--J+<5,) (11) 



where 5i is called a phase shift. The asymptotic form of (8) is: 

2i r 2i r 



r— >oo » 



Z=0 



P;(cOSl 



(12) 



In order to match this to the asymptotic form (5), we write the asymptotic form of ji(kr) in (9): 



1=0 



ifcr 



P/(cosi 



(13) 



Now, (12) and (13) must match as there are 2 different forms of the same incoming and outgoing 
spherical waves. This requirement gives an expression for o;: 
• Incoming Waves: 



2i 



, v ( —i l 

^ 2l + l h lii 



(21 + 1) 



ai = — — -i l e iSl 



(14) 



Outgoing Waves: (matching coefficients) 

HIj 5i \ _ (2i + iy (-»)' 



ai{ — e 



k 



2i 



+ fi 



(21 + l)e iSl e iSl (21 + 1) 



k 



fi 



2i 2ik 
- 1 



+ fi 



(1 + 1) r^. 

2ifc 



And the expression for the scattering amplitude in terms of the phase shifts is: 

1 oo 



(15) 



(16) 



1=0 



Remark 3: Solving the differential scattering cross-section is reduced to finding the phase shift 
for the different angular momentum components of the scattered wave. 



Remark 4 : Note that we can write the asymptotic solution (13) with a/ from (14). 
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where this phase factor shifts the phase of the outgoing partial wave relative to the incident 
plane wave-this is where the entire scattering process is captured. (Partial waves refer to 
individual terms in the spherical wave expansion.) 



Remark 5: The differential scattering cross-section is terms of phase shifts is: 

<r{0) = l/WI 2 = p J> Z + l)(2/' + l)e^- 5 <') sm(5 l )sm(5 l/ )P ll (co S e)P l (co S 9) (17) 



IM 



Integrating this over a unit sphere gives the total scattering cross-section. Noting the 
orthogonality of the Legendre polynomials: 



4tt 



a *°t = -p- 5^(2Z + l)sin 2 (<f z 



(18) 



1=0 



We can write 



1=0 



with 



47T 



a x = ^-(2/ + l)sin 2 (^) 



(19) 



and G\ gives the contribution to scattering cross-section from partial waves with angular momentum 
I. 

Remark 6: The partial wave expansion is most useful when only a few terms in (18) contribute. 
To get an idea of which terms dominate, consider a scattering (classical) problem. Here 




2r 



V(r) = for r > vq. Now the angular momentum of the incoming particle is L = bp 
and is a constant of motion. If b > tq (L > rop), then there is no scattering. Quantum 
mechanically, L — > HI (more precisely, L — > h?l(l + 1)) and p = Hk. Thus, we expect the 
strongest scattering when 

I < kr (20) 

and the weakest scattering when 

I > kr (21) 

Generally speaking, we expect ai to become unimportant for I > kr^. This holds if V{r) = 
for r > ro, but turns out to be true also if V(r) is very small (though not strictly 0) at some 
range tq. 



14.3. EXAMPLE: SCATTERING FROM A HARD SPHERE 



215 



Remark 7: Equation (8) (weak scattering) holds when there is either low energy (k — > 0), or high 
angular momentum (I — > oo). Recall that high angular momentum particles will "miss" the 
target if L > rop. Physically, lower energy corresponds to large wavelength. In the limit 
tt- = <C 1, a small scatterer will have little effect on a large incident wave. (Imagine a 
small obstacle in sound and water waves.) Any scattering that does occur will be isotropic 
since large wavelength waves cannot probe the spacial structure of the scatterer. 



14.3 Example: Scattering from a Hard Sphere 

Consider the potential: 



V(r) 



oo, < r < ro 
0, ro < r < oo 




The general solution to the radial equation (§10.2.3) is 

Ri(r) = Aij^kr) + Bi m (kr) 
where Ai,B[ are constants and the BC's are 

Aiji(r k) + Bir)i(kr ) = 
3i{kr 



B, 



■Ai- 



Vl(kro) 

R i( r ) = A {vi{kr )ji(kr) - ji(kr )r]i{kr)) 



where A', 



1 ~ m{kro) ' asymptotic form of the solution (§10.3. 1.R6) is 

n ( \ Vi(kr ) . , l . ji(kr ) , l , 
Rl\ r ) > — — — sm \kr — *y ) H — — cos [kr — 



>o kr 

Comparing this with the general form (§14. 2. Ell): 
Rl(r) > 7^ sin (kr - l f + 5 t ) 



kr 



— sin (kr — ^) cos 5i + — cos (kr — sin <5/ 

kr v A kr A 
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Vi(kr ) = cos (5/ and ji(kr ) = sin<5/ 
ji(kr ) 



tan Si 



5 1 = arctan 



Vl( kr o)J 



sin 2 (V - jf{kro) 



r)f(kr ) + jf(kr ) 

Remark 1: At low energies, the spherical Bessel functions become 

(kr) 1 



jl {kr) 



Vi(kr) - n 



r^o (2Z + 1)!! 
(21 - 1)!! 

r-+o (kr) l+l 



So that the partial cross section for scattering off a hard sphere is 

4vr (21 + l)(kr ) 4l+2 



k- 



,o k 2 [(2Z + 1)!!(2Z- 1)!! 



2 



The dominant cross section is s-wave scattering (1=0): 

ai > T^(kr ) 2 = ^r\ 

k^O k 2 

where ro is the length scale characterizing the range of the scatterer, and is sometime called 
the scattering length. 

Remark 2: At high energies, we'll use the asymptotic forms of the Bessel functions (see Messiah) 
to get: 

a(Q) - ► \rl (1 + cot 2 (|) J? (kr smO)) (f) 

k — >oo 4 



and 



k — >oo 



This limit corresponds to very short wavelengths, from which we might expect the classical 
results: 

a(Q) = \r 2 \ 

4 > Classical results 

otat = nrl J 

We do not get these results since the wave nature of the scattering process does not com- 
pletely vanish. This is because the hard sphere has "sharp" edges, so that any wavelength, 
no matter how short, is affected by the discontinuity and there are diffraction effects. The 
Bessel function in (f ) is the diffraction term. 
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14.4 Scattering from a Square Well — Resonances 

Consider the potential 



V(r) 



-V Q , r <r 
0, r > ro 



r 



-V 



Remark 1: From before we know that a ID square well exhibits resonances which are very sharp 
when the well is very deep and incident particles have low energies. Similar resonances 
occur when scattering from a 3D well. 

Remark 2: Recall from §10.3.2 that regular solutions inside the square well are 

Rl(r) = ji(nr) (r < r ) 

with 



2m 



(E + V 



and we have chosen to normalize Ri{r) for convenience. Outside of the well, the solution is 
a superposition of the Hankel functions. 



R l (r)=A l h { l + \kr) + B l h 



,(-)/ 



(r > r ) 



with 



k 



2mE 



These solutions correspond to E > since we are interested in the unbound states outside 
the well for scattering phenomena. The boundary conditions require continuity at r = tq: 



ji(Kr ) =A l h ( l +) (kr )+B,h ( i \ 



cr J 



and continuity of the derivatives at r = tq: 



Kjl( K r ) = k Ath^'ikro) + Bihf(kro) 



(-)'/ 



(1) 



(2) 
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Remark 3: We also want to match the asymptotic form of the solution outside the well to the 
scattering solution: 



Ri(r) ► 7^ sin (kr - | + 5 t ) 

r-»oo kr 

1 1 



+oo kr 2i 



e i(kr-lir/2) e i8 l 



e -i(kr-lTr/2) e iSi 



Inserting the asymptotic form of the Hankel functions: 

1 



Rl(r) 



kr 



Aie i(kr-ln/2) + Bie -i(kr-lTr/2) 



Comparing the forms gives 



-id, 



A l = e l0 \ B l = e 
Bi 

So we need only the ratio of the unknown coefficients. From (1) then, 

Ai _ 1 ji{nr Q _ h\~\kr ) 
B i~ B lh\ +) (kr ) h^\kr Q ) 

Now, taking (1) • kh\ +) '(kr ) - (2) • /i[ +) (fcr ): 



fcjiCwo)^ (+) '(fcro) - Kjl(Kr )hi +) (kr ) = A t k \h\ +) (kr )h\ +) ' (kr ) - h\ +) ' (kr )h\ + \kr ) 



+ 



+ B t k 
Bik 



h^\kr )h^'\kr ) - h\~>'(kr Q )hf ! (kr Q ) 



,(+)' 



(-)'/ 



,(+), 



2i \ 
kri) 

where the last term is from the Wronskian relation of Hankel functions (Messiah X.38). 



D 2i&i 



h\ \kr ) ji(nr ) 



k2i 



h\ +) (kr ) h\ + \kr ) (kr ) 2 kj^h^ \kr Q ) - ^[{nr^ (Ar ) 

h\-\kr ) 
h\ + \kr ) 



k 1 
1 - 2i- 



-(M 2 \^'(kr )-^h^(kr ) 



Remark 4 : For a deep potential and moderately low energies: 

fiXo S> I and k S> k 
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So 



e 



2iS, _ h l ( kr o) _ l-iji(kr )/m(kr ) 



h\ + \kr ) 1 + iji(kr )/vi(kr 

]_ _ / jljkrg) \ 2 _ 2iji(kr ) 



1 + 

= cos(2(5/) + ism(25i) 

=>- tan d z « — - 

^(^o) 

which is the same result as for the hard sphere. 



mWoj ) 



Remark 5: For certain energies, the denominator in the second term of (3) vanishes. In this case, 
note that 



jl(Kr ) Kr »Z 

For simplicities sake, let kro S> Z so that 

h [ r\kr ) 



cot (kto — lit/ 2) 



Thus (3) becomes 



where 



2iS, 



/i[ +) (A:ro) kr °» l 



» e 



-2i(kr ~lw/2) 



(kr ) 2 h { ~\kr )h^'(kr 



(+)/ 



{kr fh\ ) (kr )h^' , (kr ) ► * 



(+)'/ 



-2i(kr ~lw/2) 



2i^ 



if - cot ( K r - 



-2i(fer- -«7r/2) 



W2) cot(/cr -Zvr/2) + i* 



cot(Kr - ln/2) - i\ 
Si = -{kr - Ztt/2) + 



tan i 



cot(Kro — lir/2) 



Now ^ is generally negligible when f <C 1, unless Kro — Zvr/2 = (2m + l)n/2, where m is an 
integer. When we have this secondary case, we have resonances. 



220 



CHAPTER 14. SCATTERING 



The distance between resonances is Aktq = n. 



2m 
~W 



(E + AE + Vq) 



2m 



(E + Vq) 



r = 7r 



m AE 
1 + — — - 1 



h 2 k 2 



r = 7r 



Assuming K^k^Vo^E^K, 2 



AE 



2m\r — Jl. 

W v ° = K o- 



ir 2 h 2 K 
mro 



AE 



2nV 
mro 



Since 



7 

o 




















KT 



(n + l)ir (n + 2)7r 




where in the range (*), tanc/> goes to — > oo and — oo — > which corresponds to a change of n in 
0. So (5; = —(kro — In/ 2) + <f> increases by n each resonance. In between resonances, background 
scattering is determined by 5i = —{kr$ — In/ 2) + (p. 



Chapter 15 

Relativistic Quantum Mechanics: The 
Dirac Equation 



The Dirac equation is a relativistic description of a quantum system that includes the relativistic 
nature of some systems (like high-Z atoms) as well as spin and spin-orbit coupling. 

15.1 Spinless Free Particles - The Klein-Gordan Equation 



Recall that in deriving the Schrodinger equation, we started with the kinetic energy: 

p 2 

H = — 

2m 

and let 

> => ih 

"dt. 



f) > => ih— lib) 



In deriving a relativistic description, let 



tt2 2 2, 2 4 

H = c p + m c 



So that the substitution of operators gives 

|^> = -^(c¥ + ™V)h/>> (t) 
This is the Klein-Gordan equation. 

Remark 1: Note that we considered H 2 rather than H = c 2 p 2 + m 2 c 4 . This is because the 
former gives a symmetric description of space and time as desired by relativity. Otherwise, 
expanding H in a momentum basis: 

1 P 2 1 P A 
H = mc ( 1 + 77 — 2~2" ~ + 
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gives even and odd derivatives not on an "equal footing" . 
Remark 2: (f) is a good description for spinless particles, but is insufficient for particles with 



15.2 Dirac Equation for Free Particle 

The Dirac equation for a free particle is given by 

ih-r^ \i>) = {col ■ p + (3mc 2 ^ \ip) (1) 

where 

*-(!!)■ Mo-;) 

where 

• & is the Pauli spin matrices 

• I is the 2x2 identity 

• c is the speed of light 

Remark 1: In the derivation of the Klein-Gordan equation, we took H -> H 2 to get the space-time 
symmetry. Here, Dirac assumed that y c 2 p 2 + m 2 c 4 oc p+ something without momentum 
to get space-time symmetry in an equation of 1st order in space and time. 

2 2, 24 - , n 2 

c p + m c = col - p + pmc 
where a, (3 are determined by matching coefficients: 

c 2 {pI +P 2 y + Pi) + m2 ° A = c 2 (a- p) 2 + mc 3 (a ■ pj3 + /3a ■ pj + /3 2 m 2 c 4 

Matching coefficients gives the following relations: 

1. a 2 = p 2 = I where i = x,y, z 

2. + aj&i = OLj\, = for % j 



3. P&i + ai/3 







Remark 2: From these relationships, we conclude that a, (3 are traceless Hermitian matrices with 
eigenvalues = ±1. 



Remark 3: Traceless with eigenvalues=±l even dimensional. But they cannot be the 2x2 Pauli 
matrices (since there are 4 matrices here, and there are only 3 Pauli matrices with these 
properties) => try 4x4 matrices 
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Remark 4 : The choice of a and (3 given in (2) satisfy these requirements, though it is not unique 
(since the properties of Remark 2 are preserved under a unitary transformation). 

Remark 5: in (1) is called the Lorentz Spinor and has 4 components (to be compatible with 
4x4 ct's and /3's). 

Remark 6: {ip\ip) = a constant since the Hamiltonian is Hermitian. 



15.3 Electromagnetic Interactions 

The classical Hamiltonian for a particle with charge q is: 

ff=[(p-2A)c 2 + m V] 1/2 + # 

where $ and A are the electrodynamic scalar and vector potentials. Generalizing the Dirac equation 
so that p — ► p — - A and adding the electromagnetic potential energy: 



..dip 



ca - (p- 2 A) + /3mc 2 + q& 







(1) 



15.3.1 Electron Spin and Magnetic Moment 



To see how spin and magnetic moment drop out of the Dirac equation, it is sufficient to 
let = and keep only O(^) 2 terms. Let ip(t) = ip e - lEt / h to find energy eigenstates. 
Plugging this into (1): 



ca. ■ (p — | A) + /3mc 2 



where 



Etp = [cot • II + I3mc z ) i> 
fl = p-±A 



(2) 
(3) 



Remark 1: II is called the kinetic momentum operator. 

Remark 2: Recall that ip is a 4 component spinor and that a, j3 are given by 
a = a x (g) a and j3 = a z ® I 
which suggests that we can write 



X 



(4) 
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X 







where y an d <p are 2 component spinors. 
Plugging (4) into (2): 

_, [y] [ mc 2 I ccr ■ I 
E = ^ ~ 9 
<p\ [ccr • II —mc 

~(E - mc 2 )i -co ■ TV 
-ccr ft (E + mc 2 )/ 
which gives a set of coupled equations for y and (/?: 

(i£ — mc 2 )y = co" • n<^ 

(£/ + mc 2 )tp = ca ■ Tlx 

From (6): 





X 




0" 




JP. 








co- n 



+ mc 2 

Remark 3: In the non-relativistic limit, 



X 



(5) 
(6) 

(7) 







co • {p - \ A) 


c(mv) 1 v 






.E7 S + mc 2 + mc 2 


2rnc 2 2 c 



< 1 



where 

• E s is the Schrodinger energy <C mc 2 

• mv is the typical momentum 

So in the non-relativistic limit, ip is the small component and y is the big component. 

o fty 



2mc 



In the non-relativistic limit, (5) becomes: 
E sX = c& ■ flip = — (o ■ n 



(8) 



(9) 



— (n • n -Ho • n x n) 

2m V / 



X 



and n x n = ^s. So now (9) becomes 
2m c 7 2mc 



X = E sX 



(10) 



which describes a spin ^ particle in a magnetic field. The energy splitting (Zeemann 
effect) gives 

mc \2 / 
and the picture in Fig. 15.1: 
Here 2/j>bB = 25E = g^sB g = 2s + 1, and we see that electron spin and magnetic 
moment are accounted for in the Dirac equation. 



15.4. THE DIRAC EQUATION FOR THE HYDROGEN ATOM 



225 



HbB 

-f B B 

Figure 15.1: Zeemann splitting 



15.4 The Dirac Equation for the Hydrogen Atom 

For the hydrogen atom: 



e 2 

V = e$ = 

r 



and the Dirac equation is 



dip 



ca. ■ (ft - 1 A) + f3mc 2 - - 







We still decompose ip into the 2 component spinors and obtain equations for x an d f (§15.3 (5) 
and (6) with E -> E - V). 

(E — V — mc 2 )x = ca ■ ftp (1) 
(E - V + mc 2 )ip = ca ■ ftx (2) 

Remark 1: In absence of a B field, we can let A = so II — ► ft, and hence the above results in 
(1) and (2). 

From (2), we have 

tp = (E- V + mc 2 y l ca -ftx (3) 
which is substituted into (1), giving: 

(E -V + mc 2 )x = ca ■ ft(E - V + mc 2 )~ l ca ■ ftx (4) 

Remark 2: Note that we were careful to keep the order correct, as ft can operate of V, so order 
does matter. 

Remark 3: Solving (4) and plugging into (3) gives the solution to the full Lorentz spinor ip. 
Remark 4- Keep 0(|) 2 , (4) becomes: 

„2\-l a 



(E s - V)x = ca ■ ft(E_ - V + 2mc i Y 1 ca ■ ftx 
a- ft) {a- ft) x 



c 2 



2mc 2 
P X 



2m 

which is just the non-relativistic Schrodinger equation. 
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15.4.1 Fine Structure 

Recall that in solving the Schrodinger equation for the hydrogen atom, we resorted to 
perturbation theory to get the fine structure of the energy eigenspectrum (previous 
HW). The time structure was due to relativistic effec8ts which emerge naturally from 
the Dirac equation. Expanding (4) to 0(^) 4 (one order higher than that required 
to recover the Schrodinger equation), and letting E = E s + mc 2 in (4) with E s the 
solution for the Schrodinger equation gives: 



(E s - V + 2mc 2 )" 1 



1 



2mc 2 
1 

2mc 2 
1 



1 + 



E s -V 

2mc 2 
E K -V 



2mc 2 
E„-V 



2mc 2 4m 2 c 4 



And (4) becomes 



(15.1) 
(15.2) 
(15.3) 



{E a - V)x = ca p 



2mc 2 



E s -V 
4m 2 c 4 



ccrpx 



EsX 



p- 



X 



EsX 



p 2 

2m 



a ■ p(E s - V)(a ■ p) 



4m 2 c 4 



extra piece giving fine structure 



X 



(5) 



Remark 1: Note that (5) is not simply a more complicated Schrodinger equation as 
E s is on both sides of the equation. 

Remark 2: From §15.4.R4, the Schrodinger equation is recovered from O(^) 2 , so: 



P 2 



v\2 



( E .-^ = i^ + o(?) 4 



since 



a ■ p{E s - V)(a ■ p) 
4m 2 c 4 



-*(E, - V)x 



t>\4 



15.4. THE DIRAC EQUATION FOR THE HYDROGEN ATOM 



227 



So in the third term of the rhs of (5), it will be sufficient to let (E s — V)x 

-2 

2^X- But w e have 

(E s - V) *J? X 

awkward 



So 



(E s -V)&-px = &- [(E 8 - V)p X ] 



& ■ pEsX + & 



-Vp + pV - pV 



a ■ p(E s -v) x + o- 



p,V 



X 



which we substitute into (5): 



EsX 



2m 



+ (&-p)(&- 



p,V 



4m 2 c 4 



X 



n 2 n 4 P P> V 
— — \-V : 

2m 8m 3 c 2 Am 2 c 2 



4m 2 c 2 



X 



(6) 



Remark 3: (6) now looks like the Schrodinger equation! 



Remark 4'- On the rhs of (6), the 3rd term is the relativistic correction to the kinetic 
energy. 



Remark 5: The 5th term can be manipulated: 



\PJ(x)] 



ox Ox 



'df 

k OX 

ox 



ox ox 
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ia px p,V _^ 

4m 2 c 2 4m 2 



"7J(T • p x 



-iHV 



he 2 1 . A 
-- — 5-9 -oO" p x r 

fe 2 1 

r<r • r x p 



4m 2 c 2 r 3 



where it is ok to switch r and p since the p operator never operates on the conjugate 
coordinate. 



e 2 1 



S L 



2m 2 c 2 r 3 

■H5.0. — *• the spin-orbit interaction 



where it is ok to switch r and p since the p operator never operates on the 
conjugate coordinate. 



Remark 6: The 4rth term in (6) is not Hermitian. 



x| 2 d 3 ^7^ constant 



=^ X is not a good Schrodinger wave function 
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But what is required is 



^ipd 3 r = / [|x| 2 + M 2 ] d 3 r = constant 



, («■ • PX) t (o - • Px) 



lxl a + 
lxl 2 + 



4m 2 c 2 



d 3 r 



4m 2 c 2 



d 6 r 



^2 



X' 



X' 



X' 



1 + 



2m 2 c 2 

2 



X <i 3 r 



1 + 



1 + 



P 



8m 2 c 



2 

2 \ t 



P 2 



8m 2 c' 



8m 2 c 2 



1 + 



P 2 



8m 2 c 2 



X d 3 ^ 
Xd 3 r 



1 + 



8m 2 c 2 



X 



P 2 

8m 2 c 2 



X 



d 3 r 



J \Xs\ 2 d 3 r 



So to C(|) 4 , Xs nas a constant norm. But now we must write (6) in terms of 

Xs- 



1 + 



p 2 



8m 2 c 2 

Since E s \ = H\ by (6), we have that 



X 



E s 1 + 



P 



8m 2 c 2 



EsXs 



1 + 



1 + 



P 



8m 2 c 2 

,2 



ff 1 + 



P 



P 



8m 2 c 2 



8m 2 c 2 

2 



Xs 



P 



8m 2 c 2 



Xs 



H + 



H + 



HsXs 



1 



8m 2 c 2 
1 

8m 2 c 2 



p 2 H - ifp" 



Xs 



x* + 0(f 



P 2 ,H 



(7) 



(8) 
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H x 



H x 



H + 



H + 



H + 



1 



8m 2 c 2 
1 

8m 2 c 2 
1 

8m 2 c 2 



p 2 ,H 



(9) 



Remark 7: Equation (9) gives an extra term that comes from creating a suitable 
wavefunction with constant norm, so that (8) is a suitable Schrodinger equation 
that includes relativistic effects in the Dirac equation. Combining this extra 
term in (9) with the non-Hermitian 4rth term in (6) gives: 



1 



8m 2 c 2 



pp,V 



1 



4m 2 c 2 



p,V 



1 



8m 2 c 2 



V 



p,V 



+ 



p,V 



■p-2p 



p,V 



where we used [AB, C] = A[B, C] + [A, C]B. 



1 



8m 2 c 2 

h 2 

8m 2 c 2 
-e 2 h 2 



P 



p,V 



+ 



p,V 



V 



V 2 V 



-V - 

8m 2 c 2 V t 



2fc2 



e 2 K 
2m 2 c 2 
H D 



vr5 3 (r) 



which is the Darwin term! 



Remark 8: Putting these together (Relativistic Kinetic Energy, Spin-Orbit Interac- 
tion, and Darwin) gives the potentials that, in perturbation theory, resulted 
in the fine structure for the Hydrogen atom. 

Remark 9: Recall that the Darwin term contributed when I = and the S.O. 
coupling contributed only when I ^ such that the contribution gave the 
same result independent of I. This is in spite of the fact that Ho has nothing 
to do with angular momentum. Physically, this term comes because a particle 
can not be localized more that (the compton wavelength). 

V(r) — > V(r) = some smeared average about r 



Remark 10: Solving the full Dirac equation (rather than just to O(^) 4 ), gives the 
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energy spectrum: 



E nj = mc 



1 + 



a 



n 



U + + [U + h) 2 - « 2 



1/2 



1/2 



(11) 



Expanding in powers of a: 

• O(a ) — > rest energy 

• 0(a) — ► Schrodinger energy 

• 0(a 2 ) — > fine structure (compare with perturbation results) 



Remark 11: Equation (11) is still not 100% correct when compared to experiment 
(One can not account for Lamb shift for example). To get this correct, the 
E + M field itself must also be quantized (i.e. the field must be treated as a 
bunch of harmonic oscillators). 



15.5 Negative Energy Solutions 

In general 



t~\ i i • •! ,/". i. t\ yf i • the possible outcome of a system 

Relativity + Quantum Mechanics = p of 1 or several rticlc ^ 



particle production 
with enough energy 



arbitrarily enough 
energy for short time 



Question: How does Dirac theory account for this? 



Answer: With Negative Energy Solutions! 

For simplicity, let h = c = 1, then the free particle Dirac equation is: 





dip 






X 




JP. 





or 



X 



a ■ p + f3m ) i/j 



E — m —a ■ p 
—a ■ p E + m 



(E- m)x = & p<p 
(E + m)ip = a ■ px 



X 




"0" 


JP. 








E 2 — m 



V 



E 2 



> E 2 =p 2 + m 2 
E = ± \/p 2 + m 2 



(t) 
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Remark 1: We can not neglect the negative energy solutions, though note that if p = => E = 
±mc 2 . 

But what do we do with the negative energy solutions? 



Even if +mc 2 states are stable, quantum fluctuations should be enough to allow transitions to 
negative energy states. The problem is that this is not observed. 

15.5.1 Dirac Solution 

The vaccuum is really occupied by fermions, which don't allow positive energy particles 
to occupy negative energy states (Pauli-exclusion principle). Consider 



o 



AE > 2mc 2 



oooo#oooo 

ooooooooo 



Positive Energy Space 



Occupied by negative energy 
"particles with charge -e 



Remark 1: The filled Dirac sea is unobservable, but the hole is observable-it has 
change +e and \E\ > mc 2 , and is called a positron. 



Remark 2: When an electron jumps back into the hole left in the Dirac sea, 2 
particles are destroyed (an electron and a positron), but an energy > 2mc 2 is 
liberated in the form of a photon. 



Remark 3: Note that this explanation only works if the particles are fermions. But 
bosons can have negative energy solutions (and the negative energy solutions 
of the spinless Klein-Gordan equation are what motivated Dirac to look for a 
first order theory in the first place). 



15.5.2 Feynmann's Solution 

Negative energy solutions only travel backwards in time. 
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Remark 1: We still have the notion of an anti-particle (positron in this case). This 
is how things progress: 



t 



t d 



negative energy particle 



travels backwards in time 



To be destroyed 




What we observe is the positron (particles with positive energy) created at td and destroyed 
at t c . This implies that the creation of the positive energy positron coincides with the 
destruction of the negative energy electron (and vice- versa). 
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Chapter 16 



Feynmann Path Integral Formulation 
of Quantum Mechanics 

In the Schrodinger approach, we first find eigenvalues and eigenvectors of the Hamiltonian, then 
determine the propagator, U(t) (the full solution). In the path integral approach, U(t) is found 
directly. 

16.1 Path Integrals - Concepts and the Classical Action 

Here we consider only a single particle in 1-D (the generalization to higher dimensions is straight 
forward). 




(xo,t ) 



where the initial and final points are fixed. 
Recipe: 

1. Draw all the paths between (xo,to) and (x,t) 

2. Calculate the action S[x,t] for each path. 




235 



236 



CHAPTER 16. FEYNMANN PATH INTEGRAL 



3. U(x,t;x ,t ) =A^2 M 

paths e * j where A is a normalization factor. 
Remark 1: Recall that the action is a functional defined by 

S = [ C{t')dt' 

Jto 



(1) 



where C(t) is the Lagrangian (usually C = T — U). Also, a functional eats a function and 
excretes a number. 

Remark 2: Note that the propagator is not a weighted sum, thus, the classical action does not 
seem to be more likely (as it should in the classical limit). However, the paths are all treated 
as phases, which add together constructively near the path that corresponds to the classical 
action, and destructively for paths away from the classical action path by more than 5. 

Example: For a free particle with x = t: 




(0,0) 

Here the alternate path is x = t 2 . Letting (xo,to) = (0>0) an d (x,t) = (1 cm, 1 s): 



S d = / [^mi 2 ] dt 
Jo 

Salt = / 
J 



x=t=z-x=l 



/ dt = |m 
Jo 



T^mx 2 ] dt = 2m / t 2 dt = |m 
Jo 



For a classical particle, m ~ 1 gram, which implies that the phase changes by 



Salt — S c l 



m m 



1.6 x 10 26 » vr 



h V3 21 h 6h 

So we can ignore non-classical paths. For an electron, m ~ 10~ 27 grams, implying that 

Salt — S c i m 



and so we can not ignore non-classical paths. 



.16 < vr 
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16.2 U(t) for a free particle 

Recall that propagating U (t) as the solution to the time-dependent Schrodinger equation gave 

U(x,t;x ,t ) ^U(x,t-x ) = [^) 1/2 e m ^-^^ 2tu (1) 

It turns out that for potentials of the form 

V = a + bx + ex 2 + dx + exx 

=> U(t) = A(t)e lS -'/ h (2) 

with A(t) a normalization (HW). A free particle is of this form, so show that (2) is sufficient to 
obtain the exact answer (1). For a free particle: a = b = c = d = e = 0. So, 



S 



C(t) dt 



to 



to 



^mx 2 dt 



(3) 



(x N ,t N ) 



Imagine that the path is broken up into n-discrete steps so that as iV — ► oo, the path is continuous: 



X 




I I I 

to ti i'2 



Remark 1: In this picture, "draw all the paths" from (xo,to) to (xat^at) corresponds to using all 
possible values of x at a given t (like moving the "beads" up and down on an abacus). 
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Remark 2: "Calculate the action for all possible paths" is represented in the discretized form by 

N-l / N 2 

s= ^mfx 1±L -x 1 \ wit he = 8t (4) 

t=o ^ ' 

Remark 3: Finally, calculate the propagator as a "sum of all paths" . Symbolically, this is written 

iS[x(t)\/*T>[ Xt t] (5) 



as 

rx N 



e 

•I'D 



where 



In this case, we write 



rx N 

/ P[ai, i] =>• Integrate between all paths between xq,xn 

Jxq 



U(x N ,t N ;x ,t ) = / e^' h V[x{t)] 

Jx 

/°° f 00 f°° ( im (xj + i — Xi 

dx\ \ dx2 • • • / dxN-i exp — — > 
J-oo J-oo \2h £jj 



Scaling x by a/^t- 



2fe' 

/m\V2 , /2fe\^ (JV " 1} f°° , f°° r f°° , («i+i-»i) 2 

W = (— j = Jim (_ 1 / d yi / rfw - - - / ^_ ie - 



£^0 

Taking terms with yi implies 



-oo ■/ — oo 



dyie * e = y —e * 



Inserting this result into the next integral... 

dy 2 e~ te ~ y2)2/i e-^ 2 - y °) 2 / 2i = Y 2 e -fa3-w) 2 /3» 

Doing all the iV — 1 integrations gives: 

iV->oo \ m / iV 1 /^ 

£^0 



Jim 4 / Zl7rn£ \ N -l/2 e -(x N -x ) 2 m/2iNhe 

\ m J 



N— >oo 

£^0 



Jim A ( m ( 171 ) V2 e -{x N -x,f m ^KNe 

\ m J \2iirhNe) 



N^oo 

£^0 
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Remark 4 : Taking the limit N — > oo and e — > gives Ne = t, which is the exact answer in (1) if 

V m J 



A 

This means we can write 



Jv[x,t] = Urn £) 1 7 oo 7^/' "/ooT^W (6) 

7 e^O 7 00 l^FfeJ 7 1/00 l2i¥fej 

Remark 5: The propagator, U(xf,tf;xi,ti) is sometimes called the Feynmann amplitude. 



